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Equation | is a definition of the practical pH value 
and of the practical seale of pH. Evidently the 
function of the galvanic cell outlined above ts simply 


where is the gas constant, and 
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4. Calculation of pH be expected to display some curvature as the chloride 
molality decreased Furthermore, the different solu 
The pwHl was obtained from each emf value by bilitv of the tartrate in the three chloride solutions 


eq 2 These values at each tem wrature were plotted 
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The positions of the straight lines shown in figure 1 da 
were determined by the method of least squares 0.01 
For the other temperatures, the lines were drawn by MOLALITY OF NOC! 
inspection. The intercepts obtained by inspection 
at 25° agreed to 0.001 unit with those obtained by 
the method of least squares. The concentration of 
chloride in these mixtures was of the same order of 
magnitude as that of the tartrate, and the lines might 
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It is apparent from figure 2 that the pH of 
tions of potassium hydrogen tartrate increas 
only 0.09 unit when the saturated selution is dil 
to 0.01 m. Furthermore, it is ordinarily conver 
to prepare the unsaturated solutions on a vol 
rather than on a weight that is, to 
definite molarity instead of a definite molality 
should be emphasized that this change does 
involve an alteration of the standard state to wh 
the numerical the activity are referr 
compare footnote 3 Inasmuch as d(pH)/d 
known from the slopes of the plots shown in figure 2 


basis basis, 


values of 


only a determination of density is needed to reveal 
the difference of pH between 0.01-, 0.02 
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equal molarity. The density of the saturated solu 
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between the pH of the 0.03-molar and the 0.03-molal 
This conclusion applies as well to the 
more dilute solutions. Consequently, the pH, values 
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potassium were 
their pH 
meter calibrated 
phosphate buffer solutions 
at 24 
tained in an earlier 


values i vlass-electrode pH 


with the 


74 > =e 
Sb.00 Or 6.00 


at 25 table 3 

The change of pH suffered by 
ards on alteration of the temperature is primarily 
determined by the variation of —log A,A, for tartaric 
This quantity has been shown to pass through 
temperature near 35° [17]. Con- 
sequently, the change of pH, between 25° and 60 
than 0.03 unit, as figure 3 demonstrates 
The pH of the 0.01-M] solution is about 0.08 unit 
higher at 0 20 
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is less 


than at 25 

The growth of molds in solutions of potassium hy 
drogen tartrate, particularly in the saturated solu 
tion, sometimes causes appreciable changes of pH 
within a few days after the solution is prepared [22] 
It is easy to replace the saturated solution as fre- 
quently as but the solutions of lower 
concentration, though less susceptible to molding, 
are not as easily prepared. Hence, some experiments 
were performed to ascertain the usefulness of thymol 
A saturated solu- 


hecessary 


in retarding the growth of molds 
tion of potassium hydrogen tartrate to which a large 
crystal of thymol had been added was found to be 
free of mold 2 months after preparation, W hereas two 
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others without preservative molded in about 14 days 
Two 0.01-M solutions of the tartrate salt, to one of 
which had been added were both free of 
molds 6 weeks after preparation. Thymol is only 
slightly soluble in solutions 0.006 
mole dissolves in a liter of water at room tempera- 
ture), and pH measurements with the glass-electrode 
pH meter indicated that its presence did not alter 
the pH of either the 0.01-\/f or the saturated buffet 
solution by more than 0.01 unit. Experiments with 
the hydrogen electrode were unsuccessful; il appeared 
that dissolved thymol was catalytically reduced by 
hydrogen at the platinum surface Nevertheless, 
the usefulness of thymol in delaying the molding of 
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that a ervstal of thyvmol, 8 to 10 mm in diameter, 
be added to each bottle of standard tartrate buffer 

Potassium hydrogen tartrate is not appreciably 
is dificult to imagine that the 


hygroscopic, and it 
moisture to alter signifi 


salt could absorb enough 
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becoming noticeably wet. A 10-percent error in the 
concentration of the 0.02-.\7 solution changes the pH 
by only 0.01 unit Henee, it is usually unnecessary 
to drv the salt before it is used 

A good grade of distilled water should be used for 
the preparation of tartrate standards, but dissolved 
carbon dioxide need not be removed. The conduc 
tance of the water should not exceed 2 107° ohm 
em! at 25 Higher conductances may indicate the 
presence of small amounts of acidic or basic impuri- 
ties that have a relatively large effect on the pH. 
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The solubility of potassium hvdrogen turtrate is 


20°, 0.034 M at 25°, and 0.040 
reference to figure 2, it be 
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pH of that saturated at 25 The pH at 30° of the 
solution saturated at 30° also differs by about 0.01 
unit from the pH at 30° of the solution saturated 
at 25 Accordingly, the temperature at which the 
saturated standard ts prepared need not be carefully 
controlled It is suggested that the excess of the 
powdered salt be removed by filtration or decanta- 
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Heats of Formation of Some Barium Aluminates 
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l. Introduction 


The hvdrated barium aluminates have been studied 
by several mvestigators mm order to determine thei 
relationships [{1, 2, 3, 4],' but there 
no thermochemical data available on these 


solubility are 
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The purpose of the present investiga- 
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2. Experimental Procedure and Results 


2.1. Preparation of the Compounds 


barium 


he anhvdrous aluminates were 
obtamed by igniting appropriate mixtures of BaCO 
and ALOs3H,O (gibbsite) at 1,400° C to obtain 
S2QO. ALO. and at 1.300° C to obtam 3BaQO.ALO 
he hvdrates BaO ALO H ) and BaQ. ALO 2H ) 
were prepared hydrothermally; BaOQ.Al,O,.4H,O was 
prepared by allowing BaQO.ALO,.7H.O to stand in 
contact with an aqueous barium aluminate solution: 
7BaO.6ALOs.36H,0 and BaO.Al,O,.7H,O were pre 
pared by pre¢ ipitation from a supersaturated aqueous 
solution of Ba(OH and BaO.ALO,; and 
2BaQO.Al,O,.5H,0 was prepared by boiling gibbsite in 
a Ba(OH Full details of the preparation 
of these compounds hay elsewhere 
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The compounds were stored in glass jars for vary 
ing lengths of time prior to this investigation. The 
composition of each had been established by chemical 
analysis and petrographic examination, and it was 
assumed that no significant changes had occurred 
The uncertainty in the of the com- 
pounds was estimated to be of the order of 1 percent 


two 


solution 
bye en presented 


composition 


2.2. Heat of Solution of the Compounds 


The heats of solution of 2-gram samples of each 
compound in 640 grams of 2.00 N hydrochloric acid 
determined in an isothermally jacketed 
calorimeter [5]. Solution was complete except in 
the BaO.Al,O;.H,0 and BaQO.Al,O,.2H,O 
where about 1 and 0.1 percent, respectively, of in- 
soluble residue remained. An analysis of these res- 
idues indicated that they were uncombined gibbsite. 
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The initial 
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that the final temperature would be approximately 
25° ©, the temperature of the bath The heat-of 
data (see table 1) are the calculated iso 
thermal lhe of solution at 25 &. It be 
noted that the heats of solution per gram ol the two 
anhvdrous compounds are higher than that of any 
of the hvdrates, and that in the series Ba ALO HO 
BaQ. ALO, 2HLO saQO. ALO, 4H.O, and BaO ALO 
7H.O, the heat of solution decreases as the percent- 
age of combined water in the compound increases 
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2.3. Heats of Formation 
The heats of formation of these barium aluminates 
mav be calculated from their heats of solution and 
the heats of solution of BaO and ALO, in the acid 
used The heats of solution of these two oxides in 
2 00 N 
gibbsite dissolves very 
difficult to prepare 
Thorvaldsen, Brown, and Peake calculated the 
heat of solution of ALO HC] HO. Theu 
method is applied in the present work except that 
S| are used The heat of solu- 
HC] 26.61 H.O is obtained as 
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su Is 


were calculated 
at 25° C and 


acid 
slow ly 


hvdrochlori 


t) 
in 200 
more recent data {7 
tion of ALO, (e in 
follows 


a 


lir 


nas HCl, 2 
1eous H¢ 





2Al(e HCl, 


12.74H,0O 
All 


AIC] 


t 
; 5.67 keal 


i$ 21H.0 


45 t} 


ALO, (cla 


2Al( 


399.09 keal 


All 


20 Od 
P0495 ke al 


All 


SLO 


O71 


) 


HO 2(AlC] 


0.24 keal 


2AIC] 
HLO 


iS.21 


All 
The 


sum of the above four equations rives eq 


ALO Ca OOHCy 
9.71 H,O 


2(AIC] 


59.77 keal 


12.74 H,O 
All 


‘l he desired heat of solution however, Is re presented 


by the following equation 


ALO 


S| 


Hc} 
4 HOO 


~H 61 


All 


HO AIC |, 


59.95 keal 


ar 


Kquation 6 mav be obtained by adding to eq 5 the 
following two equations, which represent the neces 


sary dilutions 


2(AIC 1), 39.71 H,O ) 


All 


AIC] 


2 YO ke al 


S| 


4 HO 


GHC], 26.61 HO 


A//, 


HH! 


y 7s ke al 


12.74 H.O 
S 


The value of A//, was determined by the authors 
to be 0.24 keal: that of A// 59.77 keal. the 
sum of A//,, A//,, A//,, and A//, The value of A// 
was found to be —2.96 keal 

The heat-of-formation data for aqueous HCI [7] 
solutions were plotted against the square roots of 
their corresponding molalities, and from the resulting 
curve the value of A//, was estimated to be 2.78 
keal In each case, A// represents the heat effect of 
the equation as written 

The value of A//,, which is the sum of A// 
and A//, was calculated to be 59.95 keal 

The heat of solution of BaO in 2.00 N hvdrochlori« 
acid was obtained as follows 


Is 


All 


BaO( 2( HCI, 26.61 H,O 


HO A//, 


BaCl 


64.28 keal 


H.O 


BaCl,, 53.28 H.O+ H,OU BaCl, 54.28 H.O-+ A/T, 


10 
The sum of these equations ts 


BaOu 
BaC] 


2(HCI, 26.61 H,O 
54.23 H,O All 64.28 keal 1] 
which represents the solution of BaO in an exactly 
equivalent amount of 2.00 N hydrochloric acid. The 
value of A//, was obtained from heat-of-formation 
data available in the 7). The 


literature [7] value of 
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A//,) was assumed to be negligible; consequently 
value of A//,, is the same as that of A//,. 

When the equations describing the solutions of 
compounds are subtracted from the sum 
the solution of the 
equations are obtained from which the heat oly 
tion of the constituent oxides may be deduced 


ol 


equations desertbing 


OX 


bye 
culated from the heats of reaction of their constit 
oxides, using their known heats of formation 

The calculation of the heat of formation 
BaO. ALO, is presented in detail to illustrat 
method The three that represent 
solution of the constituent oxides and the compo 
are 


heats of formation of the compounds may 


equations 


BaOQ( 2(HCIL, 26.61H.0 


AH, 


BaCl 
64.28 keal 


54.23H.0 


AL,O 


bCAICI, 26.61H,0 
All, 59.95 keal 


2(AIC],, 81.3410 


taC'L, 2AICL. 150.3HCI, 4218SH.0 
ISS.3(HCL 26.61H.0) All 


BaQO. ALO 
100.49 keal 


The sum of these three equations 1s the follow 
equation 


BaQl 
1218 HO 
BaO. ALO 
BaC] 


ALOg(c.a BaClL 2AICL, 150.3HC] 
2(AICI,, 81.34H.O 

150.3(HCI1, 26.61H,0 

14.23H.0 All 23.74 keal 


It may be seen that eq 13 contains the desired eq 
tion, which 


Is 


BaO(e ALO 


A/1,, 


a sa0 AL.O4 (« 
23.99 keal 


‘ 
In order to obtain eq 4 from eq 13, the followin 
two equations may be added to eq 15 


150.3 HC, °6 61H () 
150.3HCI, 4054H.0 


BaC] 
All 


54.23H,O 
0.072 keal 


taC 


BaCl,, 150. 3HCI, 4054H,0+-2(AlC] 
BaCl. 2AICk. 150.3HCI, 4218H.O 
All 0.175 keal 


81.34H,.0 


To calculate A//,,, it was necessary to determi 
AH,, and Al/\, experimentally 

The sum of A/7,\,,—Al1,;, and — Al), is the desir 
value of A//,,. The values of A//, and Ald, ha 
been shown to be 59.95 and 64.28 keal, resp 
tively. The value of A//7,. was determined by 
authors to be 100.49 keal. The value of AH/ 
sum of A//,, Al7,,, and Aly», was then calculates 
be 23.74 keal. To this was added the values 
AH, and Ald, determined by the authors to 
+-O0.072 and 0.175 keal, respectively, to obt 

23.99 keal, the value of A/7,,. 

In an exactly analogous manner, the heat of 














mation of 3BaQ.ALO 


rr presented by eg 17, may be calculated. 


3BaO. ALOg( 
14.05 keal (17 


ALO, (c,a 
AH 


SBaQ( 


The procedure for calculating the heats of forma- 
tion of the hydrates is the same, except that the water 
The calculations 


presented in detail to illus- 


of hydration must be considered 
for BaO.ALO,.H.O are 
the method. The equations representing the 
solution in 2.00 N hydrochloric acid of the constit 
rent oxides and of the compound are 


trate 


BaOQ( 2(HCIL, 26.61H,.0 BaClL, 54.23H.0 
AH 59.95 keal 6 
ALO, (e,a) + 60HCL, 26.61H,O) = 2(AICI,, 81.34H,0 
All 64.28 keal 1] 
BaC lL, 2AICL, 161.5HC1, 4517H,0 
BaQ. ALO, HLOCe 169.5(H Cl, 26.61H.O 
All 91.02 keal IS 
H.O+ BaClL, 2AICL, 161.5HClL, 4516H.0 
BaCl, 2AICL, 16L.5HCI1, 4517H,0+ Al/ 19 
The sum of these four equations may be separated 


mite the followimeg two equations 


BaQ ALO HOG BaQO. ALO. HO 
All 32.95 keal 20 
BaC'L, 2ZAICL, 161.5HCI1, 4516HLO saC] 
54.231H1,.0+-2(AIC],, 81.34H,0 l6L.SCHCI 
°26.611LO All 0.262 keal 2] 


where eq 20 is the desired equation and eq 21 may be 


broken down into components whose heat effects 
mav be separately determined as follows 
161.5 (HCI, 26.61H,O BaClL, 54.23H.0— BaCl 
161. SHC, 4298H.O All 0.074 keal (22 
taC | 161. 5HCI, 4298H.0+4-2(AlC] 
81.34H.0 BaCl, 2ZAICL, 161.5HCIL 4516H.LO 
All O.1SS keal 23 
The values of A//, and A//,, have been shown to 
hye 59.95 and 64.28 keal, respectively, and the 


value of Al. 

91.02 keal. Anattempt was made to determine the 
value of A//\y experimentally 
so small that it 
degree of certainty with the 


was determined experimentally to be 


but its magnitude was 
measured any 
available apparatus and 
The sum of A//, 
calculated to be 


and All 


could not be with 
is therefore assumed to be 
All All, and Al/ then 
34.89 The All 
determined by the authors to be 0.074 and 
Adding these to the 
32.95 keal for A// 
the formation 
oxides of the other barium 


Z0TO 
was 

were 
O.1SS 
previous 


keal values of 
keal, respectively 
sum gives a value of 
In an analogous mannet heats of 
from their constituent 


OTOSR22 





from its constituent oxides, as 





aluminate hydrates may be calculated The corre 


sponding reactions are 


BaQ(e ALO, (c,a PHO’ BaQ) ALO, 2ELOK 
All., 16.39 keal 24 
BaOl ALO, (c,a HHLLOU BeO. ALO, 4ELOt 
All 15.36 keal 25) 
BaOQ( ALO, (ea) + 7HLOW) = BaO. ALO,.7 HOt 
All 51.29 keal 26 
2BaOQ( AL.O,(e,@) 4+-5H.O(I 
2? BaO. ALO, SHO All TS.SS keal 27 
7BaO( GALO,(e.a SOHO’ 
7 BaO.GALO,.36H.0 ( All 274.60 keal (28 


For ere hy of thes compounds 
Was written, and 
were performed 


tional equations of the type already shown 


a group of equations 
appropriate dilution 
This work required some 17 addi 
These 
equations differed slightly from those already given 
the heats of were all determined 
using 2-gram samples and a fixed 40 
of 2.00 N HC] It was felt that these equations and 
calculations shown in detail lt 
found to be to perform each dilution 
f the heat effects 
Accordingly, the values table 2 


plotted and the desired heats of dilution determined 
} 


experiments 


be« nuse olution 


amount 
need not be Wiis 
UNNECeSsSaLy 
rane ol the 


because 0 narcow 


presented in were 


Vv inte 1polation 


From the known heats of formation of the constit 
oxides, the heats of formation from their con 
stituent elements of thes 
calculated and are shown in table 1 

The authors’ value of the heat of formation of 


BaO. ALO, differs from that of Avgustinik 


uent 


barium aluminates were 


widely 
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and Nichedlov Petrossian [11] who calculated then 
value from the pressure developed by heating un 
equimolar mixture of BaSO, and ALO These 


mvestivators however. did not identify thre products 
of ther reaction, and in a ibsequent investigation 
[12] in which one of them analyzed the products 
formed by igniting BaSO, and ALO, in the appro 
priate stoichiometri ratio, he was unable to produce 
BaQ ALO 

Grube and Heimtz [13] have ealeulated the heat of 
formation of 2BaQO.ALO, from BaQO ALO, and BaO 
to be > li keal/mok The existence of 2BaQO. ALO 


howe ver, ms hot accepted by many wmvestigators and 


the present authors made no measurements with 
amples of that stoichiometric composition 
3. Summary 
Thy heats ol solution m 2.00 N hvdrochlori acid 


thy BaQ-Al,O,-H,O 
were determined in an isothermally jacketed calorim 
The of 
compounds were found to be higher than those of th 


of eight compounds u svstem 


eter heats solution of two anhvdrous 
a series of hydrates that differ only in 
of ol the heat of 
decreases as the percentage of combined 
The heats of formation of 
barium aluminates were calculated from their heats 
the heats of solution of BaO and ALO, in 


acid used in the litera 


hvadrates 
the 
solution 


amount water hvdration 


water mere these 


uss 


of solution 


the and from data available 














Lure From the data presented the heats of h 
tion of several barium aluminate hvdrates ma 
readily calculated 
The authors gratefully ne knowledge the assist 
of EK. T. Carlson, who prepared the compo 
studied in this investigation 
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Thickness of Glass Electrodes 
J. J. Diamond and Donald Hubbard 
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l. Introduction 


The effeet of the thickness of the membrane of a 
lass electrode on its electrode function ts of interest 
m connection with a continuing study of the class 
electrode [10, 12, 13, 16] Kratz IS] made the state- 
ment that the theoretical variation of potential with 
pH can be obtained only with thin-walled membranes 
Kahler and DeEds |17] found that there was roughly 
a linear relationship between the “average electrical 
thickness” of an electrode and the difference, in milli- 
volts per pH unit, between its electrode response and 
This difference up 
proached zero as the membrane thickness approached 
zero Dole [7], however, quotes von Steiger [20] to 
the effect that glass electrode potential is completely 
The pres 
ent study has been made in order to clarify the 


that of a hvdrogen electrode 


independent of the thickness of the glass 


matter 

NMacInnes ana Dole 1} determined the thickness 
of their thin membranes by weighing a measured area 
of uniform thickness, this uniformity being deter 
mined by the use of interference colors. Knowing 
the vlass density, they calculated then films to be 
about lu thick. Von Steiger [20] measured the thick- 
ness of his Cremer-Haber electrodes [5, 9] by bringing 
a glowing glass rod up to a selected distance from the 
glass surface. By observing the reflection against a 
black background, he observed interference fringes 
and then calculated thicknesses. Yoshimura [22], in 
his study of the effect of thickness on asymmetry 
potential, assumed the thicknesses of his electrodes 
were in the order of their electrical resistances and 
had no direct knowledge of the thicknesses involved 
Kahler and DeEds [17] measured the thicknesses of 
their electrodes by using a microscope with a cali- 
brated fine focus adjustment to focus on the two 
surfaces of the membrane. The difference in setting 
corrected for the refractive index of the glass, gave 
the thickness 

For the purpose of this study, a nondestructive, 
magnetic method for the measurement of the glass 
membrane thickness has been developed. It involves 
the use of the Magne-Gage? [1], developed by 
Brenner [2, 3], for the measurement of electroplated 
coatings on metals, and is also used for measure 


F igure n brackets indicate the terature reference at the er vf th paper 
Manufactured by the Ame in Instrument ¢ Silver Spring, Md 


been developed for the measurement of 

the use of the Brenner Magne-Ciage as the 
ler as the backing matena Keach of the four 
acteristh departure thickness ranging fron 
rhitnimum thiekrne Was equal to or le ttiat 
1 50 millivolts per pH unit electrode response 
d increasing departures from this theoretica 
lasses investigated were found to be a funetiorn 


ments on paint, porcelain enamel, plastic, and other 
films or coatings that are or can be backed by iron 
or steel plates 

The adaptation made by us is the use of carbonyl 
iron powder as the backing material and the use of a 
corresponding direct calibration procedure 


2. Experimental Work 
2.1. Method of Measuring Thickness 


The operation of the Magne-Gage is based on the 
decrease in magnetic attraction resulting from the 
interposition of a nonmagnetic material between a 
magnet and a magnetic material, The apparatus, 
shown in figure 1, consists of a lever arm that is 
connected to a coiled spring, which is in turn con 
Thee ted to a cial and knob Suspended from the na 
of the lever arm is one of several interchangeable 
permanent bar magnets of different weights and pole 
strengths. Two of the magnets available have been 
found to cover the range of thicknesses of interest in 
this study 

In use, the steel-backed specimen ts pla ed in post 
tion under the magnet, and the knob is rotated 
counterclockwise until the presser arm, which is slip 
mounted on the assembly, brings the rounded tip of 
the magnet into contact with the surface of the speci 
men. The counterclockwise rotation is continued 
until the built-in stop is reached, or until a trial o1 
previous experience indicates that the tension on the 
spring has been released sufficiently so that the 
magnet will remain in contact with the specimen 
when the kneb is rotated clockwise, taking the presse1 
arm off the lever arm and magnet. The knob is then 
rotated clockwise, increasing the tension on the 
spring until the force is sufficient to pull the magnet 
away from the specimen surface. The dial reading 
at this point is a function of the thickness of the 
nonmagnetic material When the instrument ts 
properly calibrated with material of known thick 
ness, the dial reading can be converted to a thickness 
reading 

A Cremer-Haber type glass electrode is typically a 
thin glass bulb about 10 to 20 mm in diameter blown 
on one end of a glass tube a few millimeters in diam 
eter Since it is impossible to use a steel plat as a 
backing material, for a thickness measurement on an 
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Fieure | 

\ 
object of this shape carbonyl iron powder Is used 
instead The powder Is pla ed in the clean, drv ele« 
trode by a standard procedure, and its open end is 
stoppered After adjusting the magnet to the ap 
propriate height by means of the rack and pinion, the 
bulb is placed mn with the 


around the magnet, and the measurement 


contact plastic guard 
is made 
Contact with the guard serves to steady the specimen 
during measurement, and, if maintained during mul 
tiple measurements made to obtain a precise value 
for the dial reading, will insure that the readings are 
being made on the same spot. This is important, as 
all such glass bulbs are quite nonuniform in thickness 
One additional pre aution must be obs rved in han 
dling the instrument when making measurements on 
powder-backed, thin membranes It is easv for the 
magnet tip to dent soft or thin metal or to puncture 
a thin glass membrane unless contact is made gently 
and no pressure is exerted on the specimen by further 
rotation of the instrument knob 
But in the usual case, the spring tension must be de 
creased to insure that the magnet will remain in con- 
tact with the specimen when the presser arm ts lifted 
This is accomplished before contact is made and 
without removing the specimen from contact with 
the guard, by rotating 


counter lockw Ine 


the knob « ounterclockw ise, 


while using the thumb to prevent the presser arm 
from also rotating 


When used for the measurement of steel-backed 


THICKNESS , MICRONS 


l ] 
100 = ss 110 i20 13s”o0 
SCALE READING 








films o1 the Magne-Gage is calibrate: 
against thickness standards supplied by the mat 


coatings 
facture These standards are squares of steel eles 
troplated with known thicknesses of a nonmagneti 
metal 
for interpreting dial readings obtained with powder 
backed samples. We use instead thickness standards 
made by cementing aluminium or platinum foils of 
known and uniform thickness * over the open ends 
of glass tubes of about 9-mm diameter 


The resulting calibration curves are useless 


These simu 
lated electrodes are filled with carbonyl iron and 
measured in the same manner as the electrodes, The 
calibration curves resulting for magnets 2 and 3 are 
shown in figures 2 and 3 

In his work on the thickness of electroplated met 
als, Brenner [2,3] found that the geometry of t! 
specimen may have an appreciable effect on the seal 
reading obtained. In order to prove that no serious 
errors are introduced by the difference in shape |! 
tween the Cremer-Haber electrode and the straig! 
tube standards, the following experiment was p 


formed. A thick-walled bulb 15.5 mm in outs 
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diameter was blown at the end of a piece of 6-mm 
tubing \ second similar electrode Was prepared 
with a 19.0-mm bulb. The ends of the bulbs were 
ground off to leave holes about 2 to 3 mm in diametet 
Pieces of aluminum foil 118 « in thickness were then 
cemented over these holes These simulated elec 
trodes were filled with powder, the foil thicknesses 
measured with the Magne-Gage, and the readings 
compared with that obtained from the 118-u straight 
tube standard Seale readings of 141, 141, and 142 
were obtained, the higher reading for the 15.5-mm 
bulb. This is of the order of the experimental errors 
involved 

From the very nature of a powder, it seemed likely 
that the dial reading obtained on the Magne Cage 
would depend on the degree of compacting of the 
powder, as well as on its particle size. Experiments 
were accordingly made on four of the five sizes 
available. They were performed on the 118- 
straight-tube thickness standard and essentially 
duplicated on a typi al glass electrode The pow de | 


Available fron Ant Diy Gg 4 i: ( 


Tante | Vagqne-Gage ¢ eading Nected | ribration 


was packed by holding the filled electrode m= the 
normal bulb-down position and running a mechan- 
ical vibrator up and down the outside of the glass 
tube. In these experiments, a Sargent Ideal electric 
marker, set for a fairly high vibration amplitude, 
was used. Preliminary experiments indicated that 
vibration time required to reach equilibrium and 
Magne Crave scale reading nt equilibrium both de 
No quantitative 
data were obtained on this point, but it must be 
emphasized that the method and energy of vibration 
selected must be the same during measurement as 


pended on the energy of vibration 


during standardization 
From the data summarized in table 1, several 


facts are evident For each of the four grades of 
carbonyl iron powder tried, vibration will result in 
packing to an equilibrium condition. The time 


needed to reach equilibrium decreases with increas 
ing particle size. Figure 4, from the data for the 
straight-tube standard in table 1, illustrates this 
tendency 
librium decrease numerically with mereasing particle 
size, as illustrated in figure 5 As can be seen from 
the calibration curves shown in figures 2 and 3, the 


In addition, the seale readings at equi 


sensitivity of measurement, as microns per scale 
division decreases with MnerTeasiny se ale reading 
Partially to keep this sensitivity as high as possible 
for any given thickness, but mainly to take advantage 
of the saving in vibration time involved, we selected 
Measure 


ments were made after vibrating for 1.5 minutes, 


20-u iron powder (grade L) as standard 


The data table 2, indicate the precision possible 
\leasurements were made on two electrodes. using 
magnet 2, filling with iron powder L, and vibrating 
for 1.5 minutes. After every reading, the electrode 
was held bulb up, still stoppered, vibrated for about 
30 seconds until the powder separated from the 
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clectrode tip and then held bulb down and repacl 
for 1.5 minutes before the next reading was taken 
With ordinary care, scale readings can easily 
reproduced to within one unit From an exami! 
tion of the calibration curves shown in figures 2 anc 
it can be seen that one seale unit amounts to abo 
2 to 3 percent of the thickness measured, for thir 
nesses down to about 30 u For thinner membran: 
this percentage mcreases somewhat but measu 
ment to | wis easily obtained 


2.2. Measurement of Electrode Response 


bene h ele trode was prepared for voltage respon 
measurement by soaking in distilled water ov 
night or longet It was then filled with mereu 
21], into which the measurement lead was dippe« 
\ Beckman model G pH meter was used to measur 
the differences in potential developed in a series of 
Britton universal buffer solutions [4] between th 
experimental electrode and a well-conditioned com 
mercial glass electrode (Beckman 1190) serving as a 
reference electrode The difference between the 
measured potentials for the buffer solutions of pH 4 
and &, divided by four, gave the departure of the 
electrode from the theoretical, in millivolts per pH 
unit. This departure subtracted from the theoretical 


H EQUILIBRIUM MIN 


59 mv per pH unit, gave the electrode response of 
the experimental electrode 


3. Results 
3.1. Thickness Profile of a Glass Electrode 
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It is obvious from the nature of the glass blowing 


° 7 Se. SOR. process that a bulb blown at the end of a glass tube 
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ments of the thickness of glass electrode bulbs 
the powder backed magnetic method, 
that only in verv unusual cases will the thin spot be 
at the tip of the bulb in 
Almost 


using 
have shown 
where it might be expected 
line with the axis of the tube 
itis Loto (0 degrees off this spot 
In order to fully 
ass electrode bulb 
Beckman 1190) was emptied, cleaned and 
and its thickness profile taken by means of 
25 points on the lower half of the 
For presentation of the data, the bulb was 


invariably, 


the range of thickness 
a typical commercial elec 


illustrat 
nag 
trode 
dried 
measurements al 
bulb 
treated as a geographir vlobe with the tip as a pole 
and thickness measurements were taken along the 
latitude longitude The great circle 
passing through the thin spot of the electrode was 
The data are plotted 
in figure 6 on an azimuthal equidistant projection [6] 
of the bulb, where the 
pont 


circles ol and 


selected as the prime meridian 


tip is taken as the central 


3.2. Electrode Response versus Thickness 


bulb 
glasses, 
They 
9-mm, thin-walled Corning 015 glass electrode 
recently purchased this investigation 

2) Corning 015 glass slab purchased about 20 years 
Kimble soft glass tubing, and (4) glass A, 
length of 9 thick-walled tubing in- 
cluded in the shipment of Corning 015 tubing, but 
differing from it in appearance and in all properties 


Che dependence ol electrode on 


thickness was investigated for 


response 
four different 
two of which were purported to be identical 
are (1 
tubing for 


a single mm 


investigated 

The electrode bulbs tubes of the 
same glass for all except the Corning 015 slab glass 
Bulbs of that glass were blown on the thin-walled 
015 tubing, the only available glass tubing from which 
it would not crack away on cooling 

The electrode for bulbs of 
various thicknesses of the four glasses are shown in 


were blow nh on 


responses obtained 


o- —~ 
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that of the thin 
Phe area of this thin spot is 


figure 7. The thickness 
of the glass bulb 
estimated to vary from perhaps 2 to 10 


vive Is 
spot 
mm *, and 
would in any case be no smaller than the 0.75 mm, 
which MaclInnes and Dole found to be large 
full with 
membrane electrodes [19] 

The data show that, for each of 
investigated, relatively thin” bulbs 
theoretical 59 mv, pH electrode response 


enough 
ther thin 


to give electrode response 


the four 


rive 


glasses 
the full 
Kor each 
glass there is a characteristic “departure thickness’ 
thickness that will permit full theo 
Beyond this point, increasingly 
increasing voltage 
difficulties and 
measurement Krom the 
evident why the two batches of Corning 015 glass 
considered 2S itl rent This 
the earlier findings of others [8] that this glass 
may be somewhat 

When the hygroscopicities of the four glasses were 
compared by the method of Hubbard [11, 16], it 
was found that the higher the hygroscopicity of the 
This 
is shown in figure 8, and is consistent with previous 
full electrode 
hygroscopicity 


the maximum 
retical 
thick 

well +t hs 


response 
bulbs show departures as 
Increasing uncertainties of 


voltage data it is alsa 


were glasses bears 
out 


variable in composition 


glass, the greater was its departure thickness 


findings that it is impossible to get 
with very low 
It is also of interest to note that the loga- 
linear 


response glasses of 
[14, 15] 


rithm of the hygroscopicity appears to be a 
function of the departure thickness 


4. Conclusions and Summary 


1. A magnetic method, employing iron powder, 
has been developed for the simple rapid, and accu 
rate determination of the thickness of glass electrode 
The procedure has the additional advantage 
destroying the The method as 
given could be used with equal facility for thickness 


bulbs 
of not specimen 
measurements on odd shaped envelopes made of any 
nonmagnetic material 
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2. A Cremer-Haber electrode bulb is nonuniform 
in thickness and has a thin spot that is usually 15 
to 30 degrees from the tip of the bulb. The tip ts 
usually not, therefore, the weakest spot on the bulb, 


as has been supposed 
}. Thin glass electrodes made of any of the glasses 


mvestigated will vive full theoretical electrode 


response 


} Kor each lass investival« ad there 
\ glass electrode whose 


is a character- 


istic departure thickness 
minimum thickness is equal to or less than this de- 


parture thickness will give full theoretical electrode 


response, whereas a thicker one will not 

5. For the classes investigated the departure 
thickness ts a function of the hygroscopicity Those 
of vreater hygroscopicity will therefore permit the 
manufacture of electrodes of thicker and hence more 


rugged construction 
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Some Characteristics of Stedman Packing in the 
Distillation of Hydrogen and Its Isotopes 


Abraham Fookson, Philip Pomerantz, and Simon Rothberg 


(\pparatus was designed and constructed 


packing in the distillation of hvdroget 


and the hold Ip at 
charge Mixtures of 
distilled at boil-up rates of 454 milli 
From the 

found to be 


as sti hydrogen 


listillation data the height 
1.0 inch for the 12-ineh 


1. Introduction 
This investigation was undertaken to determine 
some of the characteristics of Stedman [7] ' packing 
of value in engineering calculations relating to the 
distillation of hydrogen tsotopes. Specifically, the 
still hold-up at various boil-up rates and the efficiency 
of the packing in fractionating hyvdrogen—hydrogen- 
deuteride mixtures were studied, and the flood point 
determined. To obtain these data, 
were made with a Stedman still that 
for use at very low temperatures 


measurements 
was designed 


2. Apparatus and Material 


The apparatus, similar to that described previously 
1], is shown assembled in figure 1 It consisted of a 
still, A, with a 25- by 300-mm stainless-steel Stedman- 
packed section and a graduated pot of 50-ml capacity, 
shown in detail in figure 2. A eylinder, integral with 
the still, surrounded the condenser and served as a 
reservoir for the liquid hvdrogen coolant The 
vacuum jacket of this cylinderical portion was con- 
tinuous with that of the still itself. The still was 
surrounded by a dewar (not shown) filled with liquid 
nitrogen and was manometer, B, 
enabling the pressure in the still to be read. When 
the still could be operated 
under total reflux conditions, when C was open, ad- 
withdrawal of distillate 
A regular valve, D, operated by solenoid 
kK, permitted control of the rate of distillation. This 
regulator consisted of a length of small diameter 
drill rod fixed to an iron The drill rod fit 
loosely into a glass capillary, the length of the rod 
in the capillary being controllable by the solenoid, 
the annular space through which the gas had to pass 
was thereby variabli 

The regulator valve led to a graduated Toepler 
pump, F, of 500-ml capacity, and to a bypass, G 
The latter was used when charging or evacuating 
the still. The top of the Toepler pump was con- 
nected to a manometer. The readings on manom- 
eters B and H indicated the pressure differential 
across the regulator, D. A length of pressure rubber 
tubing, J, led from the Toepler pump to a mercury 


connected to a 


Was ¢ losed 


stopcock }. 


mission of charge o1 was 


possible 


core 


Figures in bracket licate the literature references at the end t} 


isotopes were 
various boil-up rates were measured using 
hvdrogen deuteride of know: 
iters per hour and | 
equivalent 
packing used in these 


characteristics of Stedmat 


average sti hold 


h sore 
measured Phe i} 
both hvdrogen and deuteriutr 
COMpoOsitlor 
10O mi 
to a theoretical plats 
} 


liters per hour respect 
Was CAICHUIATE* 


experiments 


Raising or lowering of this reservoir pet 
adjustment of the mercury level in’ the 
pump and consequently of the pressure 

This. together with the regulator valve 
constituted the method of distillation rate 
control used in this work 

The Toepler pump led to a manifold to 
were affixed the bulbs for receiving samples (about 
20 ml each), and the 5,000-ml flasks 
used for storing hydrogen deuteride and collecting 
distillate The manifold was in turn connected to a 
mercury diffusion pump. A line (not shown in fig 
1) between the Toepler pump and the diffusion pump 
bypassed the manifold, enabling the evacuation of 
the system whether or not bulbs were affixed to the 
manifold 

The 99.979 H this 
work was identical with that used in the preparation 
of liquid hydrogenat the Bureau. It was prepared at 
the low temperature laboratory of the Bureau by 


reservoll 
mitted 
Toepler 
within if 
setting 


which 


calibrated 


hvdrogen vas used in 


These rec er \ sibbratedd wit! 


Volumet Cilassware t the Bure 


TO MANOMETER 6 








MANIFOLD - 
13 SAMPLE BULBS (20 CC) & 
@ COLLECTING FLASKS (5000CC) TO MANOMETER © 
ON MANIFOLD 


MERCURY RESERVOIR 
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» through No 
n Stedman 


the electrolysis of potassium hydroxide solution ay 
pumped into a steel evlinder of convenient size 

The hydrogen deuteride was prepared by 
method of Wender, Friedel, and Orehin [2], involy 
the reaction at O° C between lithium alumin 
hvdride and deuterium oxide For this work, abo 
50 liters ol hvdrogen deuteride was prepared, Co 
sisting of about 98.5 percent of HD, 1 percent of HH 
and O.5 percent of DD This material was us 
without purification, allowance being made for t! 
hvdrogen content when mixtures to be distilled we 
prepared 


3. Technique 


3.1. Hold-up Measurements 


The apparatus was cooled by filling the appr: 
priate Dewars with liquid hydrogen and nitroge: 
\ quantity of hydrogen or deuterium, about 50 1 
liquid as measured in the still pot, Was then intro 
duced into the still as ras, and permitted to liquef 
stopcock (" (fig. 1) was closed, and the exact volun 
noted. Current was then passed through the heater 
and the volume in the still-pot and the pressurs 
the still were noted when equilibrium had been 
established The difference between the initial 
volume and the equilibrium volume was taken 
the hold-up corresponding to the boil-up rate as 
calculated from the heat input. This was repeated 
for various values of the current 

The results of these experiments had to be correet« a 
for three factors |. There existed a heat leak into 
the still from the exterior, whose magnitude was 
estimated. 2. The initial volume as measured by 
means of the graduations on the still-pot was less 
than the amount of material actually in the still 
since an unknown amount was present wetting the 
column (static hold-up 3. At each boil-up rate 
the volume as read on the still-pot graduations at 
equilibrium did not take into account the volume 
of the bubbles in the liquid 


3.2. Heat Leak and Static Hold-up 


The magnitude of the heat leak was computed 
from the results of two subsequent experiments 
The first experiment was performed by charging 
hvdrogen into the still and reducing the pressure 
by means of the vacuum pump so that the liquid 
would boil at a temperature below that of th 
condenser. The vacuum pump was isolated from 
the system by closing stopcocks G and L, and 
hvdrogen gas was withdrawn with the Toepler pump 
at such a rate that the still pressure remained 
constant. The time rate of withdrawal of gas wa 
observed, and from the known heat of vaporization 
of hydrogen, the heat leak was calculated. Thi 
method will be referred to as the “withdrawal 
method” in later discussion 

The second method employed for estimating th 
heat leak consisted of introducing into the still 


All analyses of gas mixtures reported in this work were made by the M 
Spectrometry Laboratory of the NBS 


450 





hvdrogen, reducing the 


and observing the time necessary 


quantity ol pressure to a 


re asured low value 


for the pressure to rise tots normal equilibrium value 


From the known heat capacity of the still contents 
the still itself having negligible heat 
20° K the heat necessary to produce 
rise was calculated This method 
to the “time method 

second that of the 
hold-up 
quantity of 
still 
appeared in the pot 


capacity at 
this pressure 


W“ ill be referred 


factor “zero port 


volume. or “statue was computed by 


introducing a measured hvdrogen 


into the empty and observing the 


that The 


represented the material that 


deuteride 
volume of liquid 
difference in 
wet the packing plus the amount of reflux occasioned 
by the heat leak into the pot 

Phe third factor, the augmentation in liquid volume 
because of the gas bubbles, could not be 
had to bye 


volume 


in the pot 
corrected 
higher boil-up rates 


and volumes estimated at the 


3.3. Distillation 


this work 
shown in 


were 


distillations conducted in were 


with the 


The 
enrried out 
After the 
nitrogen 
deuteride 
flasks on the 


hgure | 
with 
charge of 


apparatus 
appropriate vessels cooled 
liquid hvdrogen il 
admitted to the still from 
the storage manifold. The Toepler 
pump Was used to withdraw the hvdrogen deuteride 
after the equilibrium 
150 mm had been reached From the initial 


their volumes 


liquid and 


hvdrogen was 


from the flasks pressure of 
about 
and final pressures in the flasks, and 
the quantity of hydrogen deuteride in the still was 
calculated. The volume in the still-pot was then 
and hvdrogen rus Was introduced from a tank 


until the volume was increased to a predetermined 


noted 


amount. The assumption was made here that the 
hold upin the column cal ulated volume of hvdrogen 
deuteride introduced minus volume neasured = in 
still-pot) would remain approximately 
and hence the volume 
hvdrogen introduced 

Current was then passed through the still heater 


while 


unchanged 


mMcrease in represented 


predetermined boil-up rate 
The 
manometer B was noted from time to time 
this was constant, distillation was started 
regulator 
Toeplet pump allowed to fill up, by 
lowering the mercury 
correspond to a predetermined reflux ratio Each 
time the Toepler pump was full (500 ml), stopeock K 
was closed and L opened, thereby discharging the 
either the vacuum pump, a sample 
bulb, or a collecting flask Samples were colle ted 
at intervals suitably spaced during the distillation 
so as to give a representative curve when the com 
amount 
Toeplet 


to provide a 


stopcock (" was closed pressure as read on 
and when 
Stopeock 


the 
continuously 


(" was opened LD) was adjusted ana 


reservoir, at such a rate as to 


eontents to 


position of distillate was plotted against 
distilled as volumes in the 
pump and pressures read on manometer H (fig. 1 
All gas volumes were corrected to 0° C and 760-mm 
Analyses on the were 
determined with the mass spectrograph 


measured by 


pressure several samples 


451 


The still pressure Was noted during the cours of 
the distillation 
value 


and when it had fallen to a constant 
it was presumed that the material in the still 
then pure hvdrogen Ihis 
corresponded to the temperature of 
at the boil up rate 


Wis deutericde pressure 
condensation 
hvdrogen deutericde used im the 
function of the rate of 


Was hot the same from run 


since it was a 
to the still, it 


experiment 
heat input 
to run 


4. Results 
4.1. Hold-up Measurements 


The data obtained in the first experiment relating 
boil-up rate with hold-up hvdrogen as. still 
table | The temperature, 7, 
in the still was caleulated from the obser ved pressure 


using 
charge, are given in 


P, by means of the following equation [3] 


i L() 020537 7 


Log P?O H 166687 


1) 


From this temperature 
zation, L 


the molar heat of vapori 
was computed, using the following equa 
tion 4 


H 219.7 —0.27(7 


16.6 
The molar volume ol the liquid yY mil mole) 
at the appropriate temperature was obtaimed by use 
of the following equation [8] 

H 24.747 


0.08005 7°) 0.012716 7" 


Krom the values of ZL, (cal/mole), V,, (ml/mole), the 
heater current, J amp and the heater 
of 46 ohms, the boil-up rate of the liquid, B 
was calculated using the equation 


resistance 
ml/h 

derived from the 
definitions of the quantities cone erned 


96 rot 


The factor 0 10° imeludes the heater 


resistance value and the necessary conversion factors 


constant 


Flooding was first observed, using n-deuterium as 


th hydrogen 





at a boil-up rate of 3,745 ml/hr, at which 


This 


still charge 
rate the pressure in the still was 1,990 mm Hg 


corresponded to a Vapor velo itv of 0.26 ft/see In 
a& previous experiment using hydrogen as still charge 


induced to flood within the 
The controlling factor here 
whose upper limit was 


the column could not be 
limits of the apparatus 
was the manometer, B (fig. | 
about 2,000 mm Hg 

In the ease of deuterium 
creased from the lower values, the level of liquid at 
the edges of the plates was substantially constant, 
increasing from approximately | mm at 2,626 ml/hr 
to about 3 mm at 3.420 mi/hi But as the flood 
point was neared the amount of liquid continued to 
increase in the top third of the column. Shortly, on 
the top plate liquid completely covered the cone 
and continued to increase until a column of liquid 
started to enter the This boil up rate 
was designated as the flood point 

The caleulation of the boil-up rates when deute 
rium was used in the hold-up measurements are not 
as accurate as were those when hvdrogen was used 


as the boil-up rate in- 


condense! 


This is a consequence of the fact that there does not 
exist a reported heat of vaporization-temperature 
relationship for deuterium as there does for hvdro- 
Hence, the deuterium analogues of eq 2 and 3 
could not be emploved. An approximation was 
obtained for the heats of vaporization of deuterium 
at various temperatures (1,, table 2) by drawing a 
through the single reported value, 302.3 cal 
10.70 k parallel to the corresponding eurve 
error i 


ren 


eurve 
mole at 
for hvdrogen as 
this approximation was not estimated 

The hold up nt and in the neighborhood of, the 
flood point, was exceedingly difficult to measure by 
the technique emploved here This was due to the 
inability of accurately determining the volume of 
liquid in the pot because of the violence of the boil- 
ing An attempt was made to climinate this diffi- 
culty by momentarily switching off the heater cur- 
rent and immediately the volume. How- 
ever, this method proved unfeasible, because before 
the gas bubbles ceased rising the liquid volume was 


derived from eq 2. The 


ements 


changing rapidly due to the return of reflux fro: 
the column. The difficulty in accurately 
the volume under the circumstance 
undoubtedly the major factor in accounting for th 
discrepaney in the values for the hold-up when h 
(table 1) and the deuterium 
However, in both cases, the order of magi 


reading 
described Wis 


drogen table 2) we 
used 
tude is believed to be correct 

The results of the “static” hold-up determinatio 
are shown in table 3 In these experiments hydroge: 
deuteride was used. The results were used to cor 
rect the observed values of the hold-up in both t! 
hydrogen and the deuterium experiments I} 
figures for the hold-up in tables 1 and 2 are th 
corrected values 

In computing the the stat 
hold-up, experiment 4 of table 3 was rejected 
being evidently in error. The much higher p: 
centage of hydrogen present in the charge 5% 
opposed to 1.5% in the other three 
tend to give high results, since, through fraction 
tion, this material would be localized in the colu 
The average of the first three experiments was the 
fore taken as the correction to be applied in tab 


correction due to 


cusses wo 


land 2 


4.2. Heat-Leak Measurements 
a. Calculation from Still Dimensions 


At an atmospheric pressure of 753 mm, the p 


sure inside the still was 771 mm, with no cun 
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flowing in the heater coil This pressure differential 


corresponded to a temperature differential of 0.073 


irom eq l Since the condenser wall thickness / 


was 0.075 em, and its area, «4, was 493 em the 
rate at which heat was entering and leaving the 
system to maintain a steady state was obtained 
from the equation 
dQ dT 7 
KA-| ) 
dt dL 


l'sing a value’ for A of 3 10 cal em To deg 


this rate was 0.144 cal sec, or from eq 4, equivalent 


to a boil-up rate of 70 mihi 


This value is undoubtedly too high, since the 
heures used for the effective area ol the condensing 
surface is excessive This was indicated by the 


fact that the pressure in the still was not dependent 
level of liquid hydrogen coolant in’ the 
condenser over a range of levels there 
was no heat imput. In addition, the condensing 
iris condense! wall and condenser wall coolant liquid 


upon the 


wict when 


temperature differentials were not taken into account 


These differentials may exceed that across the con 
denser wall itself, though the latter is the one em 
ploved in the d7/dZ term in the eq 5 
b. Time Method 
he “tine method deseritbed above vielded 
values of 2.5 and 2.0 ml/hr, for the boil-up rat 


the results of 


Though these values possess 


equivalent to the heat leak as two 
table 5 


consistency 


experiments 
the merit of objections can be raised to 
the exp rimental The 
mpl it im the method is that when the pressure Was 


reduced, the entire body of liquid in the still was in 


method used assumption 


thermal equilibrium This was probably not real 


ed in practice Sines cooling through evaporation 
would occur at the surface of the liquid, and the 
re mamder would have to be cooled by eonvection 


\loreover 
1 absent 


thermal equilibrium would most likely 


throughout the entire experiment, as the 


convection would be slow 
In addition material in the column 
and its thermal state known 
considerations, the validity of the time method 


or measuring the heat leak is open to question 


pressure rose, since the 


the amount of 
Because of 


was not 


thes 


c. Withdrawal Method 


In the opinion of the authors, the best method for 
measuring the heat leak is the “ withdrawal method” 


bet the k . 


at which liquid is vaporized at 
pressure as a of the heat leak 
The tuble 4 
the question of thermal equilibrium did 


in which the rate 
constant consequences 
measured results are 


In this case 


Wiis yivelh in 


not arise, stace the entire process occurred at con 
stant temperature Also, the quantity olf material 
in the still was of no conse uence since the heat 
leak was probably independent of the height of 


liquid in the still pot 


Parue 4 Il ’ } j 
\ 
I i h 
4 ‘ ‘ 
x4 
The inconsistencies found in the values in table 4 
can probably be explained by taking into account 


First, the Toepler pump 
It was operated man 


two experimental factors 
Was an inaccurate flowmeter 
ually, the mereury level being lowered at such a rate 


still this 


that the pressure remained constant was 
quite difficult) to accomplish at a uniform rate 
Second the condition of thermal equilibrium was 
undoubtedly disturbed each time the pressure was 
reduced or permitted to thecreus¢ for the reasons 
pointed out in the discussion of the time method 


However, the values were sufficiently close together 


to safely assert that the heat leak was equivalent lo 


a boil up rate of about 10 ml hn This was shown to 


be at least qualitatively true by the fact that the 
vigor of the boiling in the first hold-up measurement 
at 33.5 ml/hr (table 1) was verv much greater than 
that at zero current input At any rate, it is certain 
that whichever of the above heat leak measurements 


is nearest the truth, it is a factor of litth 
at the higher boil-up rates 


COTS apie nee 


4.3. Distillation Measurements 


Two distillations were earried out at boil up rates 
of 465 ml/h 190 mUhr. The results of these 
are listed in tables 6 and 7 


and | 
and are shown in figures 
Sand 4, respectively. Because of difficulty 
in manual control of the level in the Toepler pump 
and the continuous adjustment of the take-off regu 
lator reflux ratios were dificult to maintain 
though magnitude was the 


In figures 


constant 
the 
both CUSCS 


order of same on 


As can be seen from figures 3 and 4, the efficiency of 


453 











3 


COMPOSITION OF OISTU.LaTE sy) 
COMPOSITION OF OISTLL ATE (% My) 


$ 

















10 20 30 = ——_ 3 


PERCENTAGE OF CHARGE DISTILLED PERCENTAGE OF CHARGE DISTILLED 


separation is greater at the lower boil-up rate 
the lower boil-up rate, the shape of the curve sho 
that the efficiency of the column was above 
minimum number of plates necessary for separatio 
Since the reflux ratio was not absolutely consta 
the fact that the points fall on a smooth curve sho 
that the still efficiency is not dependent upon ref! 
ratio at this boil-up rate 

At the higher boil-up rate, the efficiency of sep 
ation is lower. In figure 4 it is seen that the pon 
do not fall on a smooth curve, probably because 
the dependence of separation upon reflux rativ 
Thus, at a value of 41.1 percent distilled the refl 
ratio rose to 870 to 1 because of a time interval spe 
in charging sample bulbs; the fractionation of t1 
sample was thus improved noticeably by the refluxin 


4.4. Estimation of Column Efficiency 


In usual distillation practice, the efficiency of a 
column is obtained by analyzing simultaneously 
withdrawn samples from the pot and the head, an 
computing by analytical or graphical methods thy 
number of theoretical plates required to effect th: 
observed difference in composition |4] This tec! 
nique was not applied here because of the exper 
mental difficulties involved in withdrawing a repr 
sentative sample from the pot Instead, tl 
distillation data were used to calculate the efficien 
of the Stedman column employed in this work. TI 
calculation presented here is that of George Webb [6 
of Hydrocarbons Research, Inc. In brief, the metho: 
involves computing the composition of the pot liqui 
at any chosen point in the distillation by making us 
of the known composition of the material already 
distilled, together with inferences concerning tl 
composition of the column hold-up. The latter 
assumed to vary logarithmically from the pot con 
position to the head composition. From the know 
edge of the analysis of the material in the pot ar 
that of the head at the chosen point in the d 
tillation, the efficiency of the still can be caleulat 
by conventional methods 

Thus, consider the case, in the second distillatio 
(table 7 and figure 4), at the point where 8,000 n 
(STP) have already been distilled. The chars 
initially was 1.77 moles (39,700 ml STP), and tl 
hold-up was 16 ml liquid, or 12,700 ml STP. The 
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, Bottoms = 39,700—8,000— 12,700= 19,000 ml STP | Since the reflux in this experiment was not total, the 
! efficiency of the column is undoubtedly higher than 
| Top composition (from an enlargement of the | that indicated in the above calculation Caleula- 
appropriate portion of fig. 4)—89.3 percent of H,—y. | tions made for amounts of distillate other than the 
| Distillate composition (average of composition | 8,000 ml employed in the above example also gave 
] from zero to 8,000 ml distilled 94.7 percent of | about lin. forthe HETP. In view of the simplifying 
| I] If x is the mole fraction of hydrogen in the pot | assumptions made in these calculations, it is felt 
] 4 liquid, that no valid conclusions can be drawn concerning 
! 0 893—, the relation between the HETP and the othe 
19.0002 2 700 0.947)(8, 000 variables of the experiment 
0.893 
| In 
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O.S9S3 0.958 
log | 0.107 )( 0.042 ) 
log 1.46 
13.8 > o 
ol 
} 11.8 for the packing used 


in this work 
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The temperature of a flame is something extremely 
cliffie ult to define This should hee ome ey ident “ hen 
it is pointed out that liberated in the 
combustion process appears in many forms and that 
at a given instant a state of equilibrium can hardly 
be expected be Thus 
formed as products of combustion will have varving 
amounts of translational 
each form of energy leading to a correspond- 
temperature \ detailed of flame 
temperatures and of several methods of measuring 
them is to be found in a book by Gaydon [1] The 
apparent temperature of a flame will necessarily de- 
pend on the method of measurement chosen, because 


the energy 


to realized the molecules 


vibrational, or rotational 
energy 


ing discussion 


in general the method will utilize a physical process 
corresponding to a characteristic temperature. Sev- 
eral methods for the measurement of flame tempera- 
tures have been devised For flame 
temperatures 1,700° C the problem is rela- 
tively simple because high melting-point temperature 
be inserted directly into the flame 
For higher temperatures some method utilizing radi- 
ant energy must be emploved. Some of the methods 
sodium vapor, 
ratios among 
The last 
called the 


or proposed 


below 


detectors can 


generally by 


and determinations based on intensity 


are line reversal using 
the components of the emission bands 
method what 
rotational temperature 


named leads to mav be 





The present study is concerned with a special class 
of flames, namely flames 
The of such a flame depends on the 
abundance of incandescent particles in the flame, 
which may be either carbon particles originating in 





opaque Or semiopaque 







appearane c 







the combustion process or solid particles of some 
other material purposely introduced. For practical 
purposes the radiant energy emitted by a cloud of 
glowing particles may be regarded as having proper- 
ties similar to those of the spectral continuum from 











the surface of an incandescent body at the same 
A partial report of this work was given at the New York of the Optical 
Society of America (M “My 1. Opt . 30, ¢ nu 
I r n brack t } erence na 


ft paper 








Use of Radiation From Incandescent Particles as an 
Indication of Flame Temperature 
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Earle K. Plyler and Curtis J. Humphreys 


This study wa indertaken to determine if the radiation from incandescent carbon 
particles could be utilized in the measurement of the temperature of flames This is a 
preliminary nvestigation of the measurement of flan temperatures by using the intensities 
in the ultraviolet continuum at two wavelengths chosen in regions comparatively free from 
molecular bands and comparing these energies with the continuous emission at the same 

. wavelengths of a calibrated tungsten ribbon lamp containing a quartz window Fuel-rich 
oxvacetviene flames were used because they contain the desired incandescent particles 
The ratios of intensities in the chosen spectral regions were determined photographically 
and photor lectricalls Interpretation of the observed ratio for a given flame, in terms of 
temperature has bee ade, and flame temperatures from 1,900° to 2,400° © have beet 
determined 
l. Introduction iemperature. The essential difference is that tl 


over-all intensity, or total radiant flux, will depen 
on the number of particles present and will alwa 
be less than that from the surface of a solid or liqu 
at the same temperature 

The spectral distribution of energy from a hot 
body emitting continuously is a function of tempera 
ture. If the emitter is a black body, this tempera 
ture dependence may be expressed by any of a number 
of well-known formulas [2], the Planck distributios 
law being the most precise. The form of this distrib 
tion curve, in which energies are plotted as ordinates 
and wavelengths as abscissas, is somewhat similar 
a probability curve. At any given temperature 
there is a maximum of energy emission corresponding 
wavelength. The maximum of 
shifts to shorter wavelengths for higher temperatures 
according to the well-known Wien displacement law 
\7T=K, where X is the wavelength maximum 
emission and 7 the absolute temperature. It ts 
obvious that if the wavelength position of the emis 
sion maximum could be located precisely, the corr 
sponding absolute temperature would be determin 
Unfortunately, 
maximum is beset with several experimental diffic 


to some emission 


the precise location of the emissior 


ties, not the least of which is that it is broad and 
not too well defined (nother method of tempera 
ture measurement is to establish the form of th: 
energy distribution curve. This may be accomp 


plished by determining the relative values of energ 
emitted at a number of wavelengths. When visibl 
light is emploved, this is the procedure used 
determining what is called the color temperatut 

p. 23 and 382] The method depends on the con 
parison of intensities of different portions of th 
continuous spectrum, as isolated by color filters. Ii 
application to flame-temperature determinations 

discussed by Gaydon [1, p. 175]. The temperatu: 
measured is that of the glowing particles in t 
flame, which may not be completely in equilibriv 

with the flame gases. Emission 
molecules excited in the combustion 
possible source of error, which can be avoided 


characteristics 
process is 




















using filters that transmit only in wavelength regions 
free of these characterist spectra 

The representation of spectral energy distribution 
from continuous emitters by a precise formula is 
ordinarily Corrections. 
depending on variations in emissivity, 
interpreting spectral distributions encountered in 
actual glowing bodies in terms of temperature, pre 
sent no great difficulty 

For this study of flame temperatures 


restricted to black bodies 


necessary to 


the ultra- 
violet region was selected because: (1) there is no 
appreciable radiation from surrounding objects that 
have temperatures below 1,600° C, (2) the energy 
emitted increases more rapidly with temperature 
than it does in the infrared, and (3) the difference 
between true and brightness temperature is less than 
n the infrared [2, p. 380] 

Although the procedures for measuring color tem- 
peratures are well known, very little work has been 
reported previously on the measurement of flame 
temperatures by establishment of the energy distri- 
bution in the ultraviolet 


2. Experimental Procedure 


Oxvacetvlene flames were selected for use in these 
experiments, a choice based on the availability of 
equipment for producing and controling such flames 
and the fact that rich mixtures of oxvgen and acety- 
lene burn with an opaque flame owing to the large 
amounts of carbon liberated, a condition favoring 
this method. The essential units of a commercial 
welding outfit, consisting of a torch with a number 4 
welding tip, regulators for the respective gases, and 
vas c\ linders, were used to produce the flames Both 
the volume of the flame and the relative amounts of 
the two gases are controllable by manipulation of the 
valves in the mixing unit of the torch and finding the 
appropriate settings of the needle valves in the regu- 
Rotameters were used in the later stages of 
the work in order to measure the rates of flow of the 
The appearance of the flame varies greatly 
with the mixture ratio When pure acetylene is 
burned it is vellow and soot, With small amounts 
of oxygen it is a brilliant white. As the oxygen 
flow ts increased the white appearance thins out, and 
the base of the flame becomes blue in color Further 
additions of oxygen are accompanied by a complete 
disappearance of the white color, the flame becoming 
blue and essentially transparent throughout 

A ribbon-filament lamp was used as the reference 
temperature standard. These lamps have a heavy- 
duty bipost base, a silica window attached to a side 
arm by a graded seal, and a maximum current rating 
of about 40 amperes A calibration of the lamp used 
was furnished by the Bureau’s pyrometry laboratory. 
This consisted of values of filament current corres- 
ponding to brightness temperatures at 100-deg 
intervals from 1,900° to 2,400° C. Following the 
usual pyrometric practice these brightness tempera- 
tures were measured at a wave length of 0.65 «4. The 
uncertainty of the calibration at 2,400° C was 
estimated at 7 deg C. The highest temperature 


lators 


USES 


included in the series of calibrated points represents 
the highest safe operating temperature of the lamp 
Oxvacetvlene flames giving 
these limits are essentially opaque, since they occur 
for mixtures relatively rich in fuel 

The practical application of the method of measur- 
ing temperatures herein described requires the 
development of a simple portable radiometric device 
The basic component around which the device was 
actually an RCA 9385 phototube since 
the temperature evaluation depends on the deter- 
mination of the form of the spectral energy curve, 
it is desirable to obtain readings at three wave- 
lengths Readings at two wavelengths separated by 
an interval of a few hundred Angstroms are, however, 
sufficient to establish the slope of the curve over the 
interval, and provide the basis for a temperature 
calibration 


temperatures between 


built was 


In order to isolate two selected wavelength regions, 
two narrow-band-pass filters were prepared for use 
with the phototube Each consisted of a combina- 
tion of an interference filter with a glass filter. The 
interference filter, which is essentially a Fabry-Perot 
interferometer with a fixed separation, has a narrow 
transmission band at a specified wavelength. There 
are, however, secondary maxima at longer and 
shorter wavelengths. These may be eliminated by 


combining the interference filters with a= glass 
filter, which isolates the region in which 
the principal maximum of the interference filter 


occurs. The following combinations were adopted 


interference filter having a transmission maximum at 


3541 A, with a Corning glass filter, code number 
9863; interference filter, transmission maximum 
$140 A, with a Schott filter designated BG28. These 
filters were selected after examination of flame 


spectrograms showed that the spectrum was sensibly 
free of characteristic molecular bands at the wave- 
lengths coincident with the transmission maxima of 
the filters The measurements were thus confined 
to the continuous emission o1 ivwinating In the clow ny 
particles in the flame 


A small housing was constructed to serve as a 
light shield and permit placing the filters against a 
window opening in front of the cell cathodes The 
remaining components of the metering device con 
sisted of a balanced d-c amplifier described in an 
earlier publication [3] and a 
resistance network to vary the sensitivity 


calvonometet and 


The performance and reproducibility of readings of 
the phototube meter were tested by a series of ex- 
posures to the ribbon-filament lamp Readings were 
taken with each filter combination in place for each 
temperature setting between 1,900° and 2,400° C 
The results of this test are displaved in figure 1 
of the curves are plots of deflections col responding to 
photocell readings for the different temperatures for 
the respective filters The third curve represents the 
ratio of deflections at each wavelength for the two 
filters. The variation of this ratio with temperature 
provides the basis for a calibration of the meter for 
temperature readings. The data 
figure | are for one set of observations rather than for 


Tw 0 


represented in 
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several lt is to bye noted that the 
extremely regula he 


nh avernuge ol 
curves are 
deflections 


representing 


rapid increase in 
with temperature, shown by the curves 
filter readings that radiant 
energy in the ultraviolet provides a sensitive scale of 
temperature If one had 
tures lower than 1,900°C 


indicates 


been measuring tempera- 


it would have been advan- 


taveous to select two wavelengths im othe visible 
ol infrared inh ordet! to have hore radiant energy 
available 

\ series of observations were made on the flame 


using the same filters as were emploved for the tung- 
The flame was changed by using a 
constant amount of C,H 
of the mixture In firure 
radiant energy emitted through the filter at 


sten mbbon lamp 
and varving the O, content 
2 are shown the ratios of the 


S541 A 


to that emitted at 4140 A for different amounts of 
oxvgen All the observed ratios fall on a smooth 
curve except for the fourth determination Ocean 


sionally the flame becomes unsteady and some varia 
tion in the ratio The conditions for the 
reproducibility of the flame will be discussed lates 
When the temperatures for the six points in figure 2 


results 


are determined from the calibration curve of figure 1 
it is found that thes at about 1.850°C for the 
first point and reach a value of 2,500°C for the sixth 
pone 


Start 


In order to have another method of measuring the 
emitted 541 and 4140 \ the 
spectrographic method was emploved At the time of 
measuring the radiant with the photocell 
exposures were made of the flame spectra by the use 
of a Hilger E2 spec trograph On evaluation of the 
plates it found that the temperature of the 
second point of figure 2 was 1,950°C, the third close 
to 2.000 . the fifth, 2.300 ~ and the 


ratios of energies at 


energy 


was 


the fourth 2? 100 


sixth somewhat above 2? 10) 


The method of evaluation of flame temperature by 
the spectrograms ts illustrated in figure 3 which shows 
a set of spectra of the flame and the filament with 
comparison spectra of the mercury are 
reading from the left, at 2537, 2804, 3132 and 
1358 A are from an Hg are and are used for the wave- 
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length determination of regions in the flame spectra 
The spectra are displaved in groups of three. The 
flame spectrum under selected operating conditions is 
at the center of each group and is flanked on either 
side by filament exposures at the same temperature 
but with different length of exposure The three 
flame spectra were produced under three different 
operating conditions, and are compared with filament 
exposures at 2,400, 2,100, and 1,900°C, respectively 
reading from top to bottom on the plate The fila 
ment exposures adjacent to each flame spectrum ar 
in the ratio of 1 to 4. The highest filament tempera 
ture is clearly higher than the highest flame tempera 
ture recorded on the plate For the lowest tempera 
ture the match is fairly The match is good 
for the intermediate temperature when the longest 
filament exposure is compared. With some exper 
ence it 1s possible to obtain much better matches 
than shown on the plate. Only the continuous por 
tions of the flame spectra are considered in maki 
these comparisons 

The mertching of was carried 
visual examination only, this being considered sat 
factory in view of the limited objectives of the inves 
tigation An alternate method that certainly woul 


close 


intensities out by 
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be more precise because it is entirely objective would 
be to determine the plate densities with a densitom- 
eter for both the flame and filament temperatures 
and considet the criterion of equal temperatures to 
be a constant ratio of the intensities between two 
wavelengths for the flame and for the lamp How- 
ever, this process of matching exposures on photo- 
graphic plates by a spectroscopic method, which is 
tedious and detailed if carried out precisely Is not 
required for temperature measurements, but was 
used as a check on the photocell method 

The radiation from the part of the flame that was 
2 in. above the tip of the torch, was focused on the 
photocell so that the results from separate trials 
could be compared Figure 4 shows a plot of the 
ratios of the energies for different fuel mixtures as a 
function of the fuel mixture. The composite curve 
is intended to be the closest possible representation 
of the observations taken on several different days 
each point im the graph corresponding to an observa- 
tion. The same symbol is used to distinguish ob- 
servations on a given day The spread of the points 
at any one fuel ratio is produced largely by imperfect 
duplication of the fuel mixture for the separate 
determinations \ small change in the fuel ratio 
would change the ratio of energies by several percent 
\ smooth curve drawn through these points would 
give an arbitrary scale of temperature. This curve 
would also make it possible to operate the flame 
under the same condition at different times 

The temperatures estimated on the basis of the 
observed ratios and the calibration curve of figure 1 
are in terms of the brightness temperatures of the 
filament of the tungsten lamp. The true tempera- 
ture of the filament as calculated from the Wien 
equation in degrees Kelvin corresponding to 2,400 
( brightness temperature is about 2,900° K. If an 
emissivity of 0.7 is assumed for the particles in the 
ultraviolet they would have a true temperature of 
2.830° K when the brightness temperature of the 
filament is matched, and with the emissivity of 0.8 
the true temperature would be 2,810° K With 
almost equal fuel ratios the OH spectrum of the oxy- 
acetylene flame has been recently measured [4] 
It was found that the OH vibrational and rotational 
temperatures corresponded to 2.600 kK with possible 
error being up to 200 These results indicate that 


the temperature determined by the solid particles 
is in satisfactory agreement with the OH temperature 
determination Also, the spectrum of C, emitted 
by an oxyacetylene flame has been observed by the 
use of a grating spectrometer in the near infrared [5] 
The temperature deduced from the rotational spectra 
was of the order of 2,700° K. This value is in good 
agreement with results obtained on the emission of 
radiant energy by solid particles as found in the 
present work 


3. Discussion of Results and Evaluation of 
Corrections 


[ p to this pomnt in the dese ription of these « \peri 
ments no discussion has been introduced regarding 
the limitations on temperature matching at differ 
ent wavelengths when nonblack bodies are being 
compared. It is apparent that a change of emissiy 
itv with wavelength will change the form of the spec 
tral enerev distribution curve and thus prevent co 
incidence of curves from different sources, exe ept 
for very limited ranges. The change of emissivity 
with wavelength usually proceeds at a very slow 
rate so that mismatch from this cause ts slight over 
the range of wavelengths used in these experiments 
It. is to be noted also that, if one is comparing diff 
erent radiators each with constant but different emus 
sivities or if one is comparing one such radiator with 
a black body where the emissivity is always unity 
it is still impossible to get an exact energy match at 
two wavelengths This is a consequence of the fol 
lowing equation based on the Wien distribution 


formula 


l | \ log « 
Yr S, C.loge 


where 7’ is the true temperature, Sy is the brightness 
temperature at wavelength \, C, the radiation cor 
stant, and ¢, the spectral emissivits \ conversion 
by means of this equation of the calibrated values 
of brightness temperature of the ribbon-filament 
lamp, at 0.6650 yu, to the corresponding values of 
0.4140 and 0.3541 uw, based on accepted vaiues of 
emissivity of tungsten has been made. The pro 
cedure is to calculate the true temperatures corre 
sponding to the calibrated values of brightness 
temperature, then taking these true temperatures as 
the initial data, substitute in the equation to obtain 
the brightness temperatures at the wavelengths of 
the respective filter transmissions. The following 
typical example is selected to illustrate the mis 
match at these wavelengths. Given a true temper 
ature of 2,713° C, the brightness temperature of the 
filament is 2,520° C at A\—0.4140 pw and 2,547° C 
for X 0.3541 uw when ¢ 0.45 
of this is that, in comparing the intensity of the lamp 
with that of a black body, if the energies correspond 
at 0.4140 4g, the energy would be greater fos the 
lamp at 0.3541 w and make it appear brighter by 
27 deg C. This difference is less than the exper 
mental error when observations are made on the 
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flarne It can be shown from eq | that the differ os the errors introduced by the observed ratio 
ence between the true and the brightness tempera intensities 
tures decreases in the ultraviolet revion a compared I Larne temperatures determined by thus metho 
with the visible spectrum. When the emissivity in are in essential agreement with temperatures 
creases, the difference between the true and the | duced from = spectroscopic observations of flan 
brightness temperature become le While the | produced under similar conditions 
emissivity of carbon particles is not well known at 
high Lemperatures W. Ek. Forsythe ) p. 380] gives 4. References 
the following velues for the emissivity of carbon 
me (A 0465 pw), at T1,200°K, is O86, and, at | [1] A. G. Gaydon, Spectroscopy and combustion theory (( 
7’ 2.200°K.. is 0.79 man & Hall, London, 1948 
, 2) W. bk. For the, Measurement of radiant energy, chapts 
Calculations have been made on the temperature MeGraw-fHill Book Co.. Inc.. New York, N. ¥.. 1027 
of the purty les on the basis of « bemyg equal to OF 1 W. W. Coblentz and R. Stair \ portable ultraviolet 
and OS When the brightne temperature ts 2,547 tensity meter, consisting of a balanced amplifier, pl 


electric cell, and microammeter, BS J. Research 12, 2 
. 1 9 6 ‘ . 1034) RPOA7 
are 2,618° and 2,591° © for emissivities of 0.7 and | ty) wos Renediet. E. K. Plyler. and C. J Humphreys, \ 
O=.S, respective ly If the emissivity of the parti les bration-rotation lines of OH in flame spectra, Bul. A 
were 0.45 mstead of O77 at O354y the true temper Phy Soc, 26, 50 (1051 

f 1D tenedict anne art "Ivler Sear of ¢ and ¢ 
ntiuyre would bye O5 deg hivher lor thre nine brohtness W “ [a 1 barl hk. Plyl Band ! it : 


C' the true temperatures ata wavelength of 0.354 yu 


in infrared flan pectra, bul \r Phi Soc, 26, 
temperature Until more accurate values of emis 1O51 
ivity are known there will probably be an uncer 
tninty in the estimates of temperature of about 
100 deg © This is of the same order of magnitude Wasuinearon, August 29, 1950 
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Electron-Optical Bench 


L. Marton, M. M. Morgan, D. C. Schubert, j. R. Shah, and J. A. Simpson 


\ rsatile electron-optical bench ha 

| tl ¢ ent Phree irriag 
ine byective smd meshe ire arranged 

t vie for external positioning W 

l¢ ‘ re perimental fechnique are ali ( 
foca er 1a magnet el i presenter 

beneh and t licate the limits of accuracy 

l. Introduction 

For the extensive study of the properties of 


an electron opti al bench 
his bench provides for the 
of such elements from outside 


ele etron Opty al ele ments 
ed 


positions 


has by en construct 


control of the 


the vacuum system and quantitative measurement of 


the image formed on a fluorescent sereen 


Verv few uttempts have 
Among the very few electron 


been made to design an 


clectron-optical bench 


optical benches, which have been built in the past, is 
the instrument designed by Reisner and Pieard \ 
minor feature of that bench was the provision ior 
visual observation of the parts mside the vacuum 


‘ haart I 


Inasmuch as all measurements required during the 


operation Ob an electron opti al bench can be made 
the ability to see the optical system need 
For simplicity 


X-ravs the 


was made with an all 


externally 
not be considered particularly useful 
of construct nna against 
NBS electron opty al 
metal evlindrical vacuum chamber 
the 


ion protection 


bye 1h hy 
lt was designed 
The bene hy 


to satisty following requirements l 


should accommodate three carriages for lenses with 
maximum diameter of 7 m. and four holders for 
apertures test objects anal meshes used for the 


electron-optical shadow method 

>» A movement of 10 in. along the axis of the bench 
in racially by 
the svstem without breaking 


moa should possible for each eom 


ponent of the vacuum 


2. Description of Apparatus 


2.1. Vacuum Chamber 
\ rolled-steel eylinder of 15-in. inner diameter, 
in. wall and 39-in, length is part of the vacuum 


chamber used to house the electron opti al system 
The end plat 
fixed to the evlinder with bolts, the vacuum seal being 


supporting the gun is permanently 


provided by a rubber gasket. A dural plate, A 
figure 1, on wheels serves as a bed on which the 
clectron-optical system is mounted. The specimen 


B 


along the bed 


and lens carriages, C, slide 
The bed Thay be rolled 
One end of the bed is 


and aperture mounts 
ha groove 


into or out of the chamber 


H. Reis: Rr | | I 1”, ‘ MA 
L.. Mart : i Applied I 20, me) 34 
NBS 43, 400 (1040 


“4 a i al feor hig ute V« i of 

or magnetic lenses and four holders for aper 

ippropriate 1 othe rcuu chamber, while 
$} degrees of freedom for each element rt 


ed in detai 


\ set of mea 


istrate the 


i in be atta dd 
attached to the fac plate by means of a flexible 
binngre which does not encounter the stram that 
might result if a rigid jomt were used 

The face plat seals against the end of the tube by 
® heoprene rasket compressed by atmosphert pres 


sure The fluorescent screen ois coated on oa lass 
plate (54, in. in diameter and ‘-in. thick) in the face 
plate Motion of the elements is controlled by means 


of connecting shafts pPaussiny through “Wilson seals 
IK, in the face plate 


Lens Carriages 


The 
properly combining controlled vertical and horizontal 
the kor 
ment a doubl \ 


echnyaves a ve 


radial adjustment of each lens is obtained by 


motions of lens carriage vertical adjust 


2, attached 
the shaft 


figure 
ar, © 


threaded worm 


B 


toa connecting shaft 


of which rotates in the bottom dural plate 1) \ 
screw pinned to the worm gear turns mm thy top plate 
I, thus moving it vertically when the worm is 


rotated The top plate is constramed in its motion 


by sliding on four pins, F, attached to the lower 
plate. Horizontal adjustment is similarly achieved 
by a worm and gear arrangement sliding the plate, G 
ove! plat KE, the motion being constrained by a 
tongue in groove Since the vertical position of the 
worm, H, changes, two universal joints are inserted 
In its connecting shaft The axial position of the 
carriage os controlled by pushing ol pulling the 


the Wilson seals 


connecting shafts through 





chambe ‘ elect 


hanging element 
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2.3. Specimen, Mesh, and Aperture Mounts 


If one is to achieve a smooth motion of specimens 
or test objects with a minimum of jerkiness, the 
drive must be made as nearly frictionless as possible 
Pursuant to this requirement, the radial motion in 
the present arrangement is constrained by elastic 
deformation of phosphor bronze strips instead of 
sliding surfaces. A similar arrangement for posi 
tioning an electron microscope stage, which reduced 
friction to an even greater degree, was described in 
an earlier papel 


V ,, A ed Phys, 16, 151, 104 
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Radial adjustment of the specimen, mesh, o 


aperture hola rs is provided by nm put of worm and 
year mechanisms at right angles to each other \ 
worm, A, figure 3, attached at the end of a connect 
ing shaft engages a gear, B, rotating in a collar, ¢ 
fixed to the mount \ screw 
the worm gear shaft, turns in a sleeve, D, which 
attached by means of a flexible phosphor-bron 
strip, EK, to the ring, F, into which an appropriat 
mesh, or aperture holder is screwed lt 


pinned at the end of 


specimen 
this way the ring is moved along the axis of th 
worm wear B The elastic deformation of th 
phosphor bronze connec tion, because of its onmenta 
tion, provides a simple solution of the mechanica 
without th 


lokage problem of movitiy the rig 


use of friction devices, in a plane perpendicular to 
the axis of the bench while maintaiming rigidity u 
the axial direetion Because of the arrangement 
of the Wilson seals, the connecting shafts had to b 
attached to the circular plate forming the body o! 
the mount near its top, thus leading to somewhat 
jerky axial motion. This jerkiness was reduced by 
using a long, heavily weighted wheel base for th 


2.4. Vacuum System 


diagram of the vacuum system 
The lens currents are led through 


\ schemati 
shown in figure 4 
the vacuum wall by metal-glass seals in the face plat 
and the conductors are tnsulated by ceramic beads 
which produce less outgassing difficulties than would 
result if rubber or plastic insulation were use 
This is realized by means of a large fore-vacuun 
pump and a 4-tn. oil diffusion pump. This larg: 
pump brings the pressure down from the atmospheri 
to about LO microns of mercury in 5 minutes. How 
ever, as this pump causes considerable vibration 
it is turned off after mitial pumping and the quiets 
small fore-vacuum pump is used during final a 
The normal opera 


justments and measurements 
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ing pressure | about 310 microuw of mereury 


measured by an ton gauge near the gun 
The special three-way valve near the small capac 
ity pump ts cotstructe dl lo bole ‘ al nir mito tha Vacuum 
well as to tosert 


system at oa controlled rat us 


the small 


capacity pump iato or isolate it from: the 


system 


2.5. Electrical System 


{ conventional L00-kilovolt’ X-ray power supply 
unit. consisting of a voltage doubler cireuit with 
two 50-kilovolt reetifier tubs is used as the high 
voltage upply Stabilization is provided by two 
O.l-miecro-farad condensers in parallel with the eu 
euit The voltage 1 measured by the current 
passing through a |.000-megohm bleeder resistances 
connected  aeros thre high-voltage supply I ha 


normal range of operation i from 20 to 70 kilovolts 
The filament is heated by a high-frequency powe: 
supply Ihe lens current powell supply sare a modi 
fication of those deseribed in an earher paper (see 
footnote 3) and give up to 500 milliamperes with a 
stability of better than 1 part im LO,Q00 

The lenses available for use in the beneh were 


oriinally designed for use with the coils out ta th 


atmosphere rather than ina vacuum chamber; con 
wquently they were poorly suited to operation 
within the beneh, and the coils released large quan 
tities of gas to the vacuum system An extremely 
rapid prurenapretiy ystem wa needed fo tmaitamn oa 


working vacuum under such unfavorale conditions 
2.6. Optical System 

\ reat portion of the ise oof an electron optical 
bench is in the empirical approach to the develop 
ment of electron-optical leus system In such 
work the result of one measurement often determines 
what conditions should be used for the next measure 
ment Hlenee the direet method of measuring the 
Hage ts preftert d as quicker imples and cheaper 
than the photographic method of recording results 
\ low-power traveling telescope reading 0.001-in 
with a total travel of 2%j-in. was specially built for 
measurement of the image on the fluorescent screen 
Since the bench was designed for use with magnetic 


the image 
mounted on a turntable, 


with their attendant rotation of 


the measuring telescope 1 


ke nses 


which can be moved either horizontally or vertically a 


distance of 3 in. in such a wav that its center of rota 


tion may be concentme with that of the image In 
this wav all measurements are radial Variable in 
tensity side illumination of the reticule reduces the 
difficulty im seeing the cross hair when the intensits 
of the pattern on the sereen is low 


the fluorescent 


Zire 


screen 


Wis used for 
conducting it 
suggested by 


cadmium silicats 
To make the 
aluminized, using — the 
Bachman 
\ 100-ky 


scribed in an earlier publication (see 


screen Wiis 


technique 
identical to the one de 


electron gun 


footnote 3). 1s 





by hy obtamed 
three b5-volt 


ir] Th bin 


ofl electrons the 


connected acros 


thre soOuUres 
from 
batteries im 


nu pote hntiometet 


2.7. Measurements 


tuken oon tha cleetron 
graphically ») for 
The  Thewasiure 


Sample Inensturement 


optical bench are presented 


the purpose of illustrating its use 


ments of foeal length were taken with the condense 
lens, a mesh object, and the objective lens in the 
svstem The objective lens ha 10.0000 turn \ 
symmetrical cold-rolled steel pole prece of O.OS1-in 
bore diameter and O.060-in. gap width fits mito the 


ol thy len 
obtamed by 


in. immer diameter of the tron enclosure 
D>itTuse 


applying a 


Hlhimination of the object 


than film anda thin condueting 


Korn Vial 


alumimum coating on the source side of the mesh 


was used to make focusing eritical The tage 
distance was held constant throughout the measure 
ments; for a particular accelerating voltage — the 
object position was adjusted to give approximate 
focus for different objective current the final foeu 
beimy maude by a riall ndyu tment of the current 


The focal length was computed by the equation 
j= ¢ \/ | where g ts the image distance and V7 1 
the maynification 

With the assumption that the curves drawn are 


approximately as good fits as would be obtamed by 


45 
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applying the least-squares method to a “known” 
function, the rms error in the focal length is found to 
0.001 em. From the graph the maxi 


seems about O.0l em The Ta yOr 


be less than 
mum error 
limitation on accuracy seems to be the difficulty in 
setting the cross hair on a port of the mesh Hage 
Improvement of the cross-hair design and develop 
ment of greater skill in setting may reduce the erro 
from this source to less than | percent from a present 
value as high as 2.5 percent Another limitation ts 
in the precision with which the image can be focused 
With the objective lens and pole piece previously 
deseribed, a change ol lens current of the order of 

} ma, over the range of currents used, is required 
to produce a detectable change in definition from 
best focus, that is, the image appears to be in perfect 
focus over a range of about 6 ma The center of 
this range can be located within about 2 ma lt 
much improved 
equipment 


is believed that this cannot be 
without elaborate additions to the 
Errors due to small inaccuracies in alinement and to 


hysteresis ure believed to by negligible compared to 


these two limiting factors; po particular care 
taken to eliminate hysteresis errors in these p 
liminary measurements 

The electron opty al bench is quite versa til 
cause of the ease with which elements may be ad 
and removed, and because three-dimensional mo 
ment of all elements is possible from outsicdk 
vacuum wall It can be used for almost any eleetr 
optical measurement and can be adapted for uss 
electron-optical field mapping and as an elect) 
microscope 

Practically all the construction work was don 
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Random Determinants 
R. Fortet 
Determinants whose elements involve random variables are discussed and expressio 
derived for the first and second moments (Applications are made to n-dimensional geometry 
especially, to finding limiting probabilities for the event a given point lies above a randon 
hyperplane inder fairly general hypotheses. The random variable A,h— F is considered 
where A, and # are certain minors of the determinantal equation of the random hyperplane 
and A is a coordinate of the given point 


An asymptotic expression for EB} A is obtained, and it is shown that EF! FP s of tl 
order of (lin) F 1 


| KA 0) 
Also 

In this section real random variables .\,, depend- A >. 42 ¥ Te 1 oS 
ing on two indices i and 7 and having the following 5 Fk 
properties are considered Now 

lat A,, A A, be k random variables chosen 

from the \ s Then We make the assumptions i > Dy, | A \ \ 
that yy - 

a) for any 7 and ), for anv &, and for any choice 24! bis Are Ni (A hija, \ 


of the A,’s, the conditional mathematical expecta 


tion From (a) /CN\,, A,,.A \ 0: hence 
F(X 44|A A 0 


if .\ is not one of the A,’s 
b) For any i, 7, for any &, and for any choice of 


the A,’s. the conditional second order moment and from (b) this reduces to 
KUN | A J ki A aH QN(A 1) 
where o, is a constant that depends only on +) pro 
vided .\,, is not one of the A,’s | Now A; is a A which depends only on those 
Consider the following random determinant with «#1 If we assume that the formula 
\ \ \ . 
EiA n!> Jo 2) 
\ \ \ 
is true for n=—1, 2 r, (1) shows that (2) is also 
A valid for n=r+1; on the other hand (2) obviously 


holds for n l Thus ») is established by induc 


tion, and the proof of the lemma is complet 


Wi have the 
Lemma: E(A Oand k(A n! I] o d 
Proof Let A,, be the cofactor of .\ m A: we 
can write 


Nonrandom real numbers s and the determinant 
are now considered 


suppressing the & rows 1 ’ / in d,, a 
rectangular matrix is obtained with (n—s) rows; 
From (a). ECN,,/A 0: hence Tlin.te, . ted.) is called the sum of the squares 
a : ; of all determinants of (n h th order deduced from 
——=—_ ~~ —— the Ae | this matrix 
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Let J), be the following random determinant 


/ Au a \ J \ 
! \ 

iz) 
F Kes 2 \ J \ 


DD, can be written as a sum of elementary determi- 
nants, the elements of which are either z,, or .Y,, but 
never a X One of these elementary determ- 
nants is d,: all the other elementary determinants 
are random determinants Let S be one of these 
random determinants: S consists of k rows of 


a 's h >1) and n k- rows of z,,’s For instance, 


suppose that the first row of S consists of the \,,; S 


is a linear 


homogeneous form in the X,,’s, the 
coefficients of which are independent of the .X,,’s 


Thus according to nu khiS U: consequently 
Kk(D d, 
Let i,, 2 i. be the & rows of S consisting 


of some .Y,,’s, that is, tow consists of the 

row i, of the Y,,, et 

development of S in terms of these & rows 2 lo, 
», we can write Sin the form 


Employ ng the Laplace 


S=F°A.(i3 B iF 
where A,(),.) . Je) 8 a determinant of the 
preceding A type, the elements of which are those 
Aw with 1 ¢€ %, 2 = ; 6 fy » Gat 
BO je. . Je) is the algebraic complement of 
An ()iJo, de): and the summation =* ts extended 
over all combinations (),,) ») of order fk 
taken from the n integers 1, 2 — For two 


different combinations we have 


he A,() J . Ded A, . } () 3) 
since there is at least one le. Jy SAV which is not a 
jy, 80 that the product A,( Ar (i, 
jy.) 18 a homogeneous linear form in X,,, X 

X and (3) follows from (a Consequently 


ke Ss =*B J } I: A J 


and by the lemma, 


BUS S*( BY AV aye 
i Il @ a ) 
But obviously, ‘ 
=°B > Jd=T (ii, ts, , tad, 


Hence we hav c 


E(S?)=T(i,,i; i 


~ 
~ 


Consider the product of two different elementa 
determinants S,; and S There is at least one roy 
consisting of some \,, that appear in S; and not 
S,, or in S, and not in S,;. For example, suppose th 
first row m S,, consisting of the .\,,'s, does not appea 
in S,; then the product S,S, is a linear homogeneo) 
form inh the \ 's the coeflicients of which i! 
independent of the \,,’s. Consequently by (a 
E(S,)S,)=0 and E(22?) reduces to the sum of tl 
squares of all the elementary determinants. Th 
gives us the following 


Theorem 1: 


Under assumptions (a) and (b 


vied.) IL 
al 


In (5) the summation S* is extended over all comb 
nations ; ) of order & of the integers | 
2, n; for convenience we put 


j } i 


T(1,2, nid | 


An interesting feature of formulas (4) and (5) is 
that these formulas do not depend on the probability 
laws of the V,,’s. This fact remains valid even if the 
o,'s depend on }, but (5) becomes more complicated 
on the other hand it is possible to compute the highe: 
moments of 1), by using similar reasoning and suit 
able assumptions. However, the formulas seem quits 
complicated. It also appears difficult to obtain th 
probability law of oe even under such hypotheses a 
that the .\V,, are normally distributed 


Of course if o 0» Co, Co. where ossa 
constant independent of 7, (5) reduces to 
kei IP ) y Td ke - ja 
where 
1,(d.)=2°T (2, id 
From (5) certain interesting results may be deduced 


which are perhaps already known. These results ar 
connected with the geometrical interpretation of the co 
efficients of (5 Let £,, bean n-dimensional Euclidean 
space, with orthogonal coordinates. The n numbers 
rian, La, » Tm May be regarded as the coordinates 
of a point /?, in &,, and rqg+NXa, re 2, 

F in A can be considered as the coordinates of a 
random point M, in £,. That d, is invariant unde! 
orthogonal transformations is well known If 
(OP Ps P,) is the volume in F, of the paral 


lelepiped formed by the vectors OP, OP», 
OP, we have 

d, PT firs... «Sal 
interpretation holds for D, and 
(0, MM, M,): hence E(D?) is invariant unde 


an orthogonal transformation. If the Y,,’s are nor- 
mally distributed, this property also persists. The 


An analogous 
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; remain unaltered also We conclude that the 


coefficients 7%, . %&: d,) are invariant 
T(i , 3 d,) depends only on P, , 7 
consider the set of points P Fin r defined 


as follows 
oP oP 


oP OP d oP 


oP. OP, +rn0P hee 0P 


. 


where the A are such that the vectors 0P" , 0P 


0?” are mutually orthogonal. From a classical! 
property of determinants we can replace 
P, by £... ., P” without altering 7(),, te: da) 


The invariance of 7 - d,) under orthogonal! 


transformations implies that the vectors 0P; (a=1, 


_k) ean be taken as coordinate axes Then it 


becomes obvious that 
l(i,, ; d ot r, 


the volume (0, P Fr 


volume in 2 k-dimensional subspace of F,; but it, is 


being considered as a 


also obv 1OuUS that 
0, P. P OF, F r 


and we have established 


The orem 


and ¢ msequently, 
Td > 0.P P 7) 


On the other hand, considering ?;, , P,, as fixed 
points in £,, and 2 as a moving point, (Q, P 
1 
P and >* 
7) under a simple change of the origin of coordinates 
If one puts 


(Q, Py P,) are given by (6) or 


and 
=*(Q, P repye(™. 9 


where C and C* are any positive constants, then eq 
8) and (9) define quadrics Q(7., 4; C) and Q@(C* 
respectively, of ellipsoidal type. Thus it follows 
from (5) that the relation 


E{ (Q, M M, é (10) 


where ¢€ is some positive constant also defines a 


quadric ¢ which is in general an ellipsoid) 


In Buc lidean n-space with the coordinates defined 


as above consider n fixed points ?;, 2? , P,, and 
a fixed point S on the axis Or, with coordinates 
0, 0, 0, A Let M/,, AL. , 7, be n random 


points, and consider the random plane TI, (1. e., an 
n-1)-dimensional linear variety) defined by (M,, 
V/, , M The equation of II has the following 
determinantal form 


J X J Mas teen ' 


Developing this determinant in terms of the elements 
of the first row and their cofactors, one obtains 


{ y Loz B 0 11) 
For convenience we set 
) ! \ a ? i. ae n | 
Let MW be the point () ) 0) and V the 


linear variety of 7 2 dimensions defined by 


i () / 0 (12) 


and let 7/7, be the orthogonal projection of M,; or 
WV on | An elementary computation shows that, 
in absolute value, 


a > > 
MIM, |X Vi" - 
a _———— { » 
H1i,M He,M A, 

We now suppose that 

A) the distance between P, and P, is equal to 1, 


In this case the n points P, form a 


for any 7 and j 
this poly - 


regular polyhedron in the variety s,—0; 
yi - 
hedron admits (7. ) linear varieties VV, as varieties of 


symmetry, V,, being defined as that (n—2)-dimen- 
sional variety formed by the points in s,—0 which 
are equidistant from P, and P,; (#7). Now suppose 
that V is a moving (n—2)-dimensional variety, and 
consider the quantity 


max (distance from V to P 
' 
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s attains its minimum value when V isa V,,, in which 


case s=5- Consequently in all cases there is at least 


one point P, such that the distance from V to P, is at 
least a 

Now suppose that n tends to + © with condition 
(A) always holding; we make the following assump- 


tions 
B lim Pr) Max A eal < l for any «>0; 


, l 
() there exists a positive number a and inde- 


9 
pendent of n such that 
>» 5 an i 
lim Pr? Max A a l 
/ ay] \ 
If for instance the distance from P, to V is a maxi 
mum when j=1, then except in cases of small 
probability 
a we 
H,M,>.—a 
and 
Xiel<e 
Thus we have 
Theorem 3 / nder hypothe Nes A (B and i 


=, 
lim Pr \ 0 0. 


Remark: Assumptions (B) and (C) are not particu- 
larly restrictive. They are satisfied, for instance if 


the .Y,,’s are independent and normally distributed 
with 


K(X os (14 


p being< = and independent of 7, j, n; also when the 


Vv, are independently distributed with uniform 
probability density over the interior of the sphere 


of center P?, and radius p 


» 


Denote the center of yravity of P). ' P.. by G, 


with coordinates (&, f,) given by 
e. 
. = eat 7=1,2 n) (15) 
biel 
—,=0. Let A be the center of gravity of A/, 
M, with coordinates &,+ Z, (j= 1,2, 7 Put 
l wy . 
Z A (16) 
fond 
We make the following assumptions 
(D) the distance OG remains bounded as n—-+- « 


KE) the A and b) o 


satisfy assumptions (a 


section 1, with 
p 
0 
yi 
where p- isa constant independent of 7 and n The 
algebraic distance L from G to II is equal to 
} A é B 
L 17 


K belongs to II, so that 
L Z Gh 


the vector GK are the Z,’s 
it then follows that 


But the coordinates of 
defined by (16). From (FE 


Hence Gh and a fortiori L tend toward zero in prob 
Write (17) as 


ability 


fay IS B \ 
A =9 
\ AV 1+e4 
SA? 
S° A 


We have by 


the Schwartz inequality 


S ae) <q/S ay Se 
so that 
SA 
n—1 V\ > 
ay “Ajé ‘ 


4 
"4 


> & A?) A, tends to 


Now \ S> # is bounded and \ 


zero in probability. Therefore under the assump 


tions (A)—(E) we have: 
Lemma As no, |B/A,)| tends to zero in prob 
a b ity. 


The algebraic distance from S toIlis 
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hA B hA B 
A 
‘A pe. 
Vea A \ we 9 
| 
: B l 
(Ana i) S 
SA 
- 
\ 
Suppose that A 0: then S is above Il if AA O.or 
what is the same thing, (A. BA 0. We may 


write this last expression as 


above lemma (and under the 
same hypothe SCS, VIZ assumption A)- K 

Theorem 4. If h>0, the probability that S is above il 
tends to | when? fh 0, this probabil ly te nds 


lo 


Thus accordimyg to the 


ero 
The case A=O is not clear-cut But if, for in 
stance, the \7,’s are independent with the same con- 
tinuous probability law which is svmmetric about 


the plane ] 0, then clearly Pri S v l 


Ke eping the notations of section 2, we assume that 
As | for any that the .Y 
ditions (a) and (b 1 with o 0 
p/yn where p is a constant independent of 7 and 7 
The other assumptions, (B), (C), and (D) of section 
ae will not be required We wish to obtaim some in- 
random variable 





ana s satisfv con- 


ot section 


formation about the 
I] 1,4—B. 19 


Theorem 1 implies immediately that 


Ae VA, x1) 

where V,, is the determimant 

r } } l 

I y ] 

| 

J i ] 

V’,., as we have seen, is given by 
V et” ee 


formula (33) of the Appendix 


That is to sav (cf 


if 


V \ 50-1 


Let uw be 


the determinant 


} I J 0 
J J 0 
J I 0) 
For (71°) we have at once 
KH h? RCA KiB (22 
Since B is a determinant of the form J), (section 1 


formula (5a) gives the result that 


i (B?)= >> T,(w,) k= 2: 
t=1 n* 

A, is not exactly of tvpe )),, but in entirely similar 

fashion we find 


Srv, 


er. 24 
al 2 
TN s the sum of the minors of order (n—k) of 
| obtamed without s ippressing the last column of 


V Fo 


instance 


] | (P,.Fa Ff P 
where } run through all the numbers 
.. except | , } Hence 
/ V.)=V as 25 
and 
rey () ) n—k)Q-*te+t, 26 


Substituting this expression into (24) one finds 


yy 4 ) ‘ +? ; 
ECA i ) n—k)2 i = 
; . ie l / ke) ~~ 
=e 2p 27 
- fi 
If we suppose that 2p 1, (27) vields immediately 
l 7 
I A ™ - 27a 
| Ze” 3" 
as 7 
B is of the form D,; we know that @,(C*) is a 


quadric, and it must have as planes of symmetry the 


to) planes of symmetry V,, of the polyhedron 
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ire & rs Consequently the center of Q(C* 
’ > / 
is G, and 0, (¢ s a sphere sine ( oa n, and 
we can writk 
Tu h(n)OG? + Tu 28 


where 73 (u is the value of 7,(u when OG—0 


On the other hand 
Cr,P P 


where ),, ) },-» are all the numbers 1, 2 


nexcept? lf O—G., we find (ef appendix, 


formula (36 


and 


bY py w ere , > ea I a) 
—_  * a” ae Samed -} 
I ' } Na rt / 
”) 
If 2p l 29) becomes asymptotically 
: v lp 

3 T3(w,)k! 2 > 200 

> i n2 | 2p 
We can now compute the value 7%*(w,) of Tye 


for P, =O (cf formula (38 Writing 


appendix 





T*(w,)—A(n)GP Tsu 

we have an equation which determines A, (n); we find 

. nv 2) } } h 

j yi 

k—1)!2 
and 
2 p | 2 7 } | ) ? y i 
SAMA! —- > : k(2p*) 
kel U - } 
(30) 


Asymptotically this expression becomes, for 2p°< 1, 


, 9 
‘ p l -~—p . 
S*d\.(n)k! ~ 30a) 
ya i » ) 
kel i - | ~p 


We verify that 4(B 1/n) E(A2), 
which agrees with the lemma of section 2. In the 
particular case when G=O, E(B") is of the order of 
(1/n*) EA?) 


is of the ordet of 


Appendix 
points P,, P, in ar ’ 1)-dimensional 
space E, such that ?,P,<—1 for any i and j 
(that is, these points are the vertices of a regular poly- 
wish to find ar the volume 


Given n 
Euclidean 
(iF) 


hedron P.) we expression for 


) 


| 
| 


I/n') V,.) of P, in EB, Let G,; be the center of gravity 
the points P,, P P.- h, the distance between G ‘ 
P,; and d, the distance between G, and any P We hav 
| h,,' 4 
nal) G.G G.P.=0 ; 


G G,, P, are on the same “straight line’ 


‘and PG, ,P 


near variety 





i= ?P,P.= PG Gf 
It follows that 
n l 
d / 
hi} = 
fame 


Clearly hy is known: the ree 


solved for A, 


irrence 


} l 
v2, afittie 
It follows by recurrence from (31 } 
that 
if ? 
On the other hand 
; l ’ l 
\2\ 
ane 
— l ] n 
G.G Slane a\; , 
We now compute the volume G P,P 





is the orthogonal projection of G, on the 
PP P..), we have 
G., Fil’ P's Gal | 
But / ( o that 
G1 G.G G G G 
As in > 
— l l l , 
GL, 2 ( n—k ’ 2nin—hsh 
and so 
k 
G,, P P 
We can also compute (P,, P,P P's 
P PP P P.L: V . 
PL. P,G G,:G G 
l l 
9 ( n k ) 
Pherefore 
i—k 
P,, P P, —— 


From 37) it follows that 
T#*(y 
‘ Un 
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(one-dimensio 


isa right triangle 


n—-1l\n—k l 
k i) 2-2 


henes 


being equal to 3/2 


Pp If] 
linear vari 
ty 
7 
asad 
7 

















Solutions of Ax 


\ Bx’ 


M. R. Hestenes* and W. Karush 


ait pro em . to determine charact 
1 \/ vhere 1, ? are + n Hermitian 
bre i first approximatior a second appre 
approxima th appropriate alphas, c« 
1. Introduction 
Let A, B. be Hermitian matrices of order n with B 


definits The n the characteristic vectors of 


the equation 


positive 
A AM | 
are the eritieal pomts of the “Ravleigh quotient 


i Aa 
rae ’ rx, 


Ba 


corresponding values of the quotient are the 
characteristic values \. In particular the minimum 
is the least 


Our purpose is to discuss a method of 


and th 


maximum) of greatest) characteristic 
value of (1 
finding the that is based upon this 


observation and that avoids a transformation of the 


solutions ol | 


proble aa 
The 


scribed briefly as follows 


iterative one that mav be de- 
With each non-null vector 
associate a vector 1 that 
We then pass from one approx 


next ] by 


method is an 
rewe sm ina certam sense 
the gradient of wat 
the 


imation s to means of the formula 


r’ j an a 0, 
where the scalar a may depend upon zs. The gradient 
used here is determined by the equation 


tin Ag ull Br. 3 


where G is an arbitrary positive definite Hermitian 
matrix. In computational practice G would be se- 
lected so thatits inverse G gv. (dT, the 
identity matrix In section 4 it will be shown that this 
method is convergent if the scalars a(r) 
priately chosen, and in section 6 two feasible schemes 
for this choice will be described. In general, 
vergence is established only to some, possibly inter- 
mediate, characteristic value (and Under 
special hypotheses this will be the least characteristic 
value (see 


Is known e 
ure appro- 
cCon- 
vector 


section 5 


The method has several computational advantages 
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tie rin md { the prot 

rt rice \ eralized gradic ” lefined 
Ximation an detern ed SUCCESSIVE 
nverge to a solution 


problem (1 It 
minimizes round-off errors by beginning each step 
with a new initial vector The calculations at each 
simple and identical in 
The method 


It avoids a transformation of 


stage of the iteration are 
form with those of the preceding one 
is thereby particularly suited to high-speed auto- 
matic computing machines However, it appears to 
converge too slowly to be of use for hand calculation 

When one o 
the method may be modified so as to vield a new 
This is achieved 
} for the gradient 


more characteristic vectors are known 


characteristic vector (see section 7 


by appropriately altering eq 


For arbitrary complex matrices A, B it is of in 
terest to know when the problem 

Cu MDa } 

may be transformed to one of type (1 Several 


characterizations ® are section 8S, some of 


which are of computational value 


given im 


2. Preliminary Results 


In this section we shall state some definitions and 


assemble some well-known facts on matrices No 
proots will be riven 


By a 


a , & 


vector we understand an n-tuple a= (a, 
, of complex numbers. We deal with 
the space of such vectors over the sealar field of 
complex numbers. We let 

ru a h a b a h 7] 4 4 . h, 
where € denotes the complex conjugate of the scalar 
ce. Thus (2,y Yd The length of x ts [r]= (2,7)! 
If Cis an arbitrary matrix then 


rcy C2." 


meaning, and C™* is the 
A matrix // is Hermitian 


the usual 
conjugate transpose of ( 
if and only if 


where Cy has 


f° =H 


We shall sav 


In this case (27,/ir) is a real number 


Suc " ma xamy tir 1 or 
/ 4 asil che I 1 write / LL*w6 / ilar 
L* njuga Ans] {l wud cu im 
Fox, Huskey 1 Wilkir ! hut al : i 
: is equatior Quart. J. M and A i Ma ! ‘ vis 
r Its are ¢ ely ated Hw.W mcdt, Zur Al ing der 
ad te t nwerta aS endlicher I) Matt 
Nac _ t & (1040 





that two vectors v and y are H-orthogonal in case 3. The Gradient 


ry y tla 0 The direction for which the directional derivati 
of the function, « given by (2), is a maximum w 
Two sets of vectors are H-orthogonal in case each | now be calculated. This optimal direction will | 


vector of one set is H-orthogonal to each vector of the | determined relative to the inner product (2,4 


other. By orthogonality is meant /-orthoganality corresponding to an arbitrary, fixed positive defini 
with J the identity matrix , matrix G. The generality of an arbitrary inner pro 
A matrix G is positive definite in case it is Her- | uet) has computational significance as well 

mitian and theoretical interest; im practice it is limited 

matrices G whose inverses are known. The iterati 
r.Gu 0 whenever 10 method and convergence theorems that are to follo 

later depend only upon the final formula that will | 

Let G be positive definite There exist positive obtained for the maximizing direction, not upon tl 

numbers m(G) and A/(G) such that derivation of the formula; the derivation is intend: 
to suggest the motivation for the method 

m(G) r*<(r.Gr)< MIG) 2 5) For fixed vectors ++#0 and ér+0. consider tl 

function u(2+-edr) for reale. By a simple calculatio 

Also. we have the inequality we find that at e=0, 

rary <(r Grr) (yy 6 dy 2R i br,t 
de r, Ba 

the ‘ quality holding if and only ifs and yuare linearly where 

dependent Further, the matrix G is positive t= 8(2 Az—pis) Be. _” - 

definite, and there exists a positive definite matrix 

, such that G=@ and Ric} denotes the real part of ¢ We therefor 

We turn now to problem 4), where ( and D are | cook that vector ar for which 
arbitrary matrices The number \’ is a character 
istic number (root, value) of (5) in case there is a Ri (ax.¢ max br. Gar) —1 g 


non-null vector y’ such that 


n is defined by the equation 


( y" MN Dy’ 

: Gn=t=Ar—yp Ba 

We allow the characteristic value \ ; im this | Then. using (6 

case Dy’ =0. Wesavy that y’ is a characteristic vector 

belonging to A For a problem of type (1), where A Ri (ar.t Ri‘ (6r,Gn) } <|(82.Gy 

is Hermitian and B is positive definite, every char- 

acteristic value is finite and real. Let < (6r,Gbr)) (n,Gn)} n,Gn 

M<d " It is an easy matter to verify that dr=7/(,Gn 

the unique normalized vector for which equality 

be the & distinet real characteristic roots of (1), and | holds between the first and last terms above. Hen 

let L,=L(d,) be the characteristic manifold belonging | this veetor is the desired solution of (8 Introducing 


a change in normalization for convenience, » 
termed the gradient of uw (with respect to @). 

We shall have occasion to use the gradient relatiy 
to a side condition 


to A), that is, linear subspace spanned by the char 
acteristic vectors belonging to \ Then anv two 
subspaces belonging to distinct \’s are B- and 
A-orthogonal, and have only the null vectors in 
common Further, every vector has a unique 


decomposition of the form 6r,2 0, fixed lf 
el Here we wish to solve (8) relative to (10). ¢ 
defined by the equation 
For problem (1) the important extremum principle is Ge=—t+hz—Ga+hz 
mim ula r B-orthogonal to L,. = ae where A is determined so that (¢,2 0. Thus 
> A« | 7 | 
d max wis), , B-orthogonal to L . as" L bie te: (G-'2+. 2) 
Then, in light of (10), (x. é6r.Ge As before 
n particular, follows that the maximizing vector is proportion: 
to ¢, and this vector is chosen to be the gradien | 
A Sul) <d r #0 More generally, with several independent si 
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conditions 


Od 0) 


is as the gradient the vector ¢ is obtained where 
Go=—t+h h 13 


with the A’s determined so that ¢ Is orthogonal to 


each of the 2’s 


Thus the A’s are th solutions of 


h.(G ho(G n 
hia hoiG ” 14 
in which the determinant of the éA’s is nonzero. by 
the positive definiteness of G~' and the independence 
ol the 's 
The change th gp when we pass from a vector 
#() to the vector 2 an will now be computed 
where @ is some real number and n(x) Is given by (9 
Assume that s is not characteristic, i. e nx0 
Direct calculation leads to 
n Gin 
ld vee an : fir. a), rxzv 7x0, 
J Ba 
LS 
where 
> 
cr a IS Ot 
Ae r l 16 
({ a rao 
with 
n, Bn Ri(r, Br n, Bn 
Dp any ¢ : 7 >s=—pul? ulz 
' n,Gn f J Ba r, Br) ; ' 


Our iteration procedure takes the following form 


\n mitial vector zo is given. Then the sequence 
J is determined by 

j j ' n n(a 17 

where the real number a, is to be specified at each 


In order to be sure that (17) determines a 
well defined sequence we must verify that r,+0 for 


step 


every t. To this end suppose that for a given j, 
#0; notice that from (9 
r,.Gr Gri n ) 
Hence by 17 
i (ra r,xa a’(n,.Gn >, IS 
since G is positive definite Hence ] x) smce 
ro #0, the sequence Is well defined 
From (15) we have 
” Gn 
Mid uid j a / a Ae) a 
J Ba 
19 


this equation holding whenever 7,0 


970822 —52 ‘ 


4. General Convergence Theorems 


established 
general 
in the section 
two effective 
For the present, 
we assume that the sequence has the property that 


and the next will be 
theorems under a certain 
sumption on the real sequence a 
following these two we 


In this 
convergence 


section 


tis 


shall deseribe 
wavs of meeting these conditions 


there exist real positive constants 6, and ¢ such that 


0 H~a,<hb., and 0<¢ lila whenever n, #0 (>) 


To simplify our discussion we 
the trivial case in which »,—0 for 
In this instance s, is a characteristic 


wish to dispose of 
first index 7 
vector belonging 
The 
are 
pro 


and J J | 
XI 
we 


value u 
results to be given in this, and the ne 
immediately verifiable Hence 
ceed on the basis that 


to the characteristic 
section 


shall 


now 


n, +0 . | 4 


In particular (19) holds for every 7 


Theorem \ Su p pose that the sequence \a satisfies 
20 Then u Tae sw oa decreasing SEQUENCE that 
CONVETGES to a characteristic alue h of ] Also 
lim (r,, G4 d>0O; in particular the lengths xr, are 


hounded and 
ulation pont of \s 


ro levery acculi- 


in Lin 
by (19 


since 


bounded away from 
s a characteristic vector 
To make the proof first that 
2) is decreasing 


we notice 


and that the sequence uu 
B and G are positive definite 
bounded from below by the minimum characteristic 


value X\,, it that it limit 


Since this sequence is 


follows has a eall it yp’ 


By (5 19 20) there is a positive constant ¢ such 
that 
n,.Gn 
ft é(ula ul J f : , 
r,.Ga 
Hence ¢ O° in fact 
i | 1) 
From (18) we derive 
J Ga J (ra Il l act 


It is well known that the product on the right con- 
verges if Jat does The latter condition holds by 


0 
20) and (21 This establishes the existence of the 
limit d and the asserted property of 4 


It now follows that »,->0 and hence, by (9 
Ag u Ba 0) 22 
Let y’ be any limit point of {7,) ; there exists at least 
one. Then 7’ is a non-null vector which, by (22 


satisfies 


Ay’ —p' By’ =0 
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Thus yp’ is a characteristic value and y’ belongs to 
Lip’ This completes the proof 

Theorem 2. Let 20) hold If the sequence I has 
an isolated accumulation point y, then (x,} converges to 
y Consequently if the characteristic root X of the 
preceding theorem is sim ple l ‘ aun L | then 
j converges to a characteristi ector 

Lat 7] bye at isolated accumulation poim*#l and let 
P be the set of remaiming accumulation points. Let 


S, and S, be open sets with disjoint closures with y 


in S, and Pin S There is an such that for 
r, hes mn the union of the two open sets Let d’>0 
be the greatest lower bound of |x for vu in § ~ 
rm Ss since 9 0 by the preceding proof, we 
mav, by (20) and (17), choose an such that 
“iy ! Md’ for Hence if s, is in S,, 
i then j is it S lt follows that for some - 
r, is mS, for all Thus P is null This estab 
lishes the first conclusion of the theorem 

Let \ be a simple root. From Theorem 1 every 
accumulation pomtl y must sautisty (y (iy d There 
are exactly two veetors in L(A) which satisfy this 
condition By the first part of the theorem ] 


must converge to one ol the mi This completes the 


prool 

Theorem 3. Let (20) hold, and let \ be the character 
wstic value of Theorem | Then there is a sequence of 
rectors iy n LOX) such that 


lim (a 7] i) 


For the proof we utilize the decomposition 


J 7] vy, in LOA) and B-orthogonal to L(A 
Then 
The cor 
further subse 
which LIA), by 


By the above decomposition, the cor 


does not eonverge to zero 


Now 


Some subsequence e 


SUP Pose 
converges to ze () 
responding subsequence J has n 
quence J converges to y mm 
Theorem | 
responding subsequence | y 


converges, necessal ily 


to a vector yy” in LIA Hence 
lim (2, —y)= yy" 


and 2B orthogonal to this sub 
This contradiction completes 


is both in L(A 
Hence 0) 
the proof 

It is worthy of notice that insofar as the iteration 
deseribed in this paper is to be used as a practicable 
numerical method for finding characteristic 
vector of (1 then the conclusion of the preceding 
theorem is as effective as the assertion that the se 
quence {2 actually converges. For the theorem 
that the will come, and remain, 
within an arbitrarily small distance of some char 
acteristic vector, this vector possibly varving with a 


Thus 


space 


some 


asserts sequence 


5. Convergence to the Least Characteristic 
Vector 


method is a gradient pro- 
it is to be expected that 


Because our iteration 
cedure which decreases ult), 


under 
and the matrix G@ the sequence pg 


appropriate hypotheses on the problem 
will conuvergt 


\,, the minimum characteristic value. We sh 
show under a rather strong assumption that su 
convergence will take place, and further, that 

sequence {4 will converge, whether or not \ 
simple 


In passing, we remark that although for defin 
ness the iteration so as to produce a decreasing 
quence yw, has been formulated a slight modificati 
in (20) produces an increasing sequence; the change 

« a b. and f,(a,) <e 0 
The results of the previous section hold in this ca 
and under the forthcoming additional hypothesis 
this section, convergence will take place to A,. tl 
the greatest characteristic value, and to a correspone 
ny characteristic vector 
Lemma 1 Suppose that 


AG'B=—BGaG"'A 2 


Ther proble m (1) and the proble m 


(ra vB 24 
have a common comple te set of characteristi ectors 
Yo, VY, Uu ith Yi, By, 6 q Kronecker de lta 

To prove this B=// is written with //7 positiv: 
Hermitian. Then (1) and (24), rv 


definite and 


spectively, are equivalent to //-'AH AD ane 
HGH where Hla It is easily verifis 
that the condition (23) is equivalent to the commuta 
tivity of the Hermitian matrices //~'AH and 
HGH It follows by standard theory that thes 


matrices are simultaneously reducible to diagona 
form by a unitary transformation: hence thev shar 


rT] complete ortho-normal set of characteristic vector 


, ; The desired vectors Y, are now 
given by Hy. 

Theorem 1 Assume that 2°00 and 23 hold f 
a gi vn initial rector Ia. let ym he the smallest inted 
»>yu=1,2 kh) for which xo is not B-orthogona 


Then 


the characteristic manifold L 


lim yu iis lim g yXxO with yin L 


We employ the basis of Lemma | to write 
J a y oe Yor lon YU» 
It is assumed that the basis has been ordered so that 
the first r, vectors span L,, the next ryo—r, Vectors 
span Z,, ete. (We take m=0 By multiply 


each vector of the basis by | we may assume that 


lo, >, P wl a. - i 


Let 


be the characteristic numbers of (1) correspondit 
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Thus 
Furthes 
number of 24 


vectors of the basis 
A, 7 T As, ete 
th 
corresponding to ¥y 


to the SUCCESSIVE 


characteristic 


Thus 


more let » be 


/ ry 0) 
the 
(; By | j }? a 2 / 
Using (17 i) 1) and the preceding equality 


it is found that 
with 


(Also 


Since 29 is by hypothesis B-orthogonal to each of the 
subspaces L L we have a 0 for q | 
2 a By 25 a 0 for every 

Hence every x, is b-orthogonal to the same subspaces 
and by the extremum principle of section 2 we have 


\ 

Now for each q=? | — consider the 
sequence | a Fon 0, a@,>0, with the strict 
inequality holding for at least one value of g. From 
25) it follows that the sequence is nonnegative and 
nondecreasing, since the term in braces is not less 
than 1 for the present range of g by (20), »,>0, and 
the last displaved inequality The sequence is also 
bounded, since |r, is, by Theorem 1 Thus 

lim a ¢,>0 @=9 l, } 27 


with at least one limit positive 


Let y’ be an arbitrary accumulation point of {2 
there is at least one Let y be a subsequence 

By Theorem 1, ’ is a characteristic 
In addition it must belong to L,,, for 


converging to y’ 
vector of l 
otherw ise 


a y., Ba y,, By’ 0, 

contrary to (27 It now follows that a;,-°0 for p 
outside the range of g, as in (27 Thus 
y’ lim } » > ey, #~O mL 

Since y’ was an arbitrary limit point we have the 


desired convergence of x, to a vector y in L,,. Finally 


lim lim (4s, “4a y, Ay) 
mM . = = X, 
: r,, Ba y, By) 
Corollary The condition 23 is satished if l 


G BR. or 2 G 
This is easily verified 


land AB=BA 


The case G=B=T] with A 


real svmmetric was studied in detail in the 


papel by the present authors referred to earlier 


greacel 


6. Construction of the Sequence {a 


i 


We shall 


this sequence so that condition 


methods of 
x1) 


describe two coustructing 


is satisfied 


Le mma 2 Let the eal number b. sat sf 
} 2 Lt G » 
{) ) %) 
N \, MB 
where X X s the spread of the characteristic values of 
| and the other quant fics are defined hy 5 Then 
there is a constant ¢ 0 such that for every r#O0 with 
nO we have 
fir a ( " for h a h 
Let J and a@ satisfy the required conditions, 
Assume further that a#0. Then we may write (16 
in the form 
fa , Ne 
Y és / r (99) 
r I rT Ba ayn J Ba 
Now 
Vib 
Sa Nb b 2 
mG 
where we have used (5) and the extremum property 


of A, and A Thus the numerator on the right side 
of (29 the number (2—4 Our 
proof will be complete if we can show that the corre 
bounded 
is sufficient 


exceeds positive 


denominator 1s 
By 5 il 
But 


sponding 
formly 


” J 


positive ull 


in zs and a to show 


is bounded this is an immediate conse- 
quence ot (9 

As a consequence of Lemma 2 we have the follow 
ing result 


Theore Mh 5 
he such that 


Let the SE ME NCE a 7 | real numbers 


O0< b a bh (30) 
for constants h h. with the latter as a 28 Then 
this sequence satisfies condition (20 

Our second method of prescribing a-al(z) stems 


from the idea of maximizing fla f(a, a) as a fune- 
tion of a, hence choosing a@ as a zero of f’la \ 
simple calculation leads to 


108 ria DS a l 
"(a ghd f ’ rxzu 7x0 (31 
] 2qa ra 
\ function all is now defined as follows (‘hoose 
an arbitrary fixed positive constant bh Let 
first zero of f(a) on 0 <a*< h, 
aad r£~O,n FO 
b,. if no such zero exists 
This function is computationally simple; its con- 
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struction involves only the solution of a quadraty 
equation and a comparison of numbers 

Theorem 6 
initial 
means of the 


kor aqiren constant h, i) and a aie 
determine the sequence a hy 


and the equation 


rector 2,40 


iteration formula l7 


(a F; itn, <0 
a 
b, { n 1) 
where alz is given hy (2 The th sequence 


satisfies condition (20 


For the proof several properties of the function 


32) are established The coefficients of the quad- 
ratic expression in the numerator in (31) are uni- 
formly bounded in a This is a consequence of (9 
ana >) It follows that its zeros ar uniformly 


bounded away from 0. Hence, there is a number 


hb 0 such that 


Thus the fulfilled Now 


fia) is a non-decreasing function on 0<a<a(s); for 


first condition of (20) ts 


by £00 2 and (32) its derivative ts non-negative 
on this interval Select 6 hb, so that 28) holds 
Then 

flata f(b cb 
uniformly im z Krom this we see that the second 


condition of (20) also holds The com 


plete 


prool Is 


7. Obtaining Further Characteristic Vectors. 


with 


propose to 


that a characteristic vector y’ 
characteristic value \’ is known. We 
show how the preceding iteration scheme may be 
modified so as to independent char- 
acteristic vector. The procedure will be to start 
with an initial vector 2 which is B-orthogonal to y’ 
and maintain this orthogonality at each step of the 
iteration 
Thus, let 


Suppose 


secure a new 


By’ 3 
and suppose we have a vector +#0 such that 
r, 0. We wish our next approximation 2 —a¢_to 
be orthogonal to 2. i. e.. we require 


( 0 4 


In order to select the direction ¢ in an optimal man- 
ner, according to section 3, it is chosen proportional 
to the solution of (8) with the side condition (10 


as im (33) We thereby determine ¢ by 11) and 
(12 Notice that by (11 

as er c, & r. Gi 0 oD 
using i4 and ri Now suppose ¢ fir) <0 
Then a straightforward calculation using the first 


equation of (35) shows that equations (15) and (16 


are valid with » everywhere replaced by ¢ 


The iteration formula (17) is now replaced by 


r J ag jb 
with the corresponding formula (18) valid by th 


Qur present 
property 


second equation of (35 


however has the additional 


sequen ( 


r,. By’ ‘) 7 


An examination of section 4 shows that with one 
necessary verification, to be remarked on soon, the 
three theorems of that section remain valid We may 
now add, however, that the characteristic accumula 
tion vectors y in Theorem | and 2 are 3B-orthogonal to 


y’, and that the vectors y, of Theorem 3 have the 
same property. The required verification is to estab 
lish the equivalence of tir) =—0 with ale 0 and 


c,-0O with »,--0. Here, of course, 9 is given by (11 
We shall prove only the second equivalence; this will 
That ¢,-0 when 


and (12 Suppose 


suggest the proof of the first 
7, ~O is immediate from (9), (11 


r.—0 We note first that (£,,y’ 0: this follows 
from (7), (37) and the fact that y’ is characteristic 
Hence by (11 

(Gey Hh (By yy’ 
where A, of (12) has the obvious meaning. Hence 
h,-O. It follows by (11) that £&,--0 and henes 


n, 0, as desired 
The constructions of section 6 remain valid under 


the present iteration (36 It is only necessary to 
verify the uniform boundedness of ¢ r<0 
(#0. This is an immediate consequence of (11 


(12) and (9 

To maintain the validity of Theorem 4 of section 
5 in the present context it is not necessary to modify 
the iteration procedure from (17) to (36 The 
earlier method is adequate. For, by (25), we se 
that if the initial vector 2 is B-orthogonal to y’, so 


is each vector x,, and hence the limit vector y 


If several characteristic vectors \ es \ / with 
characteristic values \‘, XA", are known then 
the iteration (36) is to be used with ¢ determined by 
(13) and (14), where By’, By’ The 


resulting sequence {4 will be B-orthogonal to the 
known characteristic vectors and hence all limit 
vectors will have this property. It may be easily 
verified that the preceding remarks concerning the 
validity of the results of the previous sections remain 
in force 

The iteration r—af determined by (13) and (14 
is theoretically equivalent to the following procedure 


First form the vector 2’ =r—an=x1—aG—1£, and then 
determine k,, ko, so that 
" r’ k G kG . sS 


Is h-orthogonal to i. y"’ It is easy to verify 
that 2’ rag However the procedure just cle 
scribed has the computational advantage that th 
vector rs, formed at each stage is accurately B 
orthogonal to the known characteristic vectors. I 
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the alternative procedure (36), B-orthogonality of 4 Ther 
mav be gradually lost through round-off errors 
although this mav be remedied by the additional 

work of occasional 3-orthogonalization 


S* dim LA dim H 


It would also be extremely convenient for the It is sufficient to make the proof for the equiva 
determination of the &’s in (38) (or the h’s in (13)) to lent problem (41), where ; L(A byand L Liv;) 
have (@G 0 for j#m. This may be achieved | [t may be assumed that the characteristic values so 
by successively orthogonalizing as additional charac- | ordered that if »—0 is a characteristic value. then 
teristic vectors are accumulated. Suppose y’’ has | ,,—0. Let A. be the space spanned by Z,, L Y 
been calculated with y’ known. Then with By’, | Any two spaces L, have only the null veetor in 
define by common It follows that the statement 

By’ +l 
with / chosen so that ( () Then the G-'- >» dim L dim H 42 
orthogonal set may be used in (38) (or (13 


Suppose now that a third independent characteristi is valid for j—2. Assume that (42) holds for j—m<k 


we shall show that it is then valid for 7=m-+-1 If 


vector y’”’ is caleulated Put 

/ , ] $2) were false for ; m 1, then there would exist 

, ; a vector y ~Oof L such that 
with /,, /., chosen so that G G 0): . . . I 
4 } } y ek 
this determination is simplified by the G~'-orthog- ; : 
onality of ; The new G~'-orthogonal set /—1.2 n 
may now be used in (38 The extension 

to more vectors Is clear with not all the terms on the right null From 


Ky vk 4, is obtaimed 
8. Problems Equivalent to (1) 
We leave now the calculation of characteristi 


vectors and values and raise the questien of when vl But I y } I y since I -. 
general problem of the type 


y vy Vin Y 
Cr=\Dr, with |C—XAD £0 in X, 39 
. so that 
is equivalent to one of type (1), that is, one with A 
Hermitian and B positive definite. For the moment ; 
we impose no additional conditions on the complex ( ) T ( yu, 0, 
matrices Cand D). Clearly (39) has at most n char 
acteristic reots (including the possible real value | since ; £0. But (42) holds for ;—m. Hence 
X ©), Where n is the order of the matrices each term on the left above must vanish At least 
Consider a second problem one y;, sav yj, is not null. Then »,.,—%, lo <m 
This contradiction completes the proof 
Ma = Sa +0 From the lemma it is clear that the condition 


By asserting that problems (39) and (40) are equiva- 
lent we mean that there is a nonsingular matrix K 
and one-to-one correspondence between the distinct 
characteristic values of (39) and those of (40) with | for (39) is equivalent to the assertion that every 


>> dim L n (43) 


the following property: if \’ of (39) corresponds to | vector may be written uniquely, apart from order of 
\’’ of (40), then y’ is a characteristic vector of (39 terms, as a sum of characteristic vectors belonging 
belonging to )’ if and only if y’’ Ky’ is a character- | ¢o distinct characteristic values 
istic vector of (40) belonging to \’’. If 6 is a number Our main theorem is the following 
such that |\C—dD) +0, then (39) is equivalent to Theorem & For a problem (39) each of the following 
C—bD)z \—b) Dr, and hence to conditions implies the ‘other 
1. Problem 39 is equivalent to one of 
Py ; : 7 . l 7 ; fy pe ] 
key vFy, with E=D, I ( bD. . j,’ F #0. i th, <a a oe 
+] real and (43) holds 
111. There exists a positive definite matrix P 
Notice that » © isnot a characteristic value of (41 such that CPD* is Hermitian 
Lemma 2. Let ) a 3 . k. be the dis- That II follows from Lisa consequence of the well- 
tinct characteristic tlues of (39) and let H be the Space known fact that for problem l condition IT is 
spanned by the S Paces L L valid To show that I] implies lil we construct the 
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nonsingular matrix ) whose columns are the n 
linearly independent by Lemma characteristi 
vectors of (41 Thus /) FY A, where A is a real 
diagonal matrix comprising the characteristic roots 


of (41). Henee A= )'F'EY is Hermitian, that is 


Hence 
EPF* — FPE* p=yy* 14 
From the definition of E and F in (41) is obtained 
DPC*—CPID* 15) 


as desired 
Finally, to show that IT] implies I we observe first 
that (45) implies (44). Problem (41) is equivalent to 


EPF* vFPF* Le! a 


which is of type (1). This completes the proof 
Corollary The following condition may be added to 
those of Theorem 8 
lV There exists a positive definite matrix P 
such that C*PD is Tlermitian 
For the proof we need only observe that Il holds 
for (39) if and only if it holds for 


C*%z= xD", 


and then apply IIL to this problem 
If our theory Is limited to real vectors ovet the 
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field of real numbers, then the following specializa 
tions occur The matrices s1 and B are to be taken 
real symmetric and the matrices C and D are to be 
taken real. The matrix P of III and IV is real! 
symmetric and the condition of reality in IT is 
superfluous 

If the matrix P of II] or LV is known then the 
transformation of (39) to (1) involves only matrix 
multiplications and hence is computationally feasible 
For example, in case of IIT with || #0 we write 


CPD*2—\DPD*:, r= PD* 
From the solutions z of this problem we obtain the 


solutions z of (39) by a direct matrix transformation 
If |\C’ #0, we write 


DPC*z=v»CPC*z, N s=/PC* 


In the case of IV with, say |) #0, we write 
D*PCr=dXD*PDer 

which is of type (1) with solutions exactly those of 

(39). In the case of either IIT or LV with |C D\=0, 

we first transform to (41) and then apply the preced- 


ing technique indicated for |D| #0 


Los ANGELES, September 11, 1950 

















Proposed Revision of the Conventional Method 
of Wave-Filter Design 


Paul J. Selgin 


\ modification of the conventional ‘‘classica method of electrical filter desig 
proposed rhe change consists of the use of new parameters frequeney numbers”, whicl 
permit the design to be made from specification of peak and cut-off attenuation rather tha 
idealized cut-off frequency, as in the conventional method 


lhe method is extended to cover the dissipative, as well as the 


Suitable design curves 


l. Introduction 


Che so-called conventional method of wave-filter 
design is not as powerful nor as flexible as more 
recent methods, based on the “lattice’’ form of filter 
section [1];' nevertheless, it still has usefulness be- 
cause of its relative simplicity, and is, in fact, used 
by practically all designers of equipment to whom 
the solution of a filter problem is an occasional 
necessity rather than a constant challenge 

Because simplicity is the main virtue of the con- 
ventional approach, one is led to ask whether the 
ultimate of simplicity and effectiveness has been 
reached in this particular domain, or otherwise. 
Upon examining the question, one finds that the 
conventional method, inception (Zobel’'s 
first paper on the subject [2] was published in 1923 
has remained essentially unchanged, the attention 
of network specialists having been diverted, by the 
work of Cauer [3] and Bode [4], to other methods 
of synthesis 

The chief purpose of this paper is to make the 
conventional image method better suited to the 
solution of practical filter problems, having regard 
to the form in which the design data must be formu- 
lated In practice the cut-off frequency as such is 
not given, but minimum and maximum attenuations 
are given together with the respective frequency 
ranges (see fig. 7 The present design method is 
geared to this formulation of the problem and does 
not require assumptions regarding the cut-off fre- 
quency value. Also, the method permits the mini- 
mizing of reflection losses at the critical frequencies 
by removing the arbitrary assumption that the 
“nominal” filter impedance must coincide with the 
load impedance 

The method is further refined by suitably selecting 
the parameters that determine the values to be used 
in the elements of each filter section. In the classical 
method the cut-off frequency is one such parameter, 
the cut-off attenuation another. In the present ap- 
proach a single function of frequency, the “frequency 
number,” is used throughout, and values of this 
function corresponding to the cut-off and peak- 
attenuation frequencies are used as the design param- 


eters 


since its 


and tabulations 


nondissipative or lossles 


are included, 


2. Basic Structures 


Space does not permit a full treatment of all filter 
types in this article. We will, therefore, limit the 
discussion to the three most frequently used types: 
high-pass, low-pass, and symmetrical band-pass 

The basic “ZL” structure, or “half section,” for 
these three classes is shown in figure 1, in two equiva- 
lent configurations, su‘table for the formation of ‘7”’ 
section and ‘x’ The so-called 
“7-connected” configuration will be given particular 
attention in the discussion 
along with the two alternative configurations, 
“frequency number’ that is, the function of 
quency that should be used, instead of the frequency 
itself, in all compucations concerning the filter 


sections, respectively 


For each class is given, 
the 
L 

fre- 


3. Dependence of the Transfer Constant on 
the Arms Ratio 


The ratio Z,/Z, between the series arm and shunt 
arm impedances of the L half section determines the 
attenuation and phase shift characteristics of the 
image-terminated T o1 joining 
two such half sections. (Conventional filter theory 
is predicated on the existence of image-termination 
or “image match” [5] between the various sections 
and to the remainder of the system.) The complex 
number a+ 78, where @ is the attenuation in nepers, 
8 the phase shift in radians, is, in general, the transfer 
constant @ of an image-connected 4-pole. According 
to this definition, @ may be shown to equal the com- 
plex ratio of output to input voltage (output to 
input current) of the 4-pole [5, p. 30] 

In the following, @ will always be associated to a 
complete 7 or m# section (it makes no difference 
which), while Z, and Z, will be the 
component half-sections 

The relationship between Z,/Z 


lowing 


r section obtained by 


the values for 


and @ is the fol- 


Z ” 
7 sinh =? (1) 
or, expressing é° 
Z Z 
y "Zz, ' Vz, " 
Z Z ‘ 
. = Va 
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URE | Basic “L"’ half-sections 


These expressions and the corresponding families of 
curves will be found standard textbook on 
filters [6]. 

The dependence of a and 8 on Z, Z,. for Z, and Z 
purely reactive, may be summarized in the diagram 
of figure 2, which is schematic and not suitable for 
purposes of evaluation 


im any 


4. Arms Ratio Expressions 


From what has been said it that the 
arms ratio Z,/Z, determines the transmission charac- 
teristics of the filter section At the moment we are 
interested only in the dissipation-less case; however 
since ultimately dissipation will have to be con 
sidered, it is convenient (in spite of the added alge- 
braic difficulty) to derive at the outset the arms ratio 


is apparent 





ta 





2,/Z> 
? 


Fieure 2. 


i1 , 











od > 
Zs Rr’ 
y Zp 
c 
/ 


Figure 3 


for the case of dissipative structures, differing fro! 
those of figure 1 in that the pure inductances at 
replaced by physical coils of given dissipation facto 
d=R/wl, approximately constant over the operatin 
range. The three classes of figure 1 will be cor 
sidered in turn 


4.1. High-pass 
(for T 





connection Th 
equivalent circuit, allowing for dissipation in the coi 


Consider the first form 


is shown in figure 3. The arms ratio is 





Z 
Z, 


At this point we discard as the variable denoting 


frequency and use instead the frequency numb 
appropriate for this class of filters, that is (fig. | 
the variable: 
l 
nm ! 
on 


terms of n by 


Equation 3 can be rewritten in 
noting that for zero dissipation, the infinite attenua 
tion frequency is given by 







by 



























Since the dissipation factor d reciprocal of the 


coil Q) is always a small number, it is permissible to 
mY clect 2? in the above proy ided we do not attempt 
to use the resulting expressions in the immediate 


neighborhood of the frequency of peak attenuation 
in which the denominator of (7 
finite if d is negle ted This is not a serous limita- 


tion to the design of dissipative filters, since the peak 


bee omes it- 


case 


values of attenuation are never critical, bye ing far in 
excess of design limits 


With this simplification, (7) goes mto the form 


The same expression goes for the second alternative 


form of high-pass filter section, which is in dual 
relationship to the first (7 
4.2. Low-Pass 
Avain consider the 7’ connected form (fig. 1 The 


equivalent eireuit for 
figure 4 


dissipative coils appears in 
arms ratio is 

9 
Discarding w in favor of the frequency number n, as 


before. we let 
Ww 10) 


Again, ih the dissipation less cause we have 


and at cut-off 


L i. - 
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so that of becomes 


As in the high-pass case, we neglect d? with 
the understanding that the resulting expressions will 


may 


not be accurate in the neighborhood of the peak 
attenuation frequency Hence 
} / 
l 
} } l jd 14 
/ 
} | 
| } 
/ 


The dual low-pass form, for # connection, vields the 


same value for Z,/Z 
4.3. Symmetrical Band-Pass 


The impedance of a series resonant arm, 


l 
Z R jJ@ L T 
Jw 
may conveniently be written 
yA / d iM 5 
where 
] 5 
j Wr 
1=R 
a 
\ 16 
a . = 
l 
\ LA 4 





Using this expression for the series combinations 
of the band-pass structure (fig. 5 
the form 


Zs u 2] 
R' np 7 
7 be 


pee 


its arms ratio takes 














Figure 5 





Z ‘ 
gz =re+ Mi le pw)" r"(d" w>| . 


Design requirements inherent to the band-pass filter 
impose a definite relationship upon the coefficients 
This is best understood by examing the behavior of 
the function 7,/7, with frequency in the nondissipa- 
tive case 


In the dissipation-less case 


; 0 when VW 0, al w 
Z 
Z, 


© when W’ orW’=0, 


l 
] ’ respectively 
nou ‘ 
y 4 


In general, the band-pass structure yields two pass- 
bands. To make these contiguous, we must meet 
the “condition of congruence” {8] that rW and 
L/W’ +-1/r’?W” vanish at a common frequency 
Since UO at w w we need 


If the filter is to have a symmetrical pass charac- 
teristic (on a logarithmic frequency scale), 
(19) 


leading directly to 
(20) 


If we add the condition that each series resonant 
arm has the same dissipation factor d, (17) simplifies 


to 
. ££ ; I ] - 
+ wid+gu lachrtacr | (21) 


As before, the denominator is to be rationalized and 
cc neglected, \ ielding 


r WOW’ ‘ er w+ Ww" 177 
rr’ ww oT iwwwsaw) wir 


We now define the frequency number for the sym- 
metrical band-pass case as 


($-2)-" 


Ww Ww 


a generic value of frequency; hence the value: 


(24) 


for either frequency of peak attenuation, by (19), an 
the value 


for either cut-off frequency (again because of sym 
metry 

Equation 22 may now be rewritten in terms of thi 
above values of n im place of Ww’, W”. By using 
the identities 


we find: 


Zz rn + ti 
2, yr’ i 

(27 
We have not yet imposed the condition that Z,/Z, 
must equal—1 (for d=0) at the cut-off frequencies 
corresponding to n=n There is no ambiguity on 
this point, because the other possible value that 
7,/Z, might take at the cut-off frequencies, which 
is zero, actually occurs at the midband frequency 
wo, the cut-off frequency common to two pass-bands, 
which merge into one by virtue of the “‘congruence 
condition. 

The cut-off relation is 


and the arms ratio becomes: 


jd Fnn, +2 
V7 Lv nN. 


5. Determination of the Section Numbers 
and of the Cut-Off Frequency 


The real part of the ratio Z,/Z, is independent of 
dand is given by the same expression for all three 
cases (8), (14), (29): 


Re (Z)—(Z"), 


where n is the “frequency number” appropriate to 
each case in turn. 
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We may now substitute this value in eq 2 and 


thereby obtain expressions suitable for the construc- 


tion of universal sets of curves of attenuation and 
phase shift, valid for all the filter sections shown in 
figure 1. To avoid imaginary arguments, three dis- 
tinct cases corresponding to the three ranges of fig 2 
will be considered 
l 
X 
v= 0.80.0 Y l } ! ! 
transmitted frequency range 
_ 
) A ) a ‘ | 
6=—2 sm ; 2 sin 3 
soa ne OO 
\ i yi 
> 
0.8—0 x 0 
ay ) ; i n n 
ee 
yi 
| 
9 } 9 sin} n »” 
a2 sinh 2 sinh 32 
\ \ \ /i ; 
yi / 
3 
0,8 x l 
ax, ft _s Ne NN, 
A 
i 
l 
9 | ) | a 2 
a » Cosi « COsD > 
VA ne? 
\ n n 
Note that in all cases n n The value of a or 


8 is seen to depend on two ratios, n,./n and n,/n 
For a given ratio n./n,, a or 8 may be plotted against 


n/n. These plots may be used to select values of 


60 


the ratio n./n_, which we will call the section number 
and designate by N 
\ 34) 
/ 
If more than one section ts used in a filter, a dif- 


ferent value of N will generally be used for each sec- 


tion As a rule there will be one section having 
\ 0 (the “prototype” of the literature and one 
section having 0.6 \ 0.7 the latter is usually 


divided into two half sections which serve as termi- 
for the filter 

Another purpose of the attenuation plots 
is the determination of the cut-off 
frequencies This is generally treated as a known 
requirement of the filter, but actually it is not. 
What is required of a filter may be summarized in 
figure 7, which is drawn for the high-pass case 
Two frequencies and two values of attenuation are 
the pass band is assigned by a ‘‘boundary” 


nations 
fig. 6) 


{re quency (or 


assigned 
value of frequency, fy, 
attenuation, a; the attenuation band, by a boundary 
frequency f, and a required attenuation, A The 
theoretical cut-off frequency generally some- 
where between the two boundaries. To determine 
its exact location, the section number for the entire 
filter must be tentatively established Suppose for 
example a three-section high-pass filter were to be 
on the requirements of 
tentatively make 


and a value of permissible 


lies 


based figure 7 Let us 


N=0 for section 1 

N 0.5 fo section 2 

N=0.7 for section 3 
Adding the attenuations plotted on figure 6 for 
the above values of N over a range of values of 
n./n, we find that the over-all attenuation goes 


the value of A 
n has the value 0.74. 


through the value 55 db 
attenuation) when n 


required 


We may 

















Ficure 6 Plot of a against n./n for various values 
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fp = 23 kc 
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ke 


Fieure 7 


therefore impose the condition n—n, for n./n=0.74, 
thus causing the attenuation to be exactly equal to 
A at the “boundary frequency” f, (fig. 7 Hence 
we write 


0.74, 


from which recalling 


high-pass filters 


the equiy alene en 


l 


Ae a 


20.9 ke: n 
0 i4 


The selection of the section number may be based 
on several criteria. If the problem is given in the 
terms of figure 7. the designe! will attempt to secure 
the minimum attenuation <A in the 
example 55 db) with the least number of sections 
Charts that permit the automatic selection of the 
optimum section number, and give at the same time 
the value n/n 
have been published by the writer 
of this criterion 

Another possible consideration affecting the selec- 
tion of the section number is the phase-shift over the 
transmitted band. It may be desirable to have the 
phase shift vary as linearly as possible with the 
frequency, to eliminate phase distortion. Gener- 
ally speaking, a certain “eut and try” 
is unavoidable when designing a filter by the con- 
ventional method. The involved may be 
minimized, however, by a favorable choice of param- 
which justifies the present treatment 


necessary 


for filters of three sections or less, 
i7] on the basis 


amount of 
labor 


eters 


6. Determination of the Maximum Permi 
sible Dissipation Factor 


values of Ne=n,/n, for the differ 
now determine for each value o 
in the dissi pat 


the 
sections, we 


Griven 
may 
(that is, for each frequency), 7,/Z, 
case (70) using one of the three equations 
(14) and (29) (fig. 10 may conveniently be used 
place of eq. 29 From the values of 7,/7,, us 
eq 2 or an appropriate chart (fig. 8), the true val 
of a and 8, may be obtained 

It is not hecessary, as a rule, to obtain these valu 
for more than one or two values of frequency 
check on the attenuation at the edge of the trar 
mitted band (frequency /,, fig. 7) is usually sufficie: 
since this is the critical value: the attenuation d 
to losses Is alwavs less toward the center of the pa 
band (in the ease of figure 7, toward infinite fi 
quency) than at the edges 

As an illustration, let us determine the maxim 


PaABLE 1 Sample computation of dissipation factor 


High-} filter 


values thus obtained for sev 
by interpolatior 


The over-all attenuation 
values of d may now be plotted to find 
value of d consistent with the maximum required attenua 
at f, (3db This is carried out in figure 9, giving as a re 
d—0.046 as the maximum allowable dissipation factor, « 


sponding to a “Q” of about 22 
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allowable dissipation (minimum () for the filter o 
figure 7, using the values of .V already selected 
Equation 8 is used to find Z,/Z, (hence a and 8 
for each of the three sections of the filter, at the 
boundary frequency f, (fig. 7) and assuming for d 
dissipation factor) a series of arbitrary values. 
Then by interpolation (fig. 9) one may arrive at a 
value of d just low enough to secure the desired 
maximum attenuation @ (fig. 7) at f,. The transi- 
tion from Z,/Z,, as computed from eq 8, to @ and 
8, is best accomplished by means of charts, such as 
that of figure 8. In this chart, as in table 1, the 


notation Z,/Z,=U+jV is used. The computation 


mav be conveniently carried out in tabular form, 
as shown in table |] 
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Midseries and Midshunt Impedances 


Having determined the frequency numbers n, (for 
the filter) and » for each filter section, as well as 
the maximum allowable dissipation factor, the de- 
signer still has the problem of “image matching” the 
filter sections to each other and to the terminations 
between which the filter is to be inserted. We must, 
therefore, reexamine the filter structures of 
figure 1 from the standpoint of their image imped- 
ances 

In this discussion the Z structure, or half-section, 
has been taken as the basic constituent of the filtet 
Such a structure has two distinct image impedances 
that is, those impedance values which, connected to 
the input and output of the L network, match, or 
simulate over the entire frequency range, the input 
and output impedances of the network, respectively 
(fig. 11 In filter terminology, the image imped- 
ance on the side of the series arm (to the left of fig 
11) is called “mid-series impedance’, 7,,; the other 
image impedance is the “mid-shunt impedance’, Zo, 
The essential difference between the two ty pes of L 
structures shown in figure 1 under the headings ‘‘for 
connection in 7” *, is this 
In the first type, the midseries impedance is inde 
pendent of the resonant frequency of the L-C arm 
so that, if such Z structures are connected together 
into 7 match 
one another at all frequencies ih having 
different peak attenuation frequencies 

On the other hand, the second type of L structure 
has midshunt impedance independent of the L-C 


basic 


and “for connection in x’ 


sections, these can be designed to 


spite of 


arm resonance, and must therefore be paired off 
into mw sections before assembling nto a matching 
filter 






































_° 


Figure 11 Figure 12 


To avoid confusion, let us limit our discussion to | Classes of filters 
the first, or T-connected, type of Z structure. The 
midseries impedance of these will now be obtained, 
for the nondissipative case 

The midseries impedance is the characteristic im- 

The mid imped the cl t { 
pedance of the 7 (fig. 12 Accordingly, if 7). and 
7,, are the open and shortcircuited impedances of 

» T, respective » have f -serie x 
the 7, respectively, we have for the mid-series | we may now substitute this expression into 3 
impedance 

obtaining 


z= 722 


We must now differentiate between high-, low 
and band-pass filters. For the basic high-pass stru: 


> 


(39) ture, we have (fig. 3 and eq 4 


We have already obtained an expression for the Z 
arms ratio in the nondissipative case, for all three i 
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; 














and also eq 0 


hence combining 


Z Vv ( 


Let us now define 


high-pass R= Oo 14 


and write the mid-series impedance in the form 


A plot of Z, against the frequency number n is 
shown diagrammatically in figure 13 
transmitted band) the radical must be assigned the 
imaginary value to fit the physical facts 
then the midseries impedance is a positive real, with 
maximum value Ry for »=0, which corresponds in 
all cases to the zero phase-shift or ‘‘midband’’ fre- 
queney Hence R \ hy (' 1s the midband impedance 

It will be noted that Z does not depend upon 
n.., hence is not affected by the peak attenuation 
frequency. High-pass sections differing in the loca- 
tion of the attenuation peak may therefore be image 
‘matched to another within the same filter, 
provided they share the same values of n,. (cut-off 
frequency number and R The above plot applies 
to all three classes of filters except for the sign of the 
imaginary part; eq 45 also applies, except that R 
has different values in terms of the circuit elements 
t and eq 40 


Fon n<n 


positive 


one 


For low-pass filters see hig 


a 
Z ol - 16 
n 
\ 
n 
which can be written 
: L n 
Z \ ( , \ n —] +7 
which, if we let 
L 
\ ; R iS 
becomes 
Zu=jR Va 1. 19 


In this case also the midseries impedance becomes a 
positive real in the range n<n,, provided we assign 
a negative imaginary value to the radical Outside 
this frequency range the impedance is a positive 
reactance, whose dependence on n is the same as in 
the case of high-pass filters, except for sign (fig. 13 

Finally, for symmetrical band-pass filters (tig. 5, 
eq 40 
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The midseries impedance will have the Same form 
as in the high-pass and low-pass cases, prov ided 


im W hic h cust 


R,. as before, is the midband The mid- 


series impedance follows the plot of figure 13 for 


mpedance 
n>n the sign of yn in eq 19 1s positive if the fre 
quency is above the transmitted band (Z,, inductive 
and negative if below 

In conclusion, we have found that in all three cases 
high, low, and band-pass) the midseries impedance 
is the same function of the frequency number; in 
therefore, invariant throughout the 
filter, whatever the location of the attenuation peaks 
in the individual sections This permits the filter to 
be image-connected throughout ; which is the assump 
“conventional filter 


each case it 1s, 


tion upon which the entire 
theory is predicated. Of course, this is only true 
provided FR), the midband impedance, and n,, the 


cut-off frequency number, are constant throughout 
the filter. This is ensured by determining these two 
parameters first and selecting element values so that 
eq 44, 48 and 51, for high, low and band-pass filters 
respectively, are satisfied 

We have discussed the determination of 7 but 
not that of R This is based on the terminating 
impedance of the filter. Always assuming that the 
‘T-connected”’ type of structure is used, and there- 
fore, that the midseries impedance is invariant from 
section to could put together the filter 
a three-section filter, for example) as indicated in 
figure 14. Each ZL half-section is shown separate: in 
practice, adjoining arms are “lumped”’ together 

This arrangement would offer to the 
the midseries im- 
that recurs at each internal junction if 
which we 


section, we 


outside 


terminations the same impedance 
pedance 
the terminations are equal resistances 








hr 


considerable “mismatching” 
would result, as indicated by the plot of figure 13; 
the midseries impedance plot is actually an are of 
circle within the transmitted range when frequency, 
rather than n, 


assume to be the case 


is taken as the abscissa 


8. Terminating Half-sections 


A significant improvement may result if we connect 
the filter as shown in figure 15. One of the sections, 
suitably selected, is “‘split’’ into its component half- 
sections and these are used to “terminate” the filter 
The filter now presents to its terminations (source 
and load) not but the midshunt im- 
pedance of the terminating section 

The frequency dependence of the midshunt im- 
pedance must now be examined. This will show how 
the terminating section should be selected and how 
the parameter Ry, equal to the midband value of 
both the midshunt and the midseries impedances, 
should be with reference to the 
termunations 


the midseries, 


assigned outside 


Zs Zs 





e —@ 


Figure 16 


is the characteristic 


We have, 


The midshunt impedance Z 
impedance of the wr section of figure 16 
therefore 


Z..= JZ, ¥. 


} pe + 
Z Le 
\ T 


(7) 4 


‘may write 


The midshunt impedance, like the midseries, is 
real in the transmission range (n<n,.) and takes the 
value Ry at midband (n=6 Its 
dependence, however, may vary from section to 
section according to the value of n, Plots of Z 
against n/n. in the transmitted band are given in 
figure 17. These show that for N=n./n,=0.7, o1 
thereabouts, Z), is substantially constant over a 
large part of the transmitted band. This is the 
value of N that should be in designing the 
terminating section. 


9. Determination of R, 


The use of properly designed terminating sectior 
reduces, but does not eliminate, reflection losses du 
to mismatching at the terminations [l, p. 352ff] 
These losses can be evaluated, but if the designer 
merely has to meet some such requirement as that 
shown in figure 7, this is not necessary. If the filt 
is designed to match its terminations exactly at 
f=f,, at this frequency there will be no added los 
due to mismatching. Above and below this fr 
quency there may be some added loss, but attenu: 


frequency 


used 
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tion falls off toward midband more rapidly than 
would be required to offset the increasing reflection 
loss 

The high-pass filter previously chosen as example, 
and based upon the requirements of figure 7, in- 
cludes a which N=0.7 
Choosing this as the terminating section, we go to 
the plots of figure 17 
to} f, 1s eq 35, table | 


section for section 3) 


The abscissa corresponding 


O.S2S 


Hence, in order to avoid mismatching, we must 


have, if Ry is the terminating resistance 
Z R 
at the point where the curve for .N 0.7 crosses the 
abscissa n/n, O.S2S 
But at this point 
- 0.985, 
hi 
hence, we must make 
R 
R - 1.015 R, 
0.985 
The difference between the optimum value of R 


and that of Rp is generally small, so that Ry, is often 


taken equal to Ry. But this simplification is not 


always permissible if accuracy is desired 


10. Computation of Filter Elements 


Krom the design parameters 7 and R , common to 


all the sections of a filter, and N=n.n,, varying 
from section to section, plus, in the case of sym- 
metrical band-pass filters, the midband angular fre 
quency @o, all the filter element values may be 


obtained by straightforward computation. Table 2 
gives the expressions for all element values of the 
three structures of the “7-connected”’ type, 
in terms of the parameters throughout this 
discussion, and with reference to equations appearing 
elsewhere in the discussion the dual 
relationship between ‘7-connected” ‘“w-connected” 
filters, it is easy to obtain the equivalent r-connected 
expressions: as always in a dual circuit relationship, 
capacitance is replaced by inductance, resistance by 
parallel, and vice-versa 


“basic” 


used 


Because of 


conductance, series by 
(1, p. 246ff]. 

The so-called “prototy pe” sections are considered 
here as limiting cases of the “basic” for 
N=0. In high- and low-pass sections, the element 
values are obtained in this limiting case by simply 
making V=0 in the expressions for the basic section 
One of the elements drops out, as indicated by the 
fact that the capacitor C’ in the high-pass section 
becomes infinite and the inductor L’ in the low-pass 
In the band-pass case, it is not 


section, 


section goes to zero 


convenient to let NV go to zero, as indeterminate ex 


pressions will result. Separate expressions for the 
prototype elements, which can be obtained quite 
simply by going back to basic considerations, are 


cviven 
Finally, it should be noted that the series arms of 
adjoining sections are combined into a single arm of 
the assembled filter Thus, 
add in the low-pass case and the Series capacitances 
go together according to the usual rules in the high 
In band-pass filters, it is convenient to 


the series inductances 


pass case 
add 


tions 


together values of r= 4Z/C for adjoining sec 

This will give the value of r for the “com 
As for the resonant frequeney, 
hence also for 


bined” series arm 
this is the same for all series arms, 
their combinations 


11. Conclusion 


The princi 
pal difference between the method outlined here and 


Comparison with conventional method 


the “classical’’ method lies in the choice of the 
design parameters for the simply derived filtet 
section. In the classical method these are: / 
cut-off frequency, and m derivation factor. — In 
the modified method they are: n., value of the 
frequency number at cut-off, and n,, value of the 


frequency number at the frequency of peak attenu- 
ation. Why is this an advantage? Because the 
frequency number is that particular variable that 
permits the use of a single function, or a single family 
of curves, to express the attenuation and phase-shift 
characteristics of all filter sections (low, high, and 
band-pass) over the frequency spectrum. With the 
aid of these universal curves or the corresponding 
functions it is simple to relate the design parameters 
as used in the present method, with performance 
If any other set of parameters is used, an_ inter- 
mediate step is required 

When using the classical method the designer is 
required to choose values of f, and m based upon his 
knowledge of filter characteristics When 
is usually given 


general 
the problem is an academic one / 
but this is not the case in practice, the requirements 
of the filter being summarized in some form similat 
to figure 7 After f. and m have been 
formulas are given for obtaining f,, the peak attenu 
ation frequency; on the basis of this and /, the filtes 
performance can be checked, to see whether or not 
the choice of f, and m was correct in the first place 
With the present method the choice of parameters 
and the verification of performance are one and the 
The parameters are selected by means of the 


selec ted. 


same 
universal charts, which show attenuation and phase 
shift at the various frequencies 

An additional improvement made possible by the 
present method is the computation of “‘minimum @.”’ 
The new parameters greatly simplify the analysis of 
dissipative filter sections. Formulas given im classi- 
cal texts for the dissipative value of Z,/Z, in the 
band-pass case, for example, are much more compli- 
cated than eq 29 and do not permit the drawing of a 
single family of curves as in figure 10 
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TABLE 2 ) ” m of element values 


Symmetrical band-pass 


4 


CAIVEN 


Series Arm: from eq 67 


I—N 


Lb. 
Ve"V = 


midband ang. freq. (giver 


l 
y Lh 


VR / , ) Compute 
’ SECTION “T” SECTION 


after adding r values of adjoining arms 
|}—-® 


SHuntT Arm: from eq 44 


neyn . l 
f 
rn 


from eq 34 


PROTOTYP! BAND-PASS 


given 


By a. 
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On the Derivation and Accuracy of Certain Formulas for 
Sample Sizes and Operating Characteristics of 


Nonsequential Sampling Procedures 
Uttam Chand* 


Formulas are given that are needed for the computation of number of observations and 
operating characteristics of single sample acceptance sampling procedures based on tests of 
statistical hypotheses Some of the same results may be obtained by reference to existing 
tables and curves located in widely scattered places rhe hypotheses considered concert 
the means and standard deviations of certain populations (that is, binomial Poissor orma 
where the test is made against a one-sided alternative The comparison of two means or 
two variances as well as the test of a single mean or variance is discussed Phe accuracy 
of the formulas is considered, and where approximations are involved, the results are com 
pared with existing tables, 


This paper develops and considers the question of 
the accuracy of certain formulas needed for the 
number of observations and operating characteristics 
of single-sample acceptance-sampling procedures 
based onstatistical tests of hy potheses For purposes 
of completeness the paper also contains formulas 
that are already known 


l. Introduction 


It is now knowledge among of 
modern statistical tools that the characteristics of a 
sampling plan must be specified in terms of the risks 
of rejecting good material (Type I error: producer's 
risk) and accepting poor material (Type II error 
consumer’s risk The problem of determining a 
sample size for a given probability of Type I error 
which will minimize the probability of Type II error, 
has been considered by several authors, in particular 
| to 4] This problem belongs to the category of 
that broad class of problems in the field of planned 
experimentation in which one is asked to provide 
adequate replication to detect treatment differences 
with desired amount of protection against taking 
wrong decisions. If for a certain sample size it is 
impossible to reduce simultaneously to small propor- 
tions the risks of two kinds it will be helpful to know 
this in advance 

Answers to most of the questions raised in this 
paper can be obtained from the existing published 
tables and curves The effectiveness of these 
formulas in relation to assumptions and approxima- 
tions that have been made in their derivation also 
has a theoretical interest We shall restrict ourselves 
to the consideration of certain parametric hypo- 
theses concerning means and standard deviations of 
certain populations, mainly against one-sided alter- 
hatives 


common users 


2. Preliminaries and Notation 


In the ensuing sections Hj, denotes the null 
hypothesis, 77, any one of a set of alternative hypo- 

Revision of a paper written during the summer of 1947 when the author was 

guest worker at the National Bureau of Standards The manuscript was 
actually revised while the author was teaching at Boston University 

Preset address, c/o I Sukhatme, Indian Council of Agricultural Re 

a New Dethi, India 

Figures in brackets indicate the literature references at the end of this paper 


491 


theses, a the probability of rejecting the null hypo- 
thesis #7; when true, and g the probability of 
accepting 7/7, when some alternative hypothesis, /7/,, 
is true. In connection with the hypotheses concern 
ing the means of certain populations in which the 
standard deviations are functionally related to the 
means and consequently unspecified, the reader will 
that the 


at once recognize acceptance-rejyection 
criterion A used for a statistic 7’ is not the best in 
the sense of the likelihood ratio test [5]. This 


difficulty, however, can be avoided by transformation 
of the original variables and has been indicated in the 
appropriate sections In connection with the two- 
sample problem the formulas assume equal sample 
sizes. These formulas can obviously be extended to 
cases in which it is desired to take unequal sample 
sizes of N, and N, that 
functionally related. 


are assumed in advance to be 


3. A General Formula Concerning Sample 
Size and Region of Rejection 


Let 2 rN be 


a normally distributed variable 


with mean yz, and standard deviation (sd) ao, 
f(w,) FUN) under H, and with mean yw and sd oa 
I ( pe FUN) under HH, vi iy), where FUN) is a certain 


and is independent 
as the critical region and 


function of the sample size 
of 0 We assume J A 


obtam 
y 4 
Ke 


Vos 
and 


l . 4 

‘ 21 ] 8, 
y2n JA 8 

where A, is the standardized normal deviate exceeded 
with probability « and actually 


A Ma , 
| 
ag I \ Ke 
A— ps ‘ , 9 
AU? I N A , As “7 
Solving (1) and (2) for FCN) and A we obtain 




















Kauollu K ey Tb 
A : ; (4) 
Kf (u Ka flu 
which may also be written as 
La io - 
l = ; >) 
A ~o Kao. 
Kau 0 Kou 0 : 
A ; > ’ (0) 
K,0,+ Kge 
for later convenience in certain simplifications. We 


note in passing that in particular cases of application 
f(u) will either be a function of u (cf. binomial and 
Poisson), or it will not depend on yz, in which case it 
will be a certain function of the population sd ¢ 
(ef. normal) or a pure number (cf. transformed 
binomial and Poisson) 


4. Tests Relating to the Parameter of the 
Binomial Distribution 


4.1. Single Binomial 


Consider a random sample of \ items drawn from 
an infinite population in which a proportion P of the 
items possesses a definite attribute B and let p be 
the fraction of the items possessing B in the sample 
Then E(p)=P and a(p VPA—P)/V¥N. Our 
hypotheses are #74,:P=P, and H., P PAP; Ps 
Using the normal approximation for the binomial 
variable p, (3) and (4) imply 


Nox (eva! P. K, 


hae VP-P ) a 


+ Ka P2y PA—P 
Kay P l P. 


—KsPivyPAl P 
Kay PA—P 


A 


While (7) determines directly the sample size that 
will (approximately) guarantee a specified a and 8, 
this may also be looked upon as providing the values 
of the probability of accepting //, for different values 
of P for given N and a. For example (7) yields 


\N(P,—P,)— Kav PA —P,;) 


K3~ . ) 


vy PA1—P, 





The inverse sine transformation @=2 sin~' yp 
where @ is measured in radians [16] avoids the diff- 
culty of the dependence of the standard deviation of 
of p on unknown /’, since @ is approximately normally 
distributed about 2 sin~' ,P with sd o(@)~y(1/N). 
In terms of the transformed quantities we obtain 


ate K+ Ks 
Ne,(5, Pp) 


: 5 (10) 
vl, 
Awa (“ee 


sim 


(11) 


y+ K, sin v=) 
K.+ Ks 
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sin 





Kg~2yN (sin VP, VP) 
We have derived (7) and (10) to illustrate the 
of the results given in section 3. The comparisi 
of the two formulas raises questions of quite a co 
plicated nature. We have so far not found a 
convenient yardstick with which to compare tl] 
accuracy. In the light of the fact that the eriti 
region @>A (11) has certain theoretical advanta 
against p >A (8) there is the temptation to reco 
mend (10 As the following example will indie: 
(a=.05, P;=0.1 and P,=0.2) the recommendati 
has nothing to do with the relative magnitud 

the values of N given by (7) and (10) 


; 0. 20 0. 10 0. O5 Oo. Ol 

7 68. 9 101. 2 132. 6 202. 8 

10 76.8 106. 3 134. 4 195. 8 
4.2. Comparison of Two Binomials 


, 


Consider two binomial]processes with P and / 
as the values of true proportions and let p and p’ 
be the observed proportions in a sample of N fron 


each process We have here ,:P*—P=0 and 
H,:P*—P>0. The test function Y= p*—p und 
H], has an unspecified variance VCY) 2|,P(1—P)IV/A 


where P is the common (unspecified) value of P an: 
P* under H Under 77, we have VL.Y)=P(i—P 
N+[P*(1—P*)|/N, where P*—P=£>0 say, wit! 
the value of & specified but not the values of P* and 
P. We are then faced with the problem of com 
paring two means having unspecified and unequal 
variances under /7/, and /7, 

The only satisfactory solution to the problem of 
the comparison of two binomial means—see [4 
chapter 7]|—is usually given in terms of the trans 
formed variables (section 4.1 Under the invers 
sine transformation sin™' , p*—sin~' yp is approx 
mately normally distributed (except when N is very 
small or the P’s are close to 0 or 1) with mean sin 


yP*—sin™' yP and variance (1/2N). We now us 
results of section 3 and obtain 
2 K.,+Ksz 4 
N=3 (= ): 
“\sin7', P*—sin=' P 


K3~ y2N (sin=', P*—sin-' P)—K.,. 15 


| 5. Tests Relating to the Parameter of the 


Poisson Distribution 
5.1. Single Poisson 


Let 7 denote the mean of a random sample of si 
N from a Poisson population with parameter 
Let Hi:m=m, and H,:m=m, To app! 


(M, >mM,). 




















the results of Section 3 we assume normality for 7 
with mean m and variance m/N and obtain 


yn kK. \ mokX 3 7 ° 
N ~( 16 
i Mi , 
{ Moy K miy MAX 3 | 
~ 
- sd ‘) 
ym Ka + y moKg 
> m mi \ m K : 
K3~ : 18) 


y 


We notice that under both //, and 77, the variance 
of 7 is known given N but differs for /7; and //, 
If we make use of the well known square-root trans- 
formation [17], the variance of the transformed 
variate is approximately independent of the unknown 
mean and is approximately equal to 1/4N. We now 
obtain 

v~i( K.+ Ks ) 19 


yi y 


= K.,, yn Ks yn 
K.+ Ks / 

and 
K3~2yN yi y iy) kK. (21 


For example for m 3.0, mz 1.0. a 05. B 10, 
formulas (16) and (19) yield N=29.3 and N=29.8, 
respectively. This is just an illustration; otherwise 
remarks made in connection with the single binomial 
in the last paragraph of section 4.1 apply here as 
well 

5.2. Comparison of Two Poissons 


Let M and M* be the parameters of two Poisson 


populations Let H,::M=M"* and H,:M< M*, 
e.g. 4 M* \ VU 7 >vV Consider two independent 
random samples of size N drawn one from each of 
the two populations. Let 7 and 7* be the corres- 


ponding sample means. We may regard the quan- 


tity ¥7*— 7 as approximately normally distributed 


with mean m= M*—, M and variance 1/2.N and 


consequently obtain 


( K,+ Kg ) rw 
2\,M*— JM fit 
hes K, »3) 
y2N 
and 
Kem y2N(, M*—~+M) —Ka. (24) 


6. Tests Concerning the Mean of the Normal 
Population 


6.1. Single Mean Test (¢ known) 


Let F denote the mean of a random sample of N 
observations from a normal population (4, 0”). We 
shall assume that ¢ is known from past experience 
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; 


kh io io 


Set wo—u,—Ac. The test function Z is normally 


Our hypotheses are //,;: w=, and /7, 


distributed with mean yz and sd a/yA Using 
results of Section 3 we obtain 


, K.+ Ks - 

N=(——*) 25) 

{ Kau Kyu 26) 

; K.+ Ks sii 
and 

K,=A vN—K., 27) 


6.2. Comparison of Two Normal Means (c's known) 


Let 7, F, be the sample means of two independent 
random samples of equal size N drawn one from each 
of the two normal populations 2,(4,,07), m(u.,¢ 
respectively 
from past experience 


We assume both o? and «2 to be known 
Our hy potheses are /7,: u Me 


and 7/75: uj<.. The test function 7,—7, is normally 
distributed | ve uy). 0 Co; N] and therefore. we 
obtain 
K.+- Kz : 
\ ( d ) 28 
K, lu a 
: 2) 
"hte . 
and 
kK vNd—K 30) 
where d - me 
yo 0 
If o CG; o” we obtam 
(Ket Ke | 
N 2( A ) 31) 
{ K, lu a > 
K.+k _ 
and 
; N . 
K3;=4 V2 K. 33 


6.3. Single Mean Test (¢ unknown) 


For notation see section 6.1 In this case the well 


known Student’s statistic t=yN(7—yp,)/s where 
s* is the unbiased estimate of o?, is used to accept or 
reject //,. Tables for determining the sample size 
for detecting a given value of A with preassigned a 
and 8 are given in [1l]-and operating characteristic 
curves for the same are given in [2]. To obtain a 
convenient formula for the sample size we proceed 
as follows: 

Consider P't>k}=P{((F—pu ks ¥N)]>0 
where & is a certain constant The quantity Z 
(7 vi (kes \ V) consists of two parts J vi is 
normally distributed (0, o7/N) under //,, and 


I 








(42 —,, o?/N) under //,; and for fixed N, ks/yN is 
a constant multiple of s where s is approximately 
normally distributed (C,e, Cyo/y2CN—1)) and C; 
and (, are certain constants less than 1.' It will 
be assumed for the purpose of this discussion that 
both C, and C, are equal to unity. The dual nature 
of these assumed approximations to the distribution 
of s should be noted. Therefore, 


E(Z\H,)=Z,=—-"=¢ 
\ \ 
” | I 
E(Z\H;)=Z:=(u2—» ° 
vA 
ose o , k o 
\ (Z) N T V 5) N j 


For Z>A as the critical region we obtain (see sec- 


tion 3) 
A ( ka ) 


K, 


; 1 - : 
JN Y T3(N—1) 
Eliminating A from the above two equations we get 


A l k 
K.+Ksz JN V Fo(N—1 #4) 
where A has been prevjously defined (section 6.1). 
This is a quadratic in N and could be solved for N 
if k were known. We now determine *& from the 
consideration that A is to be independent of o@ 
After simplification we obtain 


Ak, k 
4 (eK ow) 


The right-hand side will be independent of o if 
and only if the quantity in parentheses vanishes. 


‘ The general expressions for ¢ nd ¢ 
2 (4) 
c \ N ! : ~ Lf 
! ? 2\N ( 
N ! V 
I — } 
For smal! values of N, ¢ and (> take the following values 
N c ¢ 
2 0. 708 0. 842 
ANP 27 
‘ 921 4 
a4 we 
10 Q73 vs 
D5 an) Was 


Therefore, 
I KA \ N Qn 
K.+ Ka 


whence A=0 and is not only independent of ¢ but 
actually takes the originally intended value in th 
inequality P(Z>0). Substituting this value of 
in (34) we finally obtain 


\ b \ b ta 56 
2a 
where a (x a x.) 
and b ( l+a ( 1+ — )) 


Similarly 


Kz K.4 - ( N?—! + 
VN—1 2 


If we replace N—1 by N, the above two formulas 
reduce to somewhat simpler expressions of the form 


N (; } . Ki )+1 36a) 
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K; K.+4 A( N 


For any given A, a, and 8, values of N from (36 
are compared with the Neyman-Tokarska ® Tables 
[1] are given in table 1. These values will be found 
to be approximately the same. 

As pomted out in the previous paragraph (34 
was obtained under the assumption that & is un- 
known. The classical procedure employs k=t, (n 
where n=N—1 Therefore if the probability points 
of ¢ are not available, (35) furnishes an approxima- 
tion to such points Sa, We have considered the ac 
curacy of such & points in relation to ¢, in terms of 
Pi{t>z]=a(r). Values of t., k, a(t.), and a(k) are 
given in tables 2a and 2b for different N’s and fo: 
different a and Bg. Strictly speaking a(t,.)=a, but 
when ¢, to only three decimals is used, a(t.) may 
aiffer slightly from a@ as shown in tables 2a and 2b 

We notice that & values are in general conservative 
t, estimators and that the values of a(k) are con- 
sistently greater than the corresponding value of 
a(f.). In this sense a user of our formulas is likely 
to declare slightly too many significant results 
The danger, if it can be so called, is not very great, 
but it is still there. 

The question is asked: is it possible to eliminate 
this “danger” and still utilize formula (36) for N? 
It should be noticed that we can not utilize the 
available percentage points of ¢ in the derivation of 
(36). However, (34) can still be solved for N as- 


‘ Neyman-Tokarska’s p is equal to our Ay N. Since the standard table of 
probability points of ¢, [18] table 4, gives the twoe-tail probability points of t, our 


ta(n) corresponds to the entry given there for 2a and n degrees of freedom. 
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suming k to be known, say equal to t,(n ) resulting in 


N (a+ 1 y (a ] Zatz (n) 28 

Values of N from (38) corresponding to different 
t.(n), A, a and 8 are also given in table 1. Since the 
values of N as given by (36), (38) and the Neyman- 
Tokarska Tables are approximately the same, it is 
recommended that (38) be used in conjunction with 
the t-tables. This will save the labor of calculating 
k values 

Since any symmetric two-sided test may be re- 
garded as a combination of two one-sided tests. 
values of N and k may be obtained from the cor- 
responding single-tailed formulas by substituting 
K.,/2 for K, and t,/2 for t., in which case the “8” de- 
duced from the resulting value of Kz will over-estimate 


hod wr obtaining i 
8 
% V—4 . 
‘ ‘ { 
‘ t { { ‘ ‘ 
s % 
’ ‘ 48 { 
fe 
) { 24 
if ‘ =r) % % % 
' " 10 n 
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TABLE A Comparison o the accurac of the lifferent 
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able XXV f 


the true value of 8, the probability of accepting // 
when yw2.—y,—Ac, by the amount 1—’ where 8’ is 
given by (37) or (37a) with the term in A taken 
with a negative sign. Values of N obtained from 
8] and [9] and from the formulas (36) and (38) of this 
paper are given in table 3. The corresponding two- 
tailed values of ¢ and & are also given in table 3 


6.4. Comparison of Two Means: (common o¢ un- 
known) 


For notation refer to section 6.2. Consider two 
samples of equal size N. Our hypotheses are // 
li Me; /, py <_ Me Let s* denote the unbiased esti- 
mate of the common variance o The statisti 

Y Ps -_ 
t . which under //, has “Student’s” ¢-distribu 


“VN 
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Taste 3.* Compa ons among different w 
Two-sided single me 
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‘ 4. 184 
O19 
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s 8 " 4 _ 

Tal t 5 1 (9) and ¥ 4 


tion with 2(N 1) d.o-f. is used t 


Consider P(t >k ) P( 7. 7,—/ \ : . 0 )where 


Let 


oO a cept or reject H , 


we assume &, to be a certain unknown constant 
9 
Z=(fr.-—F kh ~~ 6 
. \ \ 
We that s is approximately normally dis- 


tributed with mean o and sd, o/ y4(V—1 There- 
fore Z is approximately normally distributed with 


assume 


E(Z iH) 


E(Z H, vi vi 


’ 2 I 
o(Z "VN 2N(N-1 


of Z. 


standard deviation 
obtain 


denotes 


where @(Z the 


Proceeding as in section 3.3 we 


A 9 I - 

K.t+Ks YN 2N(N-1 or | 

and 
A Kaod ~ h (40) 


it VA 


n) where ft, the one-sided 
1) d.o.f the equation (39) 


K, 
If we assume k,=t, is 
a-point of t for n—=2 (N 
vields 


(41) 


where a has been defined in section 6.3. The reason 
for such an illogical assumption about the knowl- 
edge of t,(n) before actually N is determined has 
been indicated in the previous section. If k is un- 
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known we determine fk, from the consideration that 


A is to be independent of ¢. The relation (40 
vields 
N K.A& 
R, — : 12 
V2 K.+kK, 
We substitute this value of k, in (39) aad obtain 
| 
N. 4 \ b? Sa 12 
2a 


where a has been prev iously defined (section 6.3) and 





b 2+a ( 1+ = ): 


For determining the characteristic we 


similarly obtain 


operating 


44 


We give in table 4 values of N as obtained from 
(43) and Tables of Neyman’® and Tokarska [1] 
While (42) provides approximate values of the per- 


TABLE 4. Comparisons among different methods for obtaining 
necessary sample sizes N* 
T'wo means one-sided test; common e« unknown: A 
a= 0 == ( 
Neyman-Tokarska|) Neyman-Tokarska 
43 43 
8 Mw) 2 2 ‘ 4 
B= .20 ‘ 4 t 7 
B=-.10 5 5 8 8 
B= .05 6 t ) u 
B= .01 s 4 12 12 
« N is the size of one of the two equal samples. 
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* For this particular case p as defined in [1] is equal to ay 
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centage points of t for 2U.N—1) d.o.f, the formula (43 
may be used advantageously by utilizing the avail- 
able percentage points of ¢t as acceptance-rejection 
criteria. Formulas obtained here for single-tailed 
comparisons can also be used for a two-tailed test 
by substituting appropriate two-tailed values of the 
quantities involved. If the two populations have 
unequal variances and their ratio is known we still 
can construct a test function similar to ¢ and use the 
formulas For the the test 
10). 


above structure of 


function see 


6.5. Comparison of Sample Sizes for Known and 
Unknown o 


It need hardly be emphasized that in situations in 
which large sample sizes are required, the normal 
test and 6.2) and f-test (sections 6.3 
and 6.4) will both vield approximately the same N 
see for example table 1] To determine the rela- 


sections 6.1 


tions between the two, let N,, be the number of 
observations required when o is known and N, the 
corresponding number when o@ is not known. After 


some simplification the (gllowing asymptotic rela- 


tions between N,, and N, are obtained For the 
single-mean test 25) and (36a vield 
‘ Kk? F 
A ~Nu (1 +55 )+1 15 
and for the two-means test, (31) and (43) vield 
Kk? ’ 
Ni~Nn(1+3y )+1 (46) 


7. Tests Concerning Variances of Normal 
Populations 


Tests relating to population variances fall into two 
well defined categories In one case we assume that 
the variability of a certain product is known and it 
is desired to find out whether a new product is more 
variable than this. In the other case we are asked 
to choose between two products on the basis of their 


variability which is unknown. We discuss these 
situations in the following sections. 
7.1 Single Variance Test 
Consider a random sample of size V n+1) from 
a normal population (m, o Let s* be the unbiased 


sample estimate of o*. Our hypotheses are //;: 
o’=o, and H,: o?=doi(A>1 For a given level of 
significance a, if y*’=ns*/o2 >X%2, we reject H, and 
conclude that o > 

Let A(a.8.n)’ denote the value of o?/e2 for which 
the probability of the decision o?=o§ equals 8 when 
the test is conducted at the a@ level of significance 
with nd.o.f. It can be shown [2], [4] that the proba- 
bility error of the second kind is exactly 8 if \(a@,8,n) 
x2 (n)/x?_5(n) If we are testing /7;: o? =o? against 
H,: o&=rA8 (A<1) have A(a,B,n fs ~ef% 


Our \(a,8,n) is « 


we 


of [4 


quivalent to p(as 
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Curves for the operating characteristics of such 
testing procedures are given In [2] and [4]. Eisenhart 
[4] has also given extensive tables for \(a@,8,n 

The problem of determining a direct relation 
between large n and any given set of values of A, a, B 


was first considered by Wallis [1, footnote of p 
278]. Assuming normality of s (see section 6.3) and 
applving the results of Section 3 we obtain 
1 {K.+Ks yi ‘ 
Ti a 47) 
2 ya—] 
ny ‘ f 
y ae ¢ K, A 18) 
K+ y\Kz 
Ka 2n 
Ma, 8, n)=( = ) 19) 
\ 27 Ks 


To compare the accuracy of this formula with the 
Tables [4] consider the following situation: if a deci- 
is a serious error from the practical view- 
point when o l and it is desired to keep the 
risk of such an error below .05 when the test is con- 
ducted at the 5-percent level of significance, how 
many d.o.f. will be needed for s*?) The formula (47) 
gives n=33.8 and from [4] we find that 34 d.o.f. are 
needed Table 5 presents the calculated values of 
n from (47) and table 6 presents the calculated 
values of \ from (49). For a comparative discussion 
on the use of these formulas in relation to others see 
next section 


sion @ a 
500 @ 


7.2. An Alternative Formula for the Single Variance 
Test Based on the Distribution of Log s* 


As pointed out by Bartlett and Kendall [11] the 
distribution of log s* depends on o* only through the 
term o° in its expected value Consequently the 
choice of the critical region based on the distribution 
In this 


some 


of log s in place of s has obvious advantages 
section we explore the possibility of using 
formulas based on the distribution of log s 


The cumulant function K(f) of log s is given [11] 


)—tog 1 (1) 


which vields the following expressions for the first 
two cumulants 


results of Section 3 





applied to log s vields 


0 


log n 
Kol ni t ov 
K+ Kg 

{ K, log kK, log \ \ ( v( j ) l } ) 

= = Y — lov 

Ki+kh 2 2 2 

= 

K+ i 

A apo 7 ‘ . 2 
These formulas assume no other approximation 


except that of normality and in that sense may be 
regarded exact relative to (47 which 
dual type of approximation for the distribution of s 


assumes 


see section 6.5 While the accuracy of (50) and 
(52 does not depend upon the accuracy with which 
we estimate x« it does depend upon the compli- 


; / . 
cated expression ¥ ( = ) For a given a, 8, and X\ the 
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asymptotic expansion of ¥ ( = ) 12). 


As a first approximation if we set m~1 


we obtain 


K+ kh 
log 
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low X K,+ Kg 4 


Vin 


In tables 5a and 5b values of » 
calculated from (47) and 
of 8, and \ are compared 
values of \ B, 
approximations based 


from [4], and value 
for different value 
In table we giv 
to reemphasize the nature 
the distribution of 
In this connection for the application of 
values of xo(n) were taken from |11] 

It appears that for the customary values of th 
probabilities of errors of the two kinds a=8=.05 

8=.01, formulas (47) and provide very 
approximations to 7 for 
If the percentage points gf the 
available, (47) is preferable because it is easier to 
compute n from (47) than from Even for such 
a small value as n=5, (47) errs on the safe side in 
this sense that it gives (at least for a=8) a sampk 
size which will be always sufficient to detect this 
difference. The formula also shares this prop 
erty with (47 In the absence of the percentag: 
points of the x*-distribution it perhaps ought to b 
emphasized that on comparison of the critical regions 
for s and log s (ef. (48) and (51)) there is not much 
basis for choice The choice of the critical region 
based on the distribution of log s has certain theo 
retical advantages, but the computation of the criti 
cal region is somewhat more complicated since it 
involves the approximation of ¥(n/2 

The effectiveness of formulas (47) and varies 
when @ and 6 are not equal. It appears (see tabl 
5 and table 6) that for 8 it is safer to use (53 
because it is always likely to err on the safe side in 
the sense of the previous paragraph. However, if 
8<a@ it appears that it is safer to use (47 
53) is likely to give a value of n, which will actually 
be less than the necessary sample size 


a, a ( 


a i oO} 


S ane 
9 


log S » 


and a 
good small sample sizes 
x*-distribution ar 
53 


i 


a he) 


becaust 


7.3. Comparison of Two Population Variances 
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Let and denote the two 


normal populations and let s7 


o; o> variances 


and be 
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and calculated values of X (a. 8. n 
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pendent sample estimates based on n, and n. d.of 
respectively Our hypotheses are //,: oj=032 and 
H 0 o The statisti } s7/85 1S used to accept 
or reject H For a given level of significance a, 
if F< F, we accept /7, and accept 7/7, when F >F 
Let @ denote the true value of oa7/a It has been 
shown in [2] and [4] that the probability of an error 
of the second kind will be exactly 8 if 


The operating characteristics (¢, 8) for one-sided 
alternatives have been given in [2] and [4 We shall 
develop here certain approximate formulas for @ 
in terms O1 a, 6, 7 and 

By taking one-half the logarithm of (55) we obtain 


- log © 3\ 7 } 56 


where denotes the 
-distribution., 


v-probability point ol Fisher's 
For purposes ol approximation, 56 


x 


has decided advantage over (55) in that it enables 
us to make better use of the Cornish-Fisher approxi- 
mation [13] since the distribution approaches 
normality relatively faster than the F-distribution 
We shall employ the following approximations [14] 
for the cumulants of 


Let and denote the mathematical ex- 
pectation of under 77, and H respectively 
Consequently we have 


By applying the results of section 3 to 
that 


it is found 


l l log > ” 
2(;, s ) (UK. K) a 


which may be rewritten as 


| 2(Ka+Ks 
io” @ Os) 


yh 


WW here h 


there results 


Directly applying to (56) the Cornish-Fishe1 ap- 
proximation 13] for the 


percentage pomts of the 
-<distribution in conjunction with Cochran's A cor- 
rections [15 we obtain (58) and the following two 


expressions for log @ 


) (A kK? 





\/ N \ 
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where A, Is given by 
€ 7) 50 oF O5 Ol 
¥ 5575S 5000 5758 Q509 1. 4020 
Formulas (58 59), and (60) are not changed when 
n, and no are interchanged Formulas (60) and (61 
are identical when either a= 8 or 7 no, or both 
Parte 7 Com} , j 61) and 
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the basis of several computations 
(not given here) that formula (58) is likely to give 
values of @ that are much lower than its tabulated 
values [4], and consequently the sample sizes given 
by it will fall below the minimum desired. For 
a=B8 formula (61), then equivalent to (60), gives 
values of @ which are much closer to its tabulated 
This is also true of (61) for B>a 
then from (60 


It appears on 


values (table 7) 
If in addition to a=, n n } 


n—raty (Ket Ke) 


log @ 


We have not found any formula which will give 
an approximately correct 
freedom as small as, 


degrees of 


The question 


answer for 
sav, No=5 


of finding appropriate sample sizes for n,#n,. cannot 
ordinarily be answered without the help of tables 


In such a situation, however, the experimentalist 
has no choice in the determination of n, and n.. If 
it is decided in advance to maintain a certain ratio 
and n, (this appears to be more often 
practice) a formula would be 
than the existing tables 

formula (61) seems to be very complicated to use 
we recommend the use of (59 Values of ¢ as given 
by this formula are given in table 8. It appears 
that for a=8=.05, formula (59) will always give 
sufficient sample sizes; but for a=8=.01, it will 
give values slightly less than actually needed 
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This paper could not have been written without 
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Asymmetries of Zeeman Patterns and s-Values 


for Neutral Manganese’ 


Miguel A. Catalan ° 


Spectrograms of manganese made at the Massachusetts Institute of Technology with 
fields in excess of 84,000 oersteds show many lines that exhibit various degrees of distortion 
n both the positions and the intensities of the magnetic components. The interpretation of 
these asymmetric patterns has beet made by the approximate theory of the partial Paschen 
Back effect Ihe 7 values that have been derived for several energy levels ot Mr I are 
found to conform, in most cases, with those required for LS coupling \ few exc ptions to 
this rule have been considered in some detail General tables have been computed, which 
will permit explanation of the distortions in other spectra 


l. Introduction 


Ever since the discovery of the splitting of spectral 
lines by a magnetic field (Zeeman effect) the ob- 
servers have called attention to the fact that some of 
the patterns are asymmetrical, the distances be- 
tween their magnetic components being not rigor- 


ously equal. These asymmetries make it difficult to 


compute the g-factors that govern the splitting. 
Two early examples of asymmetric patterns are 
given by Martinez Risco [1],° and by Back and 
Landé [2] 

At that time it was impossible to determine 
which of these different separations between mag- 


netic components ought to be used as coefficients m 
the equations for obtaining the g-values of the 
atomic levels. Back adopted an empirical proce- 
dure, which is described in his paper on manganese 
[3]. This procedure has since been much used be- 
cause In Many cases it gives values in complete agree- 
ment with the LS-coupling values 

In recent years some authors, who have had the 
opportunity of measuring the beautiful plates made 
at the Massachusetts Institute of Technology by 
Harrison and his collaborators using magnetic 
fields of about 85,000 oersteds, have again called 
attention to the fact that some of the patterns are 
very asymmetrical [4]. 

The author of this paper, who is at present work- 
ing on the structure of the manganese spectrum, 
has measured some excellent plates made at MIT in 
1939. The measurements indicate that many of 
the patterns are more asymmetrical, and 
that these asymmetries affect the g-values by as 
much as a few percent depending on the rules 
adopted to calculate the average separation between 
components. On these plates many of the most 
interesting asymmetric lines appear too weak to be 
measured accurately, some of the faint components 
being absent in many It was, therefore, 
decided to obtain new, long-exposure spectrograms 
of the manganese spectrum in the magnetic field. 


or less 


Cases. 





his research was carried out during the author's visits to the United States 
1048-49 and 1950-51. It was made possible by grants from the American Academy 
of Arts and Sciences, the American Philosophical Society, Princeton | niversity, 
and Massachusetts Institute of Technology, and also through the cooperation of 
the National Bureau of Standards 
? Present address, University of Madrid, Madrid, Spain 


* Figures in brackets indicate the literature references at the end of this paper 


The large Bitter magnet of the MIT, which had 
not been in operation since the war, was again put 
to work (not, however, without considerable trouble 
by J. C. van den Bosch and the writer with the help 
of G. R. Harrison’s assistants. A beautiful set of 
plates in which the interesting asymmetric patterns 
appear strongly was then secured. The procedure 
emploved in making the spectrograms is that de- 
scribed by Harrison and his collaborators [5]. 

The electrodes for these plates were prepared by 
compressing powdered manganese with silver dust 
The manganese was prepared by an electrolytic 
process and presented to us by W F, Meggers. It 
was extremely pure. 

The accurate measurements of these 
grams by the author have permitted him to ascer- 
tain the cause of the observed asymmetries and at 
the same time to deduce the best procedure to elimi- 
nate the distortion in the computation of g-values 
These distortions in the magnetic patterns appear in 
both the positions and the intensities of the compo- 
nents. In the following we shall describe first the 
asymmetries of the positions and later we shall 
‘consider the perturbations in the intensities. 


spectro- 


2. Experimental Asymmetrical Patterns 
The pattern for the line a°D, 2*—, 24485.25 
em~' of Mn, at 4082.945 A, obtained experimen- 
tally shows unequally spaced components. Irregular 
intervals between adjacent components are found 
both in the x- and e-components. This is an example 
of an asymmetrical pattern, which contrasts with the 
symmetrical ones generally described. In table | 
the theoretical values of the wave numbers of the 
components are compared with the observed values 
[6] for this line. The experimental intervals, although 
irregular, show curious regularities. Any cause to 
which the asymmetry of the pattern may reasonably 
be ascribed will not alter the validity of the combi 
nation principle. This means that the displacements 
observed in the magnetic lines will be due to dis 
placements of the magnetic levels involved. A per 
turbed level will show the same displacement in al 
lines that originate in it. Double-array tables will 
then be suitable to represent an asymmetrical pat 
tern, and in such tables the magnetic levels wil 
appear with their displacements. 
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Figure 1 indicates the positions of perturbed and 
unperturbed components in an imaginary example. 

The wave numbers in table 2 are experimental 
values that give the adopted values for the displaced 
levels from which the wave numbers of the compo- 
nents have been recalculated. The decimal parts 
of the resulting values are listed in table 2 under 
the corresponding observed value for comparison 
The close agreement between calculated and ob- 
served values in table 2 proves the validity of the 
combination principle when applied to an asymmetri- 
cal pattern and hence the possibility of computing 
values for the displaced levels with the observed 
wave numbers 

The observed intervals between the consecutive 
magnetic levels are thus unequal, being 7.09, 7.39, 
and 7.67 for a°®D,,., and 6.43, 6.50, 6.54, 6.62, and 
6.67 for e “Di. The differences, in both cases, 
cannot be ascribed to experimental errors but are 
clearly due to an asymmetry. 

At the bottom of table 2 the pattern for another 
line a °D, 2"D>., is included. This Mn 1 line has 
the level a*D,,, in common with the other. From 


24477. 18 
ies 24478. 32 
< 24474. 2 
ie 24479. 8S 
24483. 83 
“= 244584. SS 
re 24485. 72 
. 4486. 32 
4400. 44 
ie 24491. 49 
~ 24402. 2 ; 


the observed components one may calculate the 
perturbed magnetic levels as before. The resulting 
values for the magnetic levels belonging to a*D 

are exactly the same as those for the line 24485.25 
and hence show the same perturbed intervals 7.09, 
7.39, 7.67 Many other examples could be presented 
to prove the existence of asymmetrical patterns 
All of them show the same singularities. The tr- 
regularities in the displacements are exactly the same 
for all lines having their origins in the level involved. 


3. Asymmetrical Pattern in Theory 


In 1913, while observing the Zeeman effect of the 
line 6708 A of lithium, Paschen and Back found that 
the observed pattern was in complete contradiction 
with the Preston rule It is a close doublet, and 
there is a sort of interaction between the effects of the 
two lines forming it. Similar results were obtained 
for other close doublets and triplets, although lines 
that are not organically connected do not show this 
effect, even though their wavelengths may be very 
close to each other This interaction was called the 
Paschen-Back effect In some cases the observed 
Zeeman patterns for very close multiplets exhibit a 
sort of distortion due to partial Paschen-Back effect, 
the field in these cases not being strong enough for a 
complete interaction but only sufficient to alter the 
normal pattern. 

The theory of the Paschen-Back effect has been 
developed in detail by different authors Kiess and 
Shortley [7] have recently considered the distorted 
patterns in the oxygen and nitrogen spectra, and have 
found a quantitative theoretical explanation of them, 
The multiplets 3s °S°—3p'*P and 3s8°S°—3p°P of 
O1 at 8446 and 7771 A with the narrow intervals 
0.5, 0.7 and 2.0, 3.7 ¢m show a great distortion. 
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Tante 2 Observed perturbed pattern ne 4 85.25 and 24630.08 of Mint 
P ‘ ml 1 1738. 48 
“ul La f Ma V vi 
17557.53 09 17564.62 a9 17572.01 fi 17579.68 
2 $2037.41 24479.88 
SS 
$2043.84 24486.32 i] 24479. 21 
7 22 
‘ ‘; 
0 $2050.34 24402.84 24485.72 ( 24478.32 
1 72 35 
DD f 
0 $2056.88 24492.23 4484.88 0 24477.18 
26 87 20 
‘ ‘ 
| $2063.50 24491.49 He 24483.83 
wy 82 
2 $2070.17 24490.49 
at] 
0 $2102.45 24634.93 OY 24627.84 24620.42 
Q? 83 14 
‘Dp 13.1¢ 
0 $2205.55 24640.91 24633.54 f 24625.87 
93 54 87 
The multiplet of Ni, 38 'P 3p DS” at 7423 to 7468 A, 
with greater intervals, namely 46.7, 33.8, shows J—L+S8)(J+L—S)(L+8+1+J/)(L+S8+1 My 7 


much smaller distortion. The results of Kiess and 
Shortley for N 1, in a multiplet with the relatively 
wide separations, 46.7 and 33.8, indicate clearly that 
multiplets with relatively great intervals, in 
elements, could be expected to show also measurable 
Paschen-Back interaction. This conclusion 
gested to the writer the likelihood that the observed 
distortions in Mn are due to a similar effect In 
order to prove this hypothesis the structure of the 
levels so affected was considered In some detail 
According to the second order perturbation theory 
17, p. 204] magnetic levels having the same magnetic 
numbers 7, and belonging to spectroscopic levels 
which differ by one unit in J-value, and to the same 
term, repel each other by an amount, ¢€ equal to 
€ l ’ (1) 


0 


sug- 


all quantities being expressed in Lorentz units. In 
this formula 6 is the distance between the two 
magnetic levels under consideration. Within the 
accuracy of eq | 6 represents either the perturbed or 
the unperturbed distance. The factor J is the inter- 
action element; its values depend only on the partic- 
ular values of S, L, J, and M for the levels involved. 
Equation 2 gives the dependence of J on the quantum 
numbers in Lorentz units, LU 


other 


Here the value J is the larger of the two J-values 
involved The energy perturbation is to be applied 
to each of the states always in the repulsive direction 
If more than two states interact, the perturbations of 
each pair may be supposed to operate independently 

The writer has computed with formula (2) the 
values of J? and J for the different The 
resulting values are given in table 3 with the /-values 
above and the J*-values below. In general the 
values of the interaction element increase when eithet 
S or L increases and decrease when either J or M 
increases Hence the greatest interactions are ex- 
pected in terms of high multiplicities and high L- 
values, and, within a term, the greatest interactions 
may be expected in the levels with smaller J- and 
M-values. <A very important conclusion is that the 
sign of \/ does not affect the J-values because 
appears in eq 2 as a square. The same J factor 
corresponds to level (+./) as to (—M) 

Because the measured distances bet weeen 
magnetic cemponents and, therefore, the level- 
separations are expressed in wave numbers, cm™', it 
is desirable to express the values of J in the same 
units. This is done by multiplication by the 
Lorentz unit (LU), which is a=4.669X10~ Hem 
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4 4 

975 0. 737 

O51 0. 543 

771 0. 668 0. 498 

505 0. 446 0. 248 

059 

122 

183 0. 8O4 

100 0. 800 

007 0. 850 0. 708 

203 0. 903 0. 501 

843 0. 769 0. 659 

710 0. 592 0. 434 
4 

133 

188 

193 

243 

533 0. 400 

284 0. 160 

500 

250 











U pper figu are J-valu nl wer I ‘ LI 
™ 
1. 265 1. 005 
F, and * 1. 600 1. 200 
1. 242 1171 0. 926 
, 242 
F; and *I 1. 543 1. 371 0. 857 
F, and 54 1 ool 1. O56 0. 045 0. 722 
, Are 1190 1116 0. 8403 0. 521 
F. and 55 0. 816 0. 805 0. 748 0. 653 0. 490 
wens ‘ 0. 667 0. 640 0. 560 0. 427 0. 240 
CG. and *C: 1. 195 1. 127 0 891 
1. 429 1. 270 0. 704 
1, 254 1. 214 1, O86 0. 820 
G, and *G ot 
. 1. 571 1. 473 1. 179 0). OSS 
( d + 1. 101 1. 105 1. 034 0. 903 0). 677 
aint ' 1. 273 |. 222 1 069 0. 815 0. 458 
( nd 0. 853 0. 841 0. 804 0. 739 0, 636 0. 471 
te arn ' 0. 727 0. 707 0. 647 0. 546 0. 404 0. 222 
ee , 0. 756 0. 713 0. 563 
P; and "1 0. 571 0. 508 0. 318 
P, and ?P 0. 655 0. 634 0. 567 0. 433 
_ 0. 429 0. 402 0. 321 0. 188 
1. 265 1. 005 
D, and *D 1. 600 1, 200 
1. 242 1. 171 0. 026 
), al ) r 
Ds ana ‘I 1 543 1 371 0. 857 
1. Oo1 1. 056 0. 045 0. 722 
)», anc ) 
D, and *1 1,191 1. 116 0. 893 0. 521 
D. and 7D 0. 816 0. 800 0. 748 0. 653 0. 490 
— ‘ 0. 667 0. 640 0. 560 0. 427 0. 240 
2. 000 
F, and *F +. 000 
P 7 1. 732 1. 500 
F, and "f 3. 000 2. 250 
as ; 1. 604 1. 512 1. 105 
F, and ‘Ff 2 573 2 286 1. 429 
1. 447 1. 402 1. 254 0. 957 
F, and * 2.095 1. 964 1. 571 0. 917 
a - 1. 231 1. 206 1. 128 0. 985 0. 739 
F; and 'F, 1.515 1. 455 1. 273 0. 970 0. 546 
F. and 7 0. 905 0. 892 0. 853 0. 783 0. 674 0. 500 
— 0 818 0. 796 0. 727 0. 614 0. 455 0, 250 
H being the tield strength in oersteds. where 
If the repulsion in this case Is ¢, and the separation 
is A, eq 1 can be written a?—(4.669 X10 HY = 21.80 X10 Hl? em~*. (4) 


M?. in the following the 


Since 7? is a function 0 


é -a em", (3) 
A factor [2/4 will be set as f(.M3). In this notation 
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9 


will be 


eq 7) 
é a] V/ ») 
The values of function f being the same for M) 
as for (— M), we can write 
f(M3)— f \(—_ My? 6) 
Before we preceed to the computation of the 


effects for some Mn 1 lines it is of interest to estimate 
roughly the intervals between levels necessary to 
obtain measurable distortions, with a magnetic field 
of about 80,000 oersteds, which is close to the fields 
used for our spectrograms 

From eq 3 and 4 we can write 


A a 21.8010 SO,000 l7 em 
¢ ft ‘ 
7) 
The accuracy of the measured wave numbers 


evidently depends on the region of the spectrum. 
At about 25,000 em! (4,000 A) the accuracy, in the 
MIT spectrograms, is close to0.05 em This is 
the minimum value that e can have if measurable 
effects are to be observed We can then write 


17 [?~340] 
0.05 


em 

Equation 8 is a relation between the interval A 
and the interaction factor. If we select a transition 
for which J/* is approximately one unit (see table 3), 
observable effects will be possible with intervals A 
smaller than 340 cm™'. In in which the 
/*-values are greater than one unit, such as those in 
the °D and °F terms, much wider intervals will still 
show asymmetries. Because the term intervals in 
Mn rarely exceed 200 cm™', it is possible to find 
many asymmetries in it 

The foregoing estimate shows clearly why so many 
asymmetries have been noted in the past by ob- 
servers of the Zeeman effect. The asymmetrical 
patterns belong not only to the spectra of the light 
elements that have barrow term intervals, but also 
to elements of high atomic numbers whose term 
intervals are much wider. Ir these the 
[?-values are sufficiently great for the intervals to 
show measurable asymmetries 


causes 


cases 


4. Positions of Magnetic Levels Affected by 
Partial Paschen-Back Effect 


A spectroscopic level 7,, by influence of the mag- 
netic field, is split up into 2.J+ 1 magnetic levels 7¥, 
whose values are given by 

TY = T,+-M,aq,. 

We shall consider first the simplest case in which 
only one other level, 7,_,, exists close to the level T, 
of the same term. Suppose that the relative value of 
T, is smaller than that of 7, The magnetic levels 


of 7, will be lowered in value by the repulsions due to 
the magnetic levels of 7, The amounts of these 
displacements are given by eq 5 so that the final 
values of these perturbed magnetic levels 7’ will be 


7T’'*—T, M,aq, a*f{( MM?) (9 

The whole set of magnetic levels of 7, is lowered 
in value because the repulsions are all in the same 
direction; the magnitudes of these displacements are 
different for the individual levels on account of the 
different values that the function f(M3) has for 
different M,-values. Figure 2 illustrates the relative 
positions of the perturbed magnetic levels of y 
compared with the unperturbed positions. 

Levels with the same M-value have the same dis- 
placement independent of the sign of M1; thus, all 
levels M +O% have experienced 2 displacement a, 
and all with M=—+1% have been displaced by an 
amount 6. The direction of the displacements in 
T,, is contrary to that in 7, Levels M= +2' have 
not experienced displacement because there are no 
magnetic levels with such M-values in the level 7, 

The intervals between consecutive magnetic levels, 
such as those having magnetic quantum values M, 
and M,—1, may be computed by eq 9. The result- 
ing value is 

AT”}-¥ (M,—(M, 
a*f ((M, 


M 


L)jaq, a~fi 
1)’I, 
or 


AT’? ¥ a*f(M3)+-a*f\(M, 


ad, | 1)*] (10) 

The intervals are unequal because of the influence 
of the quadratic terms, which have different values 
for the different M-values. Figure 2 shows clearly 
this inequality. The pertrubed magnetic levels 
appear to approach each other in the opposite diree- 
tion of the perturbing level. 
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M, but is independent of the sign of 
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rturbations that arise from 7 ind Ty, may be considered to act independently By influence of the * 
{7 those of 7 are splacedupward by amounts a or ¢ hese partially perturbed levels are displaced 
lownward by amounts c, d, or ¢, by influence of magnetic levels of 7 


The separations between two magnetic levels hav- 


Ai’ ?-} ag, a? f{( M?)+ a? f( M2 ,,) 


ing the same M-value but different sign are 
” a® f\(M,—1)*]—a?f ((M,.,—1)*], (14) 
AT’+#¥.-¥=M,ag; M), ag,—a’*f( M3)+ aad 
a*f {((— M)5] AT’? } ag,+four quadratic terms 
But the quadratic terms being equal by eq 6 cancel, The separations between levels with the same 
1 Beg Jy &q I 
and thus we have absolute value of \/, but contrary sign, are 
AT’+¥.-#=2M,aqgz. (11) | AT’-}#- *+#=—M,agqg, M), ag,—a*f(M3)4 
an , ,; ; 20 M2. .)+a?f! 2)_@? : 
This final equation is independent of quadratic a*f( M3 1) +a°f\(— M)i)—a*f (— Mi 11, 
terms, and hence we conclude that the separation ' ' Ht 
between the perturbed levels is, in this particular | 2U' [5 equation by eq © will by 
case, independent of the perturbation and thus equal AT’-“+M—9 Mac (15) 
to the separation between unperturbed levels. Hence Fee 
we can write 
a athe a result that shows that the separations are inde- 
, f f ‘ o a 
AT’*} ATi ~™. (12) | pendent of the quadratic terms and therefore equal 


to those of the unperturbed levels 

The magnitude of the displacements given by (3) 
decreases when the J-values increase owing to the 
influence of J and A. Table 3 shows that in a term / 
decreases as J increases. In spectral terms that fol- 
low Landé’s interval rule, A increases with J; the 
function given by eq 5 will hence decrease when J 
In general this means that in a level 7, 
T are expected to be 
_ and hence after the 
appeal close to each 


We are now prepared to discuss the general case in 
which the affected level is perturbed simultaneously 
by two levels of the same term, one on each side of it. 
Figure 3 gives an illustration of this case. The level 
72, 18 perturbed by the two close levels 73,, and 7; 

In order to obtain the final perturbed levels, it is 
necessary to consider independently the influences 
of 7T,_, and T, Each of these influences will give 
a quadratic term to the equation, but as they are 
situated on opposite sides of io the signs will be 
different for the two quadratic terms. By the sim- 
ple case explained before we know that the displace- 
ments are always In a contrary sense to the perturb- 
ing level. The final values of the magnetic levels 
will then be 


increases 
the perturbations due to 7, 
greater than those due to 7, 
final perturbation the levels will 
other in the direction of 7,, as shown in figure 3 


5. Positions of Components in an Asymmetric 
Pattern 


The wave numbers of perturbed lines formed by 


{ 
; TY = T, + Msag,—a* f( M3) +a? f( M3 .,), (13) | transitions between perturbed magnetic levels 7”¥ 
and l’°'™) are evidently given by 
a* f( M35) being the perturbation due to 7,_, and a? 
(M? 1) that due to the level Zs n’™ r i? l | 


The intervals between two consecutive levels with 


quantum numbers ./, and M,—1 by eq 13 are but introducing here the level values given by eq 9, 
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we shall have 


7 Nyy V7: aq a~?t \/72) 


Ve. -M,ag 


a* {( M73) (16) 


This equation is a general formula giving all wave 


numbers in a perturbed pattern. In the case of par- 
allel components Mj—.M,. Equation 16 gives 


M;;) 


oe te M,adqg—a* f( M:; a*f ; (17) 


n 


In the case of the normal components, with AZ’, 
VU, 1 (or M, +1 without loss of generality 


‘MM 
ny 


V3)-+ 
(18) 


n, » +M,adAq 
a’® f((M, —1)*| 


ad 


Equations 17 and 18 are closely analogous to those 
for unperturbed patterns except for the quadratic 
terms that give the perturbations 

Intervals between adjacent parallel components 
may be obtained by subtracting from eq 17 another 
equation equal to it but in which M, and M,’ have 
been replaced, respectively, by M,—1 and M,’—1. 
The result is 


a f( MG )+ 
a*f\(M,—1)"|—e?*f[l 


adg a*f( M3) 


or briefly 
adg+four quadratic terms, (19) 

which shows that the intervals between parallel 
components are unequal because of the quadratic 
terms 

The intervals between adjacent normal com- 
ponents may be obtained by subtracting from eq 
18 another equation equal to it but in which M, 
and M,—1 are changed to M,—1 and M,—2 
respectively. The result will be 


adg tags ( ads) a’t(M, 
a® f |((M,, —1)*]+-a?f|(M, —1)"]—a?f ((M, —2/', 
or briefly 
adAg + 2aq,-+ four quadratic terms, (20) 


which shows that here also, in the normal components, 
the intervals are unequal owing to the influence of 
the quadratic terms 

The separations between components of the same 
order on both sides of the line without field, and 
having the same M-value but contrary sign, are 
given for the parallel components by 


a* f( M3) + a? f( M3.)-4 
My } 2M,adq. 


M),aAq 
M)3) 


M,adg { 


a’y \( a*t \( (21) 


normal | 
| 


The corresponding separations for the 


components are given by 


M),aAq 
M)3|—a?f | 


M,aAq+ag,—( a* f( M3)- 
a* f( M3,)+ a*f |( 


2M,adg- 


(@aqgy") 

M);,) 
2aq, (22 

The separations in these cases, either between 
parallel or normal components, are independent 
of the quadratic terms and hence equal to those 
found in patterns without distortions. This last 
conclusion is very important in calculating g-values, 
as we shall explain later. 


6. Asymmetric Positions of Components in 
Some Zeeman Patterns of Mn! 


We have predicted by theory the position of the 
components in an asymmetric Zeeman pattern. It is 
interesting now to compare the predicted positions 
with those found by measurement. For this purpose 
we have selected some lines of Mn that are strong 
and thus easily observable, and for which the 
predicted distortions are very great (see fig. 4). 

The selected lines belong to multiplets formed 
by combination of the low term a *D with the middle 
terms 2°D° and z°F°. The relative values of these 
terms and their intervals, A, are collected in table 4 
The theoretical displacements of these terms by 
partial Paschen-Back interaction have been com- 
puted by means of eq 3 and 5, by using the /*-values 
of table 3 and the value 3.955 cm~' for the Lorentz 
unit. The resulting values of the displacements 
are collected in table 4 under their respective 
M-values 

Table 4 shows that very great displacements are 
expected in levels with small J- and M-values, 
especially in the case of 2 °F°. 

Table 5 shows the calculation of the values of the 
perturbed magnetic levels belonging to the spectro- 
scopic levels of table 4. Only the levels with small 
J-values have been considered, because they show 
the greatest displacements. The computation has 
been made by use of eq 13; the values of the dis- 
placements of table 4 have been rounded to two 
decimals for this calculation. 

In the first column of table 5 are given the spectro- 
scopic levels, 7,, and their theoretical g-values. LS- 
coupling is supposed to be valid. The quantum 
numbers ./, of the magnetic levels and their respee- 
tive displacements ./,ag, are listed in the second and 
third columns. By addition of these displacements 
to the respective 7,-values the magnetic unperturbed 
levels T¥ are obtained as given in column 7. An 
LU 3.955 em™~' has been used in order to make the 
results comparable with the observed values. The 
values of the perturbed levels 7’¥ have been com- 
puted by adding the Paschen-Back perturbations, 
(cols. 4 and 5) taken from table 4, to the unperturbed 
levels TY. As these calculated perturbed levels do 
not differ much from the unperturbed ones, only 
their decimal parts are given (col. 8). 

Table 6 shows the lines formed by combining the 
magnetic perturbed levels of table 5. Under the 


resulting calculated wave numbers are given two 
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groups of decimals 


to the observ ed 


numbers 


These 


Those at the bottom correspond 
observed 


wave numbers were derived by measuring the wave- 
lengths of the lines of patterns on six MIT spectro- 


grams with dispersion about 0.8 A/mm 





The p 


ro- 





that described by Harrison and his collaborators [5] 


The magnetic field for two of the spectrograms was 


81,700 oersteds and for the others 84.700 


measurements have been reduced to the field 84,700 
and accordingly 


oersteds. 
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ments. Lines marked with an asterisk indicate 
poor measurements because of the close proximity of 
another line. 

With the observed 
observed levels was derived. The values of the 
observed magnetic levels are listed under their 
respective theoretical values in table 6. To check 
the accuracy of the measurements, these observed 
values have been used to recalculate the wave num- 
bers. The resulting values (only the decimal 
parts) are given between calculated and observed 
values in table 6. With very few exceptions the 
values agree to | or 2 units in the last decimal, 
thus indicating that the measurements are very 
consistent The close agreement between the ob- 
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served wave numbers and those derived for the 
observed magnetic levels is strong proof that the 
combination principle applies exactly to asymmetric 
patterns. 

The observed magnetic levels thus computed have 
also been listed for comparison purposes in table 5, 
where for simplicity, only the decimal parts of their 
values are given in column 9. The observed inter- 
vals between consecutive levels, listed in column 10 
are evidently unequal, as expected from theory. If 
the calculated and observed values for the per- 
turbed levels given in table 5 are compared with 
each other, it will be observed that in most cases the 
agreement is extremely good rarely exceeding the 
differences by more than 0.02 em in spite of the 
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great differences that exist between values for per- 
turbed and for unperturbed levels (compare cols 
7 and 8 of table 5 This proves that the approxi- 
mate theory of the Paschen-Back effect applies 
closely to the asymmetrical patterns of the spectrum 
Mini. There are, however, a few cases in which 
the differences between observed and calculated 
values, for the magnetic levels, exceed the expected 
errors of measurement. The most important cases 
are the following (see table 7). 

All the components having origins in these levels 
show, as may be seen in table 6, very great differ- 
ences between calculated and observed wave numbers. 
It is to be noted that all these components have origin 
in levels with magnetic quantum numbers 1/= +0! 
We shall consider this problem in the next paragraphs 


aR 7 Vaagnetii eve t/ ong perturbations in Mn I 
LI ) 
La 

P Ob i Obser | 

. slculated 

Des" $2192. 41 15 0. 04 
1) $2205. 59 55 0. 04 
I : 13641. Sl 7 0. 16 
I $3646. 03 5. 95 0. OS 
I ‘ 13675. 50 32 0. 18 
} 13672. 86 O4 0. OS 


7. Computation of g-Values From the Level 
Values 


In the preceding discussion we have seen that, 
with few exceptions, the calculated values are closely 


coincident with the observed, showing that the 
theory is valid for these cases of manganese. It may 


be recalled now that eq 15 of this theory, which gives 
the separation between pairs of magnetic levels with 
the same absolute /-value but different sign, is inde- 
pendent of the quadratic terms, and hence that equal 
separations are obtained for perturbed and_ for 
unperturbed levels. This means that it will be 
possible to eliminate distortion due to weak Paschen- 
Back interaction by using eq 15. Solving eq 15 we 
have 

ee alter 

Ms 2Mya 


which gives the g-value corresponding to each pair of 
magnetic levels T’~“% and T’*¥. This formula indi- 
cates that to obtain a g-value it is only necessary to 
divide the separation between the two levels in- 
volved by 2M,a. As an example the g-value corre- 
sponding to the observed magnetic levels 17435.18 


M 2%) and 17467.97 (Vf 2's), a being equal to 
3.955, may be calculated as 
7467.97 —17435.18 . 
17467 ‘ ‘ : 1.658. 
2 2 3.955 


In table 5, last column, the g-values computed in 
this way for each pair +.M are listed in front of the 
corresponding (—./)-values 


8. Computation of g-Values From the 
Observed Components of Patterns 


It is interesting now to compute g-values directly 
from the observed patterns and to compare the 
resulting values with those derived from levels that 
were given in table 5. Table 8 shows details of the 
calculation. In the first column are given the wave- 
lengths and transitions of the spectroscopic lines. In 
the second are listed the polarizations of the mag- 
netic components whose wave numbers are in column 
3. Lines marked with (*) are poorly measured 

In order to eliminate the distortions due to partial 
Paschen-Back effect the separations are taken, as 
indicated by eq 21 and 22 between components of 
the same order on both sides of the line without 
field These separations are written after the 
components of smaller wave number. Thus, the 
difference between the lines 24477.18 and 24492.84 


is written after the former The separations, 
which are given in cm“, are divided by 2a=2 >» 
3.955=7.910 to obtain the components of the pat- 


tern in Lorentz units. These resultant components 
are listed in column 5. As usual, parallel com- 
ponents are given in parentheses 

The next step is to calculate the mean values of 
the intervals between adjacent components of a 
pattern. The intervals are given in column 6, and 
at the bottom of each set of intervals is written their 
mean value. The last interval in each set is the 
double value of the smallest parallel component, 
which itself is an interval 

In table 8 the intervals in each pattern are so 
closely equal that accurate mean values have been 
obtained. Each normal component and the mean 
interval provide two equations for calculating the 
g-values of the two levels involved 
g-values are written after the respective normal 
component. Thus, the first component in table 8, 
namely 1.980, with the mean interval 0.2/2 gives 
the two g-values 1.874 for a°®D and 1.662 for 
2 DD, 

There are three lines 3843, 3833, and 3816 A whose 
intervals, within each pattern, are very different and 
hence have not been included in table 8. The line 
3816 A, for example, has the following components 
and intervals 


The resulting 


2. 544 |1 


Intervals 0. 58 


Components 987 1.334 0.775 (0, 885) (0. 326 


0.853 0. 559 0.559 0.6452 


The impossibility of deriving an accurate mean value, 
for the interval, with values such as 0.557 and 0.653 
is evident 

Recently, Catalan and \ elasco have shown. iS] and 
[9], that in cases such as these it is, nevertheless, 
possible to obtain g-values from the components if 
these are treated in separate groups, one for each 
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TABLE 8&8 lsymmetric patterns in 
I $4,700 oersteds; LU 3.9 
a*D1 z°D 
24477. 18 15. 66 1. 980 1. 874 1. 662 
; 24478. 32 13. 91 1759 | 0 22! 1. 865 1 653 
o 24479.21 | 12.28 | 1.552 | % <0 1. 870 1. 658 
: 24479.88 | 10.61 1.339 | % #19 1. 869 1. 658 
24483. 83 2. 49 0. 314 n 
9 g 208 
aD D 24484 88 0. 84 0. 106 ® ete 
9 O45 $40" : 
{082.0454 2 24485. 72 
a 24486. 32 Mean 0. 212 
° ; 24490. 49 
, 24491. 49 
24492. 23 
a 24402. 84 a*D DD 
0 24495, 58 20. 57 2. 600 - 3. 331 1. S60 
24501. 77 8.93 | 1.129 | 1-4/4 3. 322 1. 860 
“DD D x 24503 00 5. 75 0. 727 
079.415 oe 5 
1079.4 154A | - 24508. 75 M l a 4 
o 24510. 70 noe 1. 460 
24516, 15 a*Dp, D 
24548. 33 26. 31 3. 325 3. 3257 3. 325* 
~ co | 
aD D mee 13 0. 00 0. 000 
070.2804 cee” ae 
107 - 24561. 43* 
24574 64 
a*D Ly | 
J 24620. 42 20. 49 2 590 1. 868 3. 312 
g 24625. 87 9.06 | 1.145 | 144 1. 867 3. 311 
a*D ID r 24627. 84 5. 70 0. 721 > 
1058.936A 2/830. 08 1. 442 
- 24633. 54 
g 24634. 93 Mean | 1 
o 24640. 91 
1°D D 
7 24684. 54 15. 52 1. 962 p ox 1. 650 1. S58 
o 24685. 02 13. 90 1. 757 > ow 1. 653 1. 861 
o 24685. 79 12. 25 1. 549 U. SUS 1. 653 1. 861 
g 24686. 84 10.60 | 1.3941 |°™ 1. 652 1. 860 
x 24691, 02 2.48 0. 314 
a*tD ) 7 24691. 55 0. 84 0. 106 . 
1048.747A 24692. 05 + 
7 24692. 39 0. 212 
24603 50 Mean 0. 208 
@ 24697. 44 
a 24698. 04 
. 24698, 92 
e 24700. 06 
a*D F, 
) 26018. 10 16. 97 2. 145 Oe 1. 867 1. 311 ' 
s 26020.62 12.53 1.584 261 1. 862 1. 306 
, | 2602293 | 820 | 1.037 | 724 1. 871 1. 315 
. oe men ; Ht) = 
a 26025. 18 3. 74 0. 477 1. S67 1. 311 
x 26023. 46 6. 59 0. 833 — 
a*D ; " 26025. 75 223 |(0, 282) | % 2?! 
3841. O74A 8027. 03 tat 
" 26027. 98 = 
" 26030. 05 Mean | 0. 566 
o 26028. 05 { 
; 26031. 13 
o 26033. 15 
e 26035. 07 
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TABLE 8 isymmetric patte 


n NIn I 


1. 369 


1. 963 


3. 105 
0. 631 
1. 236 


Mean 


2 190 
1. 485 
2. 149 
0. 791 
0. S58 


0. 167 


Meat 


Continued 


0. 342" 
0. 3453* 
3a 
0. 358 
0 * 
0. 3L16* 
0. 334 
O. SAF 


D 
S6S 
SO7 


671* 


658 
HOS 


H56* 


504 


606 
605 


effect; but it must be applied with caution because 


m 


S prec ! Polariza Or Separat 
’ wa ! I pa I 
| a 26031. 20 10. 83 
a*D I 7 26028. 83 15. 53 
3839.779A 26035. 81 
” 26044. 36 
o 26042. 03 
| a 26003. 27 24. 56 
o 26103. 30 1 99 
a*D I } a 26100. 92 4. 78 
3820. 670A } PAIOL. 47 
7 26110. 70 
o 26108. 29 
\ o 26117. 83 
a 26135. 17 17. 32* 
a 26136. 47 
o 26137. 89 11. 75 
J 26139. 36 9. OU 
og 26141. 15* 6. 26* 
n 26140. 70 6. 79* 
© 26141. 80 
a*D I } r 26143. 05 1. 32* 
S823. 801A 26143. 98 
- 26144. 37* 
7 26145 
7 26147. 49* 
| a 26147. 41 
| a 26148. 45 
| a 26149, 64 
o 26150 
’ 26152. 49* 
pair of magnetic levels Vf. Table 9 shows the 
calculation of g-values by that method for the three 
lines mentioned above. The general arrangement of 


table 9 is similar to that of table 8, but the com- 
ponents have been placed in two groups. It has 
been possible to obtain the g-values for the two 
levels involved in all cases except one, that of the 
component —0.064, because in this case the interval 
is unknown. The g-value of one of the levels in- 
volved, namely 3.328 for a°Do., has been adopted, 
and then the g-value for the other level has been 
derived Table 9 shows that the g-values obtained 
for a®°Da»,., a®Dyy, and a*D qa, are closely coincident 
with those obtained in table 8, thereby justifying 
the method of Catalin and Velasco. 

It is to be noted that the g-values obtained for the 
same spectroscopic level, in different patterns, are 
practically identical. Very accurate mean values 
may be derived from the data of tables 7 and 8. 
The resulting values have been compiled in table 10 

The agreement between values derived from leveis 
or from components is extremely good. This is a 
proof of the validity of both methods of deriving 
g-values. The method of levels, although a little 
more difficult to apply, is the more general, for it 
gives results in all cases. The method of the com- 
ponents is, however, much more simple and is the 
one usually emploved by all observers of the Zeeman 


curate g-values 


may serve 
values 
Most of the levels in table 10 have g-values equal 


to 


in 


some asymmetric 


in such 


cases 


patterns 
The method of Catalan and Velasco 


those of LS-coupling 


table 7 


constitute 
ese exceptions are just those derived from levels 
The two 
*Dse* have been displaced unequally 


clear 


magnetic 


it 


may give 


obtain more 


“Dan, 2 °F: 
exceptions to the rule 


less ac- 


reliable 


and 


levels 2*Ds°" and 


The level 


with negative \/-value has received a positive ad- 


ditional impulse of 0.04 em 


positive value has received 


push of 0.04 em 
levels compared with the distance between the un- 
perturbed levels is shortened by 0.08 em 


the corresponding g-factor Is altered 


and that 
additional 
The distance between these two 


The 


with the 
negative 


Hence 


q-value 


for LS-coupling is 3.333, and the real value is 3.312 
the 


We 


cannot 


suggest 


displacement 


other 


two pairs M 


The other cases are closely 
as Catalin and Velasco have shown 


p 200] 


different 
but the pair .V/ 0 
smallerthan the LS-value 
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for this anomalous 


connected with each 


[9] and 


The g-values that are obtained for the 


I's and M 
The first pair gives exactly the LS-value 
gives 1.006, which is 0.061 LU 
The additional impulses 


0 


for z °F, 


are quite 








PARLE 0 Zeemar patterns of lew 6°. of Mn 








Magnetic field 84,700 oersteds; LU 3.955 er 
a ie : , _— . separa — 
POCeTUSCOpat _ ; . . . _— thor — Inter i Vaile 
hand wavelengt! me t 
a*Dp ‘ I ‘ ° 
o 26007. 65 0. 51 0. 064 [3. 328] 1. 067 
HOO}. 30 
, 26007. 14 
| Fag 
‘ ‘ 7 25997. 85 17. 12 2. 165 3. 328 1. 002 
"Dog Pus ffs 6001, 82 9, 20 1. 163 
3843. OSSA - os 
‘BOD, if 
7 26011. 02 
7 26014. 97 
a*Diy | a 
0 26074. 00 5. 26 0. 665 1. S65 1. O65 
- 26071. 57 9. 49 1. 200 
"O75. 97 
cd 26081. 06 
o 26079. 26 
Fie 
‘1D 6} 0 26066. 36 18. OS 2. 285 0&5 1. 857 1. 000 
a . . « ‘ , a4 ‘ 
323% 865A" o 26069. 96* 11. 32* 1. 431 1. 860 1. 003 
a r 26073. 94 3.40 (0. 430 
> » Oo S60 
98OP A. GDS 
- 20077. 34 Mean 0. 85? 
; 26081. 28* 
o 26084. 44 
a D., | “ bg 
o 26183. 29 20. 12 2. 544 1 789 1. 659 1. O69 
o 26190. 42 6. 13 0. 775 1. 660 1. 070 
© 26189. 94 7. 00 0. 855 , 
» ] 0 
AID? 98 
7 26196. 94 Meat 1.770 
oO 26196. 55 3 Mear 0. 590 
6F =0% 
a*D.. Fis : @ 26203. 41 
3816. 746A a 26184. 65 15. 72 1. O87 O Bae 1. 660 1. OO7 
o 26187. 25 10. 55 1. 334 = 1. 660 1. 007 
, 26191. 24 2. 58 0. 326 . 
Oh 4A? 
Yt 192? 93 = 
- 26193. 82 Mean 0. 653 
o 26197. 80 
o 26200. 37 
] 
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The 


respectiv ely 


amount of 0.024 em 
of 0.061 LU in the g-values 
levels 6F > 


received extra displacements 


magnetic 


and 


> 8+ 
0.18 and 0.08 em 


0 
0% 


, which accounts for the change 


have 


In this particular case of *Fo., on ac- 


count of the fact that the g-value is negative, the 


TaRLE 10. Mean g-values for levels of Mut same. The sum in LS-coupling amounts to 1,067 
- 1 84,700 oersteds} 0.667)=0.400 and to 1.001 0.602) =0.399 
for the perturbed levels A clear explanation of 
Ke computed with this case has been given by these guthors (7, p 200) 
ne Scoupli 9. Theoretical Computation of Line Strengths 
7 —— in Zeeman Patterns 
a *Dorg 1. 658 1. 658 1. 657 The strength, S, of a line is defined as a quantity 
a *Diy 1. 867 I 866 1. 867 that must be multiplied by the fourth power of the 
; - ‘ ; a 1 pee 1 an? frequency and by the number of atoms in any one 
D . , REE 1 861 , 867 of the initial states, in order to obtain the radiated 
: energy \7, p 200}. There is a possibility of com- 
1) . 312 3. 312 ». 335 paring relative strengths with observed intensities be- 
* Vv 1! 067 068 i cause all the lines of a Zeeman pattern have closely 
6} VV 0 , 006 1. 004 1. 067 the same frequency and are produced by close states 
F 0. 602 0. 602 0. 667 that have practically the same excitation 
The strength S” of a component in a Zeeman 
pattern depends on the strength S of the parent 
0.16 and —0.08em~', which 2*F-° and 2 *F*°: have | line through the equation 
received, make their distance shorter by a total 


§* (24 


SK, ( 
in which A is a constant whose values depend on 
the J- and .-values of the initial levels as expressed 


in the following formulas 


Magnet 








. Spectroscopic on 
magnetic level M 0’ is smaller than M= -+0%. transiti tra K-valu 
Hence the distance between these magnetic levels 
has been shortened by a total amount of 0.26 em M—M 2.4 (J?— M2) /J(2J—1 
and this accounts for the change of g-value from J—J—1 ; 2d + I 
0.667 (LS-coupling value) to —0.602. eres ) as hast I ‘ili 
There is a close connection between the changes in — VM QAMJ(d +1) (24-41 
g-values experienced by 2°F},, and z°Fo,. As Catalan VM +1 1/2) AJ = M) (J 4 M1 
and Velasco have shown, the sums of the g-values for JJ +1) (2541 
LS-coupiing and for the observed levels are the 
PARLE 1] K- and (kK -values for components in a Zeeman patter? 
Upper figures are A-values; lower (A values 
Odd multiplicities, parallel 
V 3 2 l 0 l 2 3 
J J 
V 3 2 | 0 2 3 
50. 00 
I 0 7. 07 
a ' 00 20. 00 15. 00 
= 87 1.47 3. 87 
3 ° 7.14 13 12. 86 11. 43 7.14 
; 4 2. 67 38 3. 58 3. 38 2. 67 
' ; 1. 16 7.14 93 9, 52 8. 93 7.14 1. 16 
| 2. 04 2. 67 99 3. 09 2. 99 2. 67 2.04 
J J 
, 00 0. 00 25. 00 
00 0. 00 5. 00 
e - 20. 00 00 0. 00 5. 00 20. 00 
1 2 1.47 24 0. 00 2. 24 1. 47 
3 3 16. 07 7. 14 79 0. 00 1. 79 7. 14 16, 07 
: ; 1. O01 2. 67 29 0. 00 1. 29 2. 67 1.01 











PaBLe 1! K- and (K)"*-values for components in a Zeeman pattern—Continued 
 pper figures are A-values; lower (K value 
Odd MuUultipheitic normai component 
V 2 l 0 l 2 $ } 
j J] 
V1 3 2 l 0 l 2 3 
25. 00 
{ ua 
, 5. 00 
. 2. 50 7. 50 15. 00 
” 1. 58 2. 74 3. 87 
, . 0. 71 2.14 4. 29 7. 14 10. 71 
” 0. 84 1. 46 2.07 2. 67 3. 27 
' ; 0. 30 0. 90 1.78 2. 98 1. 49 6. 25 S. 33 
, 0. 55 0. 95 1. 33 1.73 2. 12 2. 50 2. 89 
j J 
; 12. 50 12. 50 
3. 54 3. 54 : 
. 9 5. 00 7. 50 7.50 5. 00 
. 7 2. 24 2. 74 2.74 2. 24 
, , 2 68 1 46 5. 35 5. 35 1. 46 2 68 : 
7] ‘ 
1. 64 2. 11 2. 31 2. 31 2. 11 1. 64 : 
; 
I n multiplicities, parallel 
V ; 2 | 0 0 | 2 3 
j i | 
V ; 2 l 0 0 | 2 3 
0 25. 00 25. 00 
5. 00 5. 00 
. 10. 00 15. 00 15. 0O 10. 00 
r 3. 16 3. 87 3. 87 3. 16 
. . 5. 36 8. 93 10. 71 10. 71 8. 93 5. 36 
‘ , 2. 32 2. 99 3. 27 3. 27 2. 99 2. 32 
{ . 3. 33 » 83 7. 50 8. 33 8. 33 7. 50 5. 83 3. 33 
; 1. 82 2. 41 2. 74 2. 89 2.89 2.74 2. 41 1. 82 
j j 
25. 00 25. 00 
{ 
0 , 5. 00 5. 00 
; 22. 50 2. 50 2. 50 22. 50 
1. 74 1. 5S 1. 58 1. 74 
. a 17. 85 6. 43 0. 71 0. 71 6. 43 17.8 
= - 4. 22 2 54 0. 84 0. 84 2. 54 4. 22 
3 . 14. 60 7. 45 2. 68 0. 30 0. 30 2. 6S 7. 45 14. 60 
: : 3. 82 2.73 1. 61 0. 55 0. 55 1. 61 2.73 3. 82 
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TABLE 11 K- and (K 


Upper figures are 


-values for 


The value J is the larger of two J-values involved 
For these equations it has been assumed that the 
observations have been made perpendicularly to the 
magnetic field. A-values have been computed by 
means of these equations and are given in table 11 
The value 75 has been assigned to constant A so 
that the sum of all parallel components in a pattern 
may be 50 and that the corresponding sum for each 
of the two groups ol normal components may be 25 
Thus the sum of all components in each pattern will 
be 100 

In calculations related to the strength of lines the 
root (§ of S is more convenient than the 
strength itself for comparison with the observed 
intensities. Accordingly, we have computed the 
square roots of the A-values. In table 11 the square 
given under their S-values 
Taking the square root of both members of eq 24, we 
have 


square 


roots are respective 


S™ Ss K 


This equation means that to obtain (S“)*-values 
it is only necessary to multiply the values of A” by 
the strengths of then respective spectroscopic lines 
The quantities (S” S*, and AK” are either positiv: 
or negative. The sign of (S”)* depends on the 
sign of S* and K The sign of S* depends on the 
changes of the L- and J-values in the spectral line 
The sign of K”® depends on the changes in the J- 
and M-values for the magnetic levels of the com- 
ponents Tables 12 and 13 give the signs for the 
different cases 
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com pom nis 


K-values wer 


Continues 


PARLI 


VJ 


of Com- 


oF° 


10. Theoretical (S™)” 


Strengths 
ponents in Patterns of a Multiplet °D 


The strengths mav be computed bv formula (25 
The S-values for lines'in a multiplet "D—*F 
deduced by the classical formulas of Kronig, of 
Sommerfeld and Héal, and of Russell [10] The 
square roots of the values obtained are given in table 
14. The signs have been computed by table 12 

In order to compute the final S™”)*-values we have 
multiplied, for each pattern, the strengths of its com- 
ponents, given in table 11, by the strength of the line 


can be 








ARLE 14 Va S ) F° bin 
; D 
I 18. 07 
} 7. 45 15. 63 
| 211 9, 02 » 42 
I 3. 21 9.13 , 30 
I 3. 79 8% 26 §. 11 
I § 65 6. 83 
without field, given in table 14 The resulting 


values are given in table 15 


11. Theoretical Computation of Line Strengths 
in a Zeeman Pattern Perturbed by Partial 
Paschen-Back Effect 


The partial Paschen-Back effect, in LS-coupling, 





produces intensity distortions in the component 
lines of the patterns These distortions may be 
calculated by the 30-called second-order pertur- 
bation theory In this theory two magnetic levels 
PaBLe 15 Ca i ngth n multiplet a *I 
Magnet 47M Uy iy are aluc 
VJ 
> 38. 5 17. 3 
7 37. 1 17. 4 
17.3 23. 2 23. 7 
17.8 20. 6 24. 1 
23. 7 7.7 23. 2 
0 4 0 1 6 23. 2 
F} 
0 23. 2 2 23. 7 
23. 1 10. 7 23. 3 
23. 7 23. 2 
, 23. 2 25. 8 
) 17.3 
™ 16. 6 
18. 4 12.0 12 
11.8 2 4 >». 6 
0 10. 4 l4 7 7.4 
r 8. 3 15. 0 9.0 
0 7.4 14.7 10. 4 
5. 7 14.3 12.4 
| 1.2 12.0 
2. 9 11.0 
le 
k3 
ih 
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T¥ and T*_, belonging to consecutive spectroscopi 
levels 7, and 7,_, of a term TJ, and having th 
same \/-value, interact with each other, perturbing 
mutually the strengths of their combinations wit! 
another third level If T¥<T¥ the 
third level being inappreciably lower or higher thar 
the other two the perturbed strengths of the 
transitions are given by eq 26a and b, which we 


we set 


take from Kiess and Shortley [7, p. 204] 
pert .(S}" unpert.($} 
(al /A) unpert.(S¥ , 6a 
pert.(S} unpert.(S} 
(al /A) unpert.(S} (26b 
The first members of these equations are the 


strengths of the perturbed transitions. The quantity 
a is the LU in em™! for the magnetic field, J is th 
interaction factor given in table 3, and A the distance, 


in em between the two interacting magneti 
levels. The remaining factors are the strengths of 
) k f Mnit th and without perturbation 
with t Lows with perturbat 
D ) 
32.9 
33. 5 
20.4 25. 5 
32. 4 24.3 
10. 4 36. 0 18. O 
12. 3 37. 2 15. 9 
18.0 36. 0 10.5 
20. 2 34.8 8 5 
17. 3 25. 5 20.4 
16. 7 26. 7 26. 4 
38. 5 32. 4 
10. 0 32. 4 
34. 2 22. 6 26. 5 
36. 7 23. 4 27. 1 
22. 6 13. 1 26. 1 30. 6 15. 3 
24.3 7.8 25. 8 36. 7 10. 1 
26. 1 13.0 22. 6 15.3 30.6 
26. 6 18. 3 20. 9 21.5 24.4 
13. 4 22. 6 39. 2 26. 5 
15. 1 21.9 $1. 6 24. 8 
15.8 18. 3 9. 1 34.2) 34.2 
11. 1 18. 5 12. 1 29. 1 38. 0 
9 1 18. 3 15.8 34.2 34.2 
6. 2 18. 0 20.5 30.3 39. 2 


























! U raBLe 16 Va J /A for + 
it LI 
if 
i 
I Ml 
ie 
J D 1°*D 169, 52 1. 050 0. 959 
; I) a*D 116. 06 1. 200 0. USO 
2) 1°D 68. 67 1. 247 
I 71. 42 1. 355 1, 258 
I I 18. 95 1. 470 1. 200 
. ; } 8. 51 1. 490 
i the unperturbed transitions For a component 
) : perturbed by two others the changes are additive 
Special care must be given in these calculations to 
t ' the signs of the strengths 
‘ 12. Perturbations in the Strengths of the 
; Components in Some Patterns of Mn! 
i 
vf ; We are now prepat “l to compare the observed 
: intensities in some patterns of Mn1 with the theo- 
i retical strengths computed by eq 26a and b. For 
this comparison we have selected the multiplet 
a*D—z°F° for two reasons. First, this multiplet 


has an outstanding intensity, making it easy to 
observe all the components of the patterns, even 
the faintest Second, the J factors for such terms 
have high values (see table 3 especially for the 
lower J-values of z2°F°, and the separations A 
are relatively small, so that the magnitudes of J/A 
are very important 

First we have computed the values of a//A for 
the different transitions in this multiplet. Only the 
values corresponding to small J-values have been 
considered. Within the accuracy of this calculation 
the separation betwee n spectroscopic levels ma be 
used instead of that between magnetic levels belong- 
ing to the spectroscopic levels. Values of separations 


j D F° of Mnt 

Vw 

0. 742 0. O25 0. 022 0. O17 
0. O41 0. 033 
0. O72 

0. O58 0. O75 0. O70 0. 053 
0.119 0. 007 
0. 207 

the levels of a®°D and F A few examples will 


illustrate in detail how these calculations have 
been made 
We select first a very simple case, that of the com- 


ponent a°D$ FS Its initial state a°®D®* is 
perturbed only by a°D?%, and its final state 2 °F® 
only by 2 ; The strengths of the unperturbed 


transitions involved, taken from table 15 (in the table 
the values are rounded) are 


a ®Dos—»2 °F 30.55 

a°*D$ 2 OR 18.25 a*D? >2 *ke 34.15 

Factors a//A taken from table 16 are 0.207 for 

the influence of 2°F°®, 2°F8®, and 0.072 for that of 

a®*D§e, a®Dke. The values of the perturbations will 
then be 

0.207 * 30.55=—6.32 and 0.072 18.25 1.31 


The first is to be added eq 26a) because it is due to 
the influence of level 6F on another higher level 
: The second perturbation also is to be 
adda because it 1s due to the influence of level 
1°*D,. on another higher one, a*D Hence the 
final value of the perturbed strength will be 34.15 
6.32—1.51 39.16 


have been taken from table 4 and J-values from The line a*D?3—>2°F¢S is a more complicated ex- 
table 3. For the LU, a=3.955 has been used. The | ample. Its initial state is affected by both a *D$ 
resulting values are given in table 16 and a°D§ Its final state is perturbed by 2°F8 
We shall proceed now to compute independently and 6Fe The — — pe rturbation of 
the strength perturbations under the influence of | the components (see table are 
a *—° $F? 35.99 
a *D' »2 °F 14.67 a °° oFe 13.05 a *D? ‘Fe? 30.55 
a’ pD° Fe 18.25 
The values of aJ/A, as given by table 13 are 0.04) | will be = 13.05 0.60) —2.20+-4.38 3.78 


and 0.072 for the influences of Da,, Dy, and D, 
D,.., respectively; and 0.119 and 0.207 for those of 
Puc, Fry and Fo. Fyy, respectively. 

The values of the perturbations are thus 0.041 
14.67 0.60, 0.072 (30.55 2.21 0.119 
35.99) =4.28, and 0.207 > 18.25)=3.78. These 
perturbations will be added for Dy, Dy, and for 
Fy... Fy; and subtracted for Day, Dy, and for Fo., F, 
because of the relations Dy. < Dy.< Dog and Fac< 
Fic<Fo.. The final strength of the perturbed line 


18.31 Thus the line with strength 13.05 has been 
distorted, and its strength has been changed by 
18.31. 

With the same procedure we have computed the 
perturbations for some other patterne of multiplet 
1°*D—z*F°® of Mni. The resulting values of this 
calculation are given in table 15. A comparison 
that we have made of the observed intensities with 
the values given in table 15 shows that in all cases 
the intensity perturbations are in the right direction 
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and of the right order of magnitude to account for 
the observed asvmmetries We do not intend to 
give here our own intensities, because they have 
little quantitative value, representing only visual 
estimates. However, it must be noted that par- 
ticular care was taken in these estimates, to note in 
unperturbed com- 


each pair of two symmetrical 
perturbed 


ponents which was the stronger 
Readers of this paper may check for themselves the 
accuracy of these results by making visual estimates 
of the component intensities in figure 4 and com- 
paring their resulting values with those given near 
which represent the calculated 


when 


the components, 
perturbed strengths. 
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